Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 97,176 G→C 8.5% V30V (GTG→GTC)  murD → UDP‑N‑acetylmuramoyl‑L‑alanine‑‑D‑glutamate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091397,1760GC 8.5% 104.0 / 2.5 42V30V (GTG→GTC) murDUDP‑N‑acetylmuramoyl‑L‑alanine‑‑D‑glutamate ligase
Reads supporting (aligned to +/- strand):  ref base G (14/24);  new base C (2/2);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

TTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGA‑AGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCGCCAGTCCCGGTA  >  NC_000913/97035‑97309
                                                                                                                                               |                                                                                                                                      
tttCGCTGAGGCTGGTTCTGATTAGTCTGGCGTGGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGa                                                                                                                                       <  2:743987/143‑1 (MQ=255)
       gATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCcgcg                                                                                                                                <  1:511692/143‑1 (MQ=255)
        ctgttGGTGCTGATTGGTCTGGCAACGCTGAAGGTCCGTTAATCATGGCTGATTATCAGGGTAAAACTGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGC‑GTCACACcgcga                                                                                                                               <  2:663388/139‑2 (MQ=255)
          gcgggTTCTGAGTGGTCTGGCGACGGTGAAGGTTCGTTAATCCTGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCGCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCcgcgagc                                                                                                                             <  2:662561/140‑4 (MQ=255)
               ttCTGATGGGTCGGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCGTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATgg                                                                                                                        <  2:457722/143‑1 (MQ=255)
                   gtttGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATAc                                                                                                                    <  1:990336/141‑1 (MQ=255)
                   gATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATAc                                                                                                                    >  1:374085/1‑143 (MQ=255)
                      tGGTCTGGCCACGCTGGAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTGTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATAcgcg                                                                                                                 <  2:1044083/143‑1 (MQ=255)
                       ggTCTGGCAGCGCTGAAGGTACGTTAATCAGGGCTGATTATCAGGGTAAAAATGTCGTCATTATGGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGt                                                                                                                <  2:859235/143‑1 (MQ=255)
                           tGGCAACGCTGAACGTTCCTTAACCACGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTCCCTGC‑GTTGACTTTTTCCTCGCTCGCGGT‑GTGACCCCGCGCGTTATGGATACGCGTATGa                                                                                                            <  1:447095/143‑1 (MQ=255)
                               aaCGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGGGTTATGGATGCGCGTACGACAcc                                                                                                        >  2:592706/1‑143 (MQ=255)
                                      aaGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCgcctggc                                                                                                 <  1:59404/143‑1 (MQ=255)
                                        ggTACGTTAATCATTGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCgcctggcct                                                                                               <  2:572737/143‑1 (MQ=255)
                                         gTACGTGAATCATGGCTCATTATCAGGGTAAGAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCgcctggcctg                                                                                              <  1:267173/143‑1 (MQ=255)
                                                aaTCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTAATGTATGCGCGTATGACACCGCCTGGCCTGGATAAAt                                                                                       >  1:537753/1‑143 (MQ=255)
                                                        cTGGTTATCAGGGTGAATTTGCCGTGATTATCGGCCTGGGCCTGACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGa‑a                                                                              <  2:503829/143‑1 (MQ=255)
                                                        cTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACCCCGCCTGGCCTGGATAAATTACCCGa‑a                                                                              >  1:816056/1‑143 (MQ=255)
                                                           aTTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GCGACGCCGCGCGTTATGGATACGCGTATGTCACGGCCTGGCCCGGCAAAGCTACCCGA‑Aggc                                                                           >  2:914587/1‑141 (MQ=255)
                                                              aTCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGA‑AGCCGTa                                                                        <  2:264131/143‑1 (MQ=255)
                                                                  gggTAAAAATGCCTGCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGCCTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGA‑AGCCGTAGAAc                                                                    <  2:1052744/143‑1 (MQ=255)
                                                                     tAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGGCGCCGCGCGTTATGGATACGCGTATGACAACGCCTGGCCTGGATAAATTACGCGA‑AGCGGTAGAACggc                                                                 >  1:423450/1‑141 (MQ=255)
                                                                                  aTTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGA‑AGCCGTAGAACGCCACACGGGCAGTCt                                                    <  1:763121/143‑1 (MQ=255)
                                                                                  aTTATCGGCCTGGGCCTCACCGGGCTTTACTGC‑GTGGACTTTTTCCTCGCTCGTGGT‑GTGACGCCGCGCGGTATGGATACACGTAGGACACCGCCTGGCCTGGATAAAGTACACGA‑AGCGGCAGATCGCCACACgcgcggcag                                                    >  2:50586/1‑138 (MQ=255)
                                                                                   ttATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTTTTTTTTTTTGGT‑GTGACGCGGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGGTAGATTACCCGG‑AGCCGTAGAACGCCACACGGGCAGTCTg                                                   >  2:273583/1‑143 (MQ=255)
                                                                                     agcGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGA‑AGCCGTAGAACGCCACACGGGCAGTCTGaa                                                 <  1:709864/141‑1 (MQ=255)
                                                                                      tCGGCCTGTGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGA‑AGCCGTAGAACGCCACACGGGCAGTCTGAat                                                <  2:1001468/143‑1 (MQ=255)
                                                                                      tCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGA‑AGCCGTAGAACGCCACACGGGCAGTCTGAat                                                <  1:856024/143‑1 (MQ=255)
                                                                                       cGGCCTGGGCCTCACCGGGATTTCCTGC‑GTTGACTTTTTCCTCCCTCGCGGT‑GTGGCGCCGCGCGTTCTGGATACGCGTATGACACCGCCTGGCCTGAACAAATTACCAGG‑AGCCGATGCACGCCACACCGGCACACGGAatg                                               >  2:441879/1‑143 (MQ=255)
                                                                                                ccTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GCGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATACATTACCCGA‑AGCCGTAGCACGCCACACGGGCACGCTGACtgctgaccggc                                      >  2:275938/1‑138 (MQ=255)
                                                                                                 ctgcccGGCCTTTCCTCC‑GTGGACGTTTTCCGCGCTCGCGGT‑GTGACGCCGCGGGTTATGGATACGCGTATGACACCGCCGGGCCTGGATAAATTACCCGA‑AGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCt                                     <  2:148251/139‑1 (MQ=255)
                                                                                                        ggCTTCCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGA‑AGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATggcg                              <  1:580697/143‑1 (MQ=255)
                                                                                                         ttctgccTGCGGTGGACTGTGTCCTCGCGCGGGGG‑GTCGCCCGGCGAGTGATGGATACGCGTACGACACCGCCTGGCCGGGATAAATTACCCGA‑AGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATggcg                              <  2:735219/138‑1 (MQ=255)
                                                                                                              ccTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGA‑AGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGAt                        >  1:939151/1‑143 (MQ=255)
                                                                                                               cTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCATGGCTAAAGTACCCGA‑AGCCGTAGAACGCCACACGGGCAGGCGGAAGGATGAAGGGCTGATGGCGGCAGATc                       >  1:10284/1‑143 (MQ=255)
                                                                                                                  c‑GTGGACTTTTCCCGCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGA‑AGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGa                    <  1:934519/143‑1 (MQ=255)
                                                                                                                     tGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGCCTCAATTACCCCACGGTCGCAG‑AGGCGACGCGCGCCGCCTGAATGATGAATGGCTGCTGGCGGCAGATCTGAtt                  >  1:488053/1‑143 (MQ=255)
                                                                                                                     tGGACTCTTGCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGA‑AGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGAtt                  <  1:927106/143‑1 (MQ=255)
                                                                                                                                 taacgcgcgGC‑GTCACACCGCGCGTTATGGATACGCGTATGACCCCGCCTCGCCTGGATAACTTACCCGA‑AGCCGTAGAACGCCACACGGGCAGTCTGGATGATGACTGGCTAATGGCGGCAGATCTGAGTGCCGCCAGTccc      >  2:56470/6‑143 (MQ=255)
                                                                                                                                  cgctcgcGGT‑GGGGCGCCGCGCGTTATGGATACGCTTAGGACGCCGCCTGGCCTGGATAAATTACCGGA‑CGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCGCCAGTCCCg     <  2:311096/143‑1 (MQ=255)
                                                                                                                                   gctcgcGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGA‑AGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCGCCAGTCCCgg    >  2:147886/1‑143 (MQ=255)
                                                                                                                                     gcgcgGC‑GTCACACCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGA‑AGCCGTAGAACACCACACAGGCACTCTGAATAAAAAATGACTAATGGCGGCAGATCTGAATGTCGCAAGTCCCGGTa  >  1:933861/2‑143 (MQ=255)
                                                                                                                                      tgcgGTGGTGACGCCGCGTGTTAGGGAGACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGA‑AGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCGCCAGTCCCGGTa  <  2:589676/142‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                      
TTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGC‑GTGGACTTTTTCCTCGCTCGCGGT‑GTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTGGATAAATTACCCGA‑AGCCGTAGAACGCCACACGGGCAGTCTGAATGATGAATGGCTGATGGCGGCAGATCTGATTGTCGCCAGTCCCGGTA  >  NC_000913/97035‑97309

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑