Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 673,612 T→C D391G (GAC→GGC)  leuS ← leucine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913673,6120TC100.0% 112.2 / NA 33D391G (GAC→GGC) leuSleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/20);  total (13/20)

ACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGTCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCAGAGCCGTCAGCTGCCAGGATAACCGGTTTGATGTTCAGGCCGTAT  >  NC_000913/673480‑673731
                                                                                                                                    |                                                                                                                       
aCCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGAcc                                                                                                                  >  1:296570/1‑140 (MQ=255)
aCCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGAcc                                                                                                                  >  1:463008/1‑140 (MQ=255)
        aaTCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACt                                                                                                          <  2:550605/140‑1 (MQ=255)
           cGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGc                                                                                                       <  2:15264/140‑1 (MQ=255)
             gcgcgcCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCa                                                                                                     <  1:402490/140‑1 (MQ=255)
              cgcgcCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCag                                                                                                    >  1:553772/1‑140 (MQ=255)
                      aGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACCCCCATCGCCGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGTCAAGACCGTTGAACTCGCCAGAGTTGGAc                                                                                            >  2:312727/1‑140 (MQ=255)
                            gCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACg                                                                                      <  2:116527/140‑1 (MQ=255)
                             cTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGc                                                                                     >  1:254597/1‑140 (MQ=255)
                             cTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGc                                                                                     <  2:183699/140‑1 (MQ=255)
                             cTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGc                                                                                     >  1:656072/1‑140 (MQ=255)
                             cTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGc                                                                                     >  1:142670/1‑140 (MQ=255)
                             cTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGc                                                                                     >  1:394502/1‑140 (MQ=255)
                                  gggAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCtttt                                                                                <  2:553772/140‑1 (MQ=255)
                                          ccccAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTtcagtcag                                                                        >  2:360467/1‑140 (MQ=255)
                                             cAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGc                                                                     <  2:254597/140‑1 (MQ=255)
                                              aGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCt                                                                    >  2:544552/1‑140 (MQ=255)
                                                 cgcgCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGc                                                                 <  2:656906/140‑1 (MQ=255)
                                                       ggCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGaa                                                           <  1:360467/140‑1 (MQ=255)
                                                                     ttACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCaga                                             >  2:377111/1‑140 (MQ=255)
                                                                         gCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGACTCAGAGCCg                                         >  2:448117/1‑140 (MQ=255)
                                                                             gCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCAGAGCCGTCAg                                     <  1:337267/140‑1 (MQ=255)
                                                                                       aTCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCAGAGCCGTCAGCTGCCAGGAt                           <  2:296570/140‑1 (MQ=255)
                                                                                              tcagtTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCAGAGCCGTCa                                      <  1:638747/122‑1 (MQ=255)
                                                                                              tcagtTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCAGAGCCGTCa                                      >  2:638747/1‑122 (MQ=255)
                                                                                              tcagtTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCAGAGCCGTCAGCTGCCAGGATAACCGGt                    <  2:394502/140‑1 (MQ=255)
                                                                                                agtTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCAGAGCCGTCAGCTGCCAGGATAACCGGttt                  <  2:463008/140‑1 (MQ=255)
                                                                                                 gtTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCAGAGCCGTCAGCTGCCAGGATAACCGGTTTg                 <  1:309109/140‑1 (MQ=255)
                                                                                                      tCGGCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCAGAGCCGTCAGCTGCCAGGATAACCGGTTTGATGtt            >  1:411597/1‑140 (MQ=255)
                                                                                                        ggCGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCAGAGCCGTCAGCTGCCAGGATAACCGGTTTGATGTTCa          <  1:448117/140‑1 (MQ=255)
                                                                                                          cGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCAGAGCCGTCAGCTGCCAGGATAACCGGTTTGATGTTCAgg        <  2:411597/140‑1 (MQ=255)
                                                                                                          cGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCAGAGCCGTCAGCTGCCAGGATAACCGGTTTGATGTTCAgg        <  1:664796/140‑1 (MQ=255)
                                                                                                          cGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCAGAGCCGTCAGCTGCCAGGATAACCGGTTTGATGTTCAgg        <  1:544552/140‑1 (MQ=255)
                                                                                                          cGATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCAGAGCCGTCAGCTGCCAGGATAACCGGTTTGATGTTCAgg        >  2:582975/1‑140 (MQ=255)
                                                                                                           gATGGCGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCAGAGCCGTCAGCTGCCAGGATAACCGGTTTGATGTTCAGGc       <  2:656072/140‑1 (MQ=255)
                                                                                                                cGTTGAAGGCCGCTTCATGGCCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCAGAGCCGTCAGCTGCCAGGATAACCGGTTTGATGTTCAGGCCGTAt  >  1:469605/1‑140 (MQ=255)
                                                                                                                                    |                                                                                                                       
ACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGTCAAGACCGTTGAACTCGCCAGAGTTGAACAGCACGCCTTTTTCAGTCAGGGCTTGCTGAGAAAGATCTGGCTCAGAGCCGTCAGCTGCCAGGATAACCGGTTTGATGTTCAGGCCGTAT  >  NC_000913/673480‑673731

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