Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 96,645 T→C V215A (GTG→GCG)  mraY → phospho‑N‑acetylmuramoyl‑pentapeptide transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091396,6450TC97.5% 135.1 / ‑6.2 40V215A (GTG→GCG) mraYphospho‑N‑acetylmuramoyl‑pentapeptide transferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (20/19);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGATGCCGCAGCTGGGGCTGTTCTACATTCTGCTGGCTTACTTCGTCATTGTGGGTACTGGCAACGCGGTAAACCTGACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGTGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTGGGCTTCCTGTGGTTTAACACCTATCCGGCGCAGGT  >  NC_000913/96503‑96789
                                                                                                                                              |                                                                                                                                                
gtgATGCCGCAGCTGGGGCTGTTCTACATTCTGCTGGCTTACTTCGTCATTGTGGGTACTGGCAACGCGGTAAACCTGACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTggg                                                                                                                                          <  1:612792/151‑1 (MQ=255)
      ccGCAGCTGGGGCTGTTCTACATTCTGCTGGCTTACTTCGTCATTGTGGGTACTGGCAACGCGGTAAACCTGACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCg                                                                                                                                    <  1:52551/151‑1 (MQ=255)
                 gCTGTTCTACATTCTGCTGGCTTACTTCGTCATTGTGGGTACTGGCAACGCGGTAAACCTGACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAAc                                                                                                                         >  1:575735/1‑151 (MQ=255)
                      tCTACATTCTGCTGGCTTACTTCGTCATTGTGGGTACTGGCAACGCGGTAAACCTGACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGc                                                                                                                    >  2:134050/1‑151 (MQ=255)
                         aCATTCTGCTGGCTTACTTCGTCATTGTGGGTACTGGCAACGCGGTAAACCTGACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAg                                                                                                                 <  2:47286/151‑1 (MQ=255)
                               tgctgGCTTACTTCGTCATTGTGGGTACTGGCAACGCGGTAAACCTGACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACtt                                                                                                           >  1:453905/1‑151 (MQ=255)
                                   ggCTTACTTCGTCATTGTGGGTACTGGCAACGCGGTAAACCTGACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCat                                                                                                       <  2:558230/151‑1 (MQ=255)
                                      ttACTTCGTCATTGTGGGTACTGGCAACGCGGTAAACCTGACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCatat                                                                                                     >  1:140426/1‑150 (MQ=255)
                                           tCGTCATTGTGGGTACTGGCAACGCGGTAAACCTGACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGt                                                                                                <  1:134053/150‑1 (MQ=255)
                                                   gtgGGTACTGGCAACGCGGTAAACCTGACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGTCTGTGGTTTTGCGCTGGTGGCGTGGGCGACCGGCAATATGAACTTTGCCCGCTACTTGCATATACCGTATCTGCgc                                                                                        <  2:74635/150‑2 (MQ=255)
                                                      ggTACTGGCAACGCGGTAAACCTGACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACg                                                                                    <  1:449308/151‑1 (MQ=255)
                                                           tGGCAACGCGGTAAACCTGACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCggg                                                                               >  1:259318/1‑151 (MQ=255)
                                                           tGGCAACGCGGTAAACCTGACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCggg                                                                               >  1:352082/1‑151 (MQ=255)
                                                               aaCGCGGTAAACCTGACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAAc                                                                           >  1:182626/1‑151 (MQ=255)
                                                                    ggTAAACCTGACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGtt                                                                      <  1:615818/151‑1 (MQ=255)
                                                                             gACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGt                                                             <  1:134974/151‑1 (MQ=255)
                                                                              aCCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGt                                                             >  1:168107/1‑150 (MQ=255)
                                                                                cGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTAcc                                                          >  2:354861/1‑151 (MQ=255)
                                                                                        tcGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGt                                                  >  2:221769/1‑151 (MQ=255)
                                                                                         cGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAg                                                   <  2:449975/149‑1 (MQ=255)
                                                                                         cGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTc                                                 >  2:81030/1‑151 (MQ=255)
                                                                                             ggCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCgggg                                             <  1:556007/151‑1 (MQ=255)
                                                                                                   gCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGAc                                       <  2:215105/151‑1 (MQ=255)
                                                                                                     aaTTATGCCGACCGAATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTg                                     >  2:548361/1‑151 (MQ=255)
                                                                                                      aTTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTgg                                    >  2:267186/1‑151 (MQ=255)
                                                                                                        tATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTggg                                   <  2:607945/150‑1 (MQ=255)
                                                                                                           gCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTGGGCtt                                <  1:544683/150‑1 (MQ=255)
                                                                                                             cGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTGGGCTTCCt                             >  2:193510/1‑151 (MQ=255)
                                                                                                             cGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTGGGCTTCCt                             <  2:188879/151‑1 (MQ=255)
                                                                                                                 cGTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTGGGCTTCCTGTgg                         <  1:81030/151‑1 (MQ=255)
                                                                                                                  gTATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTGGGCTTCCTGTGGt                        >  1:369050/1‑151 (MQ=255)
                                                                                                                   tATTTGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATTTGAACTTTGCCAGCTACTTGCAAATCCCGAATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCCCGATAGCCGGGGCAGGACTGGGCTTCCTGGGGtt                       >  1:404554/1‑151 (MQ=255)
                                                                                                                      ttGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTGGGCTTCCTGTGGTTTa                     >  2:621603/1‑150 (MQ=255)
                                                                                                                       tGTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTGGGCTTCCTGTGGTTTAac                   >  2:349361/1‑151 (MQ=255)
                                                                                                                        gTCGCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTGGGCTTCCTGTGGTTTAac                   <  2:515593/150‑1 (MQ=255)
                                                                                                                           gCCGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTGGGCTTCCTGTGGTTTAACACCt               <  2:575720/151‑1 (MQ=255)
                                                                                                                            ccGGTGGTTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTGGGCTTCCTGTGGTTTAACACCt               >  1:343333/1‑150 (MQ=255)
                                                                                                                              ggtggtTTTGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTGGGCTTCCTGTGGTTTAACACCTATc            <  2:596449/151‑1 (MQ=255)
                                                                                                                                    tttGCGCTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTGGGCTTCCTGTGGTTTAACACCTATCCGgcg       <  2:493027/150‑1 (MQ=255)
                                                                                                                                        cgcTGGCGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTGGGCTTCCTGTGGTTTAACACCTATCCGGCGCAGGt  >  1:265690/1‑151 (MQ=255)
                                                                                                                                              |                                                                                                                                                
GTGATGCCGCAGCTGGGGCTGTTCTACATTCTGCTGGCTTACTTCGTCATTGTGGGTACTGGCAACGCGGTAAACCTGACCGATGGTCTCGACGGCCTGGCAATTATGCCGACCGTATTTGTCGCCGGTGGTTTTGCGCTGGTGGCGTGGGCGACCGGCAATATGAACTTTGCCAGCTACTTGCATATACCGTATCTGCGACACGCCGGGGAACTGGTTATTGTCTGTACCGCGATAGTCGGGGCAGGACTGGGCTTCCTGTGGTTTAACACCTATCCGGCGCAGGT  >  NC_000913/96503‑96789

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