Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 665,576 C→T 14.8% L246L (CTG→CTA)  mrdB ← cell wall shape‑determining protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913665,5760CT14.8% 71.4 / 3.9 27L246L (CTG→CTA) mrdBcell wall shape‑determining protein
Reads supporting (aligned to +/- strand):  ref base C (14/9);  new base T (2/2);  total (16/11)
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

GAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGTAATCATGCATCAGGAAGAACCAC  >  NC_000913/665458‑665720
                                                                                                                      |                                                                                                                                                
gAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGc                                                                                                                   >  1:244957/1‑150 (MQ=255)
gAATGTAGAGAGCGAGCAGAATCAGAATCCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGc                                                                                                                   >  1:23244/1‑150 (MQ=255)
 aaTGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCTGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGcc                                                                                                                  >  2:86348/1‑150 (MQ=255)
 aaTGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGcc                                                                                                                  >  2:86340/1‑150 (MQ=255)
       gagagCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGc                                                                                                            <  2:103890/150‑1 (MQ=255)
                cagaatcagaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGa                                                                                                   <  2:48174/150‑1 (MQ=255)
                cagaatcagaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGa                                                                                                   <  2:27654/150‑1 (MQ=255)
                cagaatcagaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGa                                                                                                   <  2:53048/150‑1 (MQ=255)
                cagaatcagaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGa                                                                                                   <  2:53447/150‑1 (MQ=255)
                cagaatcagaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGa                                                                                                   <  1:192743/150‑1 (MQ=255)
                   aatcagaatGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTg                                                                                                >  1:229572/1‑150 (MQ=255)
                           tGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAAtgt                                                                                        >  1:87490/1‑150 (MQ=255)
                                            tctTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGt                                                                       >  1:82887/1‑150 (MQ=255)
                                              ttCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTgg                                                                     >  2:172994/1‑150 (MQ=255)
                                              ttCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTgg                                                                     >  2:39616/1‑150 (MQ=255)
                                              ttCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAAGTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTgg                                                                     >  2:172991/1‑150 (MQ=255)
                                                  gCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGTAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCt                                                                 >  2:88139/1‑150 (MQ=255)
                                                  gCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGTAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCt                                                                 >  2:127462/1‑150 (MQ=255)
                                                                        cAGTATGGCGTTCGGGGAGAAATTCATGCTGTGACTGTGTGCCGTGTAGCAAGCCTTTGCCGCGTAATCCGCCGGAGCATATAGCAATTTTAGTCTGATTAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGcttca                                           <  1:88139/150‑1 (MQ=255)
                                                                        cAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGTAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGcatca                                           <  1:127462/150‑1 (MQ=255)
                                                                           tATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTa                                        >  2:156006/1‑150 (MQ=255)
                                                                                  ttCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACgcgctg                                 <  2:93149/150‑1 (MQ=255)
                                                                                  ttCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACgcgctg                                 <  2:91505/150‑1 (MQ=255)
                                                                                  ggcGGGGAGACAGTCCAGCTGGGCCGGAGAGCCGTGCAGCCAGCCTTCGCCGCGTCATCCGCCGGAGCCACTAGCACTTTTCGCCTGAAGCCTGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGACCAGGAGCATCATTACgcgctg                                 <  2:91507/148‑1 (MQ=255)
                                                                                    cGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACgcgctggc                               >  2:56462/1‑150 (MQ=255)
                                                                                           aaaTTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGTa                        >  1:60509/1‑150 (MQ=255)
                                                                                                                 ccGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGTAATCATGCATCAGGAAGAACcac  >  1:164822/1‑150 (MQ=255)
                                                                                                                      |                                                                                                                                                
GAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCAGTATGGCGTTCGGGGAGAAATTCAAGCTGTGACTGAGTGCCGTGCAGCCAGCCTTTGCCGCGTAATCCGCCGGAGCCAATAGCAATTTTAGACTGAATAATGTGATAGCCCGCGCCGAGTGGGTCTGATTCCGGGTCCAGGAGCATCATTACGCGCTGGCGCTGGTAATCATGCATCAGGAAGAACCAC  >  NC_000913/665458‑665720

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