Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,389,895 G→A *322* (TAG→TAA
P235L (CCT→CTT) 
ycjG →
mpaA ←
L‑Ala‑D/L‑Glu epimerase
murein peptide amidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,389,8950GA97.0% 95.8 / ‑4.4 33*322* (TAG→TAA
P235L (CCT→CTT) 
ycjG
mpaA
L‑Ala‑D/L‑Glu epimerase
murein peptide amidase A
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (15/17);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CCTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTGCCGCAGGTCAGTTTCGCCGATCTTGACGGACCGACCTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAGGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGCTTCGTCGGAGGAGATGGGCGGAAACTCTGCGGTGATGCAATGCAAGTTCAGATCTGCACACCAACTGCCGAAAGAGCCAGGCGTTT  >  NC_000913/1389765‑1390030
                                                                                                                                  |                                                                                                                                       
ccTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTGCCGCAGGTCAGTTTCGCCGATCTTGACGGACCGACCTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAg                                                                                                                                 <  1:462258/139‑1 (MQ=255)
     cGTGCCATTAGCGCCGCTTTACCGCTGGTGCCGCAGGTCAGTTTCGCCGATCTTGACGGACCGACCTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCgcag                                                                                                                            >  1:387764/1‑139 (MQ=255)
             tAGCGCCGCTTTACCGCTGGTGCCGCAGGTCAGTTTCGCCGATCTTGACGGACCGACCTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTg                                                                                                                    <  1:406683/139‑1 (MQ=255)
                     cTTTACCGCTGGTGCCGCAGGTCAGTTTCGCCGATCTTGACGGACCGACCTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCa                                                                                                            >  2:383798/1‑139 (MQ=255)
                                    cGCAGGTCAGTTTCGCCGATCTTGACGGACCGACCTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCg                                                                                             <  2:42783/139‑1 (MQ=255)
                                         gTCAGTTTCGCCGATCTTGACGGACCGACCTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGc                                                                                        <  1:165178/139‑1 (MQ=255)
                                                            aCGGACCGACCTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGCTTCGTCGGAGGAGATGGGc                                                                     >  2:321058/1‑139 (MQ=255)
                                                                accgaccTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAgcgc                                                                                                                              >  2:2419/1‑78 (MQ=255)
                                                                accgaccTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAgcgc                                                                                                                              <  1:2419/78‑1 (MQ=255)
                                                                accgaccTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGCTTCGTCGGAGGAGATGGGCGGaa                                                                 >  2:1268857/1‑139 (MQ=255)
                                                                 ccgaccTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCgca                                                                                                                             <  1:479560/78‑1 (MQ=255)
                                                                 ccgaccTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCgca                                                                                                                             >  2:479560/1‑78 (MQ=255)
                                                                 ccgaccTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCaaa                                                                                                          >  1:496933/1‑97 (MQ=255)
                                                                 ccgaccTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCaaa                                                                                                          <  2:496933/97‑1 (MQ=255)
                                                                  cgaccTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGc                                                                                        <  2:1560336/114‑1 (MQ=255)
                                                                  cgaccTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGc                                                                                        >  1:1560336/1‑114 (MQ=255)
                                                                              cGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGCTTCGTCGGAGGAGATGGGCGGAAACTCTGCGGTGatg                                                   >  2:441992/1‑139 (MQ=255)
                                                                                    aTGTGGAACCGGCGCTTCAGATCACGACGGGCGAAATGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGCTTCGTCGGAGGAGATGGGCGGAAACTCTGCGGTGatgcaatgc                                             <  1:441992/139‑1 (MQ=255)
                                                                                        ggAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCgca                                                                                                                             >  2:1045655/1‑55 (MQ=255)
                                                                                        ggAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCgca                                                                                                                             <  1:1045655/55‑1 (MQ=255)
                                                                                         gAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGCTTCGTCGGAGGAGATGGGCGGAAACTCTGCGGTGATGCAATGCAAGtt                                        <  2:387764/139‑1 (MQ=255)
                                                                                               gcgcTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCaaa                                                                                                          >  1:1506229/1‑67 (MQ=255)
                                                                                               gcgcTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCaaa                                                                                                          <  2:1506229/67‑1 (MQ=255)
                                                                                               gcgcTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGCTTCGTCGGAGGAGATGGGCgg                                                                   <  1:843510/106‑1 (MQ=255)
                                                                                               gcgcTTCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGCTTCGTCGGAGGAGATGGGCgg                                                                   >  2:843510/1‑106 (MQ=255)
                                                                                                    tCAGTTCACGACGGGCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGCTTCGTCGGAGGAGATGGGCGGAAACTCTGCGGTGATGCAATGCAAGTTCAGATCTGcac                             >  2:174418/1‑139 (MQ=255)
                                                                                                                  gCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGCTTCGTCGGAGGAGATGGGCGGAAACTCTg                                                           >  2:473992/1‑95 (MQ=255)
                                                                                                                  gCGAATTGCATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGCTTCGTCGGAGGAGATGGGCGGAAACTCTg                                                           <  1:473992/95‑1 (MQ=255)
                                                                                                                          cATCTTTAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGCTTCGTCGGAGGAGATGGGCGGAAACTCTGCGGTGATGCAATGCAAGTTCAGATCTGCACACCAACTGCCGAAAGAGCCAgg       <  1:174418/139‑1 (MQ=255)
                                                                                                                               ttAGGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTTGCTTCGTCGGAGGAGATGGGCGGAAACTCTGCGGTGATGCAATGCAAGTTCAGATCTGCACACCAACTGCCGAAAGAGCCAGGCGttt  >  2:531287/1‑139 (MQ=255)
                                                                                                                               ttAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCg                                                                                             <  2:197289/48‑1 (MQ=255)
                                                                                                                               ttAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGCTTCGTCggag                                                                              <  1:526207/63‑1 (MQ=255)
                                                                                                                               ttAAGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGCTTCGTCggag                                                                              >  2:526207/1‑63 (MQ=255)
                                                                                                                                  |                                                                                                                                       
CCTCTCGTGCCATTAGCGCCGCTTTACCGCTGGTGCCGCAGGTCAGTTTCGCCGATCTTGACGGACCGACCTGGCTGGCGGTAGATGTGGAACCGGCGCTTCAGTTCACGACGGGCGAATTGCATCTTTAGGATGCCAGCGCAGCAAGTTTGCCATCGCAAAGAGATATTTTTCGCTGGCTTCGTCGGAGGAGATGGGCGGAAACTCTGCGGTGATGCAATGCAAGTTCAGATCTGCACACCAACTGCCGAAAGAGCCAGGCGTTT  >  NC_000913/1389765‑1390030

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