Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007779 3,542,911 G→A G168D (GGC→GAC)  rhaD → rhamnulose‑1‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077793,542,9110GA100.0% 152.6 / NA 38G168D (GGC→GAC) rhaDrhamnulose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/21);  total (17/21)

GATCGGGTGATCATGCACTGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGGCAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATCGGCCAGGCGACCGCACAAGAG  >  NC_007779/3542814‑3543008
                                                                                                 |                                                                                                 
gATCGGGTGATCATGCACTGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGacag                                                                                                <  1:1401471/101‑1 (MQ=255)
 aTCGGGTGATCATGCACTGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGacagc                                                                                               <  1:1178192/101‑1 (MQ=255)
       tGATCATGCACTGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAg                                                                                         >  3:250630/1‑101 (MQ=255)
              gCACTGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCtgg                                                                                  >  1:195305/1‑101 (MQ=255)
                aCTGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCtggtg                                                                                >  3:514338/1‑101 (MQ=255)
                 cTGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCtggtgg                                                                               >  3:714505/1‑101 (MQ=255)
                  tGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTggtggt                                                                              <  1:1345748/101‑1 (MQ=255)
                     cacgccacCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTAtt                                                                           >  1:1089312/1‑101 (MQ=255)
                      acgccacCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTc                                                                          <  1:1509199/101‑1 (MQ=255)
                            ccaACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGAt                                                                    <  3:1427942/101‑1 (MQ=255)
                             caACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATg                                                                   <  3:1173373/101‑1 (MQ=255)
                                 cTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGt                                                               >  3:865712/1‑101 (MQ=255)
                                  tGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGtt                                                              >  1:62605/1‑101 (MQ=255)
                                    aTCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTgg                                                            <  1:1204530/101‑1 (MQ=255)
                                       gCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCAt                                                         >  3:816521/1‑101 (MQ=255)
                                        cccTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCAtt                                                        <  3:1426617/101‑1 (MQ=255)
                                          cTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCAtttt                                                      <  1:1265785/101‑1 (MQ=255)
                                           tCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTg                                                     <  3:499427/101‑1 (MQ=255)
                                            cACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGc                                                    <  3:710642/101‑1 (MQ=255)
                                               cTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGt                                                 <  3:871105/101‑1 (MQ=255)
                                                  tGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGa                                              <  1:473629/101‑1 (MQ=255)
                                                  tGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGa                                              >  3:722042/1‑101 (MQ=255)
                                                   gTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGAt                                             <  1:1390111/101‑1 (MQ=255)
                                                     aCTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGAGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATgg                                           >  3:1319913/1‑101 (MQ=255)
                                                      cTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGt                                          <  3:689276/101‑1 (MQ=255)
                                                            aaCGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCgg                                    >  1:753483/1‑101 (MQ=255)
                                                              cGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCa                                  <  1:1017232/101‑1 (MQ=255)
                                                                  acCGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGa                              >  3:590071/1‑101 (MQ=255)
                                                                   ccGCGGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGAc                             <  3:1009417/101‑1 (MQ=255)
                                                                    cgcgGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACg                            <  3:175857/101‑1 (MQ=255)
                                                                     gcgGTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGa                           >  1:275319/1‑101 (MQ=255)
                                                                       ggTCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTACCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGaaa                         >  3:175705/1‑101 (MQ=255)
                                                                        gTCTTCACTCGCCAACTGTGGGAAGACAGAACCGGGTGTCTGGTGGAATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGCGCCCGGCACGGACGAAAt                        >  1:179771/1‑101 (MQ=255)
                                                                         tCTTCACTCGCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATc                       <  3:1407015/101‑1 (MQ=255)
                                                                                  gCCAACTGTGGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATCGGCCAGGCg              >  1:1130675/1‑101 (MQ=255)
                                                                                        tgtgGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATCGGCCAGGCGACCGca        <  3:63354/101‑1 (MQ=255)
                                                                                         gtgGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATCGGCCAGGCGACCGcac       <  1:580159/101‑1 (MQ=255)
                                                                                         gtgGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATCGGCCAGGCGACCGcac       >  3:245761/1‑101 (MQ=255)
                                                                                          tgGGAAGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATCGGCCAGGCGACCGcaca      <  1:102955/101‑1 (MQ=255)
                                                                                              aaGACAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATCGGCCAGGCGACCGCACAagag  <  3:1494787/101‑1 (MQ=255)
                                                                                                 |                                                                                                 
GATCGGGTGATCATGCACTGCCACGCCACCAACCTGATCGCCCTCACCTATGTACTTGAAAACGACACCGCGGTCTTCACTCGCCAACTGTGGGAAGGCAGCACCGAGTGTCTGGTGGTATTCCCGGATGGCGTTGGCATTTTGCCGTGGATGGTGCCCGGCACGGACGAAATCGGCCAGGCGACCGCACAAGAG  >  NC_007779/3542814‑3543008

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: