Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 96,008 G→A *453* (TAG→TAA
V3I (GTT→ATT) 
murF →
mraY →
UDP‑N‑acetylmuramoyl‑tripeptide:D‑alanyl‑D‑ alanine ligase
phospho‑N‑acetylmuramoyl‑pentapeptide transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091396,0080GA97.0% 96.9 / ‑6.6 33*453* (TAG→TAA
V3I (GTT→ATT) 
murF
mraY
UDP‑N‑acetylmuramoyl‑tripeptide:D‑alanyl‑D‑ alanine ligase
phospho‑N‑acetylmuramoyl‑pentapeptide transferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (20/12);  total (20/13)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAACTGCGTTAATTACGCGTCTTAAATTACTGATTGCTGAGCAACAGGTAATTACGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCAT  >  NC_000913/95876‑96142
                                                                                                                                    |                                                                                                                                      
tAAAACTGCGTTAATTACGCGTCTTAAATTACTGATTGCTGAGCAACAGGTAATTACGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTtggc                                                                                                                                  <  2:109905/139‑1 (MQ=255)
         gTTAATTACGCGTCTTAAATTACTGATTGCTGAGCAACAGGTAATTACGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAATTTGGCTGGCCGAAc                                                                                                                         >  2:23632/1‑139 (MQ=255)
           tAATTACGCGTCTTAAATTACTGATTGCTGAGCAACAGGTAATTACGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAATTTGGCTGGCCGAACAt                                                                                                                       <  2:55427/139‑1 (MQ=255)
                aCGCGTCTTAAATTACTGATTGCTGAGCAACAGGTAATTACGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGGAATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTc                                                                                                                 <  1:285838/139‑1 (MQ=255)
                 cgcgTCTTAAATTACTGATTGCTGAGCAACAGGTAATTACGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTc                                                                                                                 <  2:93846/139‑1 (MQ=255)
                 cgcgTCTTAAATTACTGATTGCTGAGCAACAGGTAATTACGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTc                                                                                                                 <  1:826582/139‑1 (MQ=255)
                     tCTTAAATTACTGATTGCTGAGCAACAGGTAATTACGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAat                                                                                                             >  2:886492/1‑139 (MQ=255)
                               cTGATTGCTGAGCAACAGGTAATTACGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCgg                                                                                                   >  1:472340/1‑139 (MQ=255)
                                                   aaTTACGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTc                                                                                                      <  1:151140/115‑1 (MQ=255)
                                                   aaTTACGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTc                                                                                                      >  2:151140/1‑115 (MQ=255)
                                                        cGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAATTTGGCTGGCTGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCATAACTGACg                                                                          >  2:932225/1‑139 (MQ=255)
                                                         gATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGt                                                                         >  2:682325/1‑139 (MQ=255)
                                                                      ggTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTa                                                                                              >  1:245613/1‑105 (MQ=255)
                                                                      ggTTCACGTAGTGCCGCCAGGGAAGATGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTa                                                                                              <  2:245613/105‑1 (MQ=255)
                                                                          cACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGc                                                        >  2:444860/1‑139 (MQ=255)
                                                                                 tGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGa                                                 <  2:472340/139‑1 (MQ=255)
                                                                                   ccgccATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGAcc                                               >  2:808337/1‑139 (MQ=255)
                                                                                          ggAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAATTtggc                                                                                                                                  >  2:159157/1‑49 (MQ=255)
                                                                                          ggAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAATTtggc                                                                                                                                  <  1:159157/49‑1 (MQ=255)
                                                                                                   agtaCGCGCTTTACAGGAGAATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCAt                               >  1:593330/1‑139 (MQ=255)
                                                                                                    gtaCGCGCTTTACAGGAGAATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCAtt                              <  1:682325/139‑1 (MQ=255)
                                                                                                           cTTTACAGGAGAATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCg                                                                                                    <  2:210273/62‑1 (MQ=255)
                                                                                                           cTTTACAGGAGAATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCg                                                                                                    >  1:210273/1‑62 (MQ=255)
                                                                                                              tACAGGAGAATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGc                    <  2:900913/139‑1 (MQ=255)
                                                                                                                   gagaATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGc                    <  2:52169/134‑1 (MQ=255)
                                                                                                                   gagaATGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGc                    >  1:52169/1‑134 (MQ=255)
                                                                                                                      aaTGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTAt            >  1:845410/1‑139 (MQ=255)
                                                                                                                        tGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGa          >  2:233553/1‑139 (MQ=255)
                                                                                                                        tGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGa          >  2:906296/1‑139 (MQ=255)
                                                                                                                        tGGGACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGa          >  2:958584/1‑139 (MQ=255)
                                                                                                                          ggACATGTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGAtt        >  2:888854/1‑139 (MQ=255)
                                                                                                                                gTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCAt  >  1:134261/1‑139 (MQ=255)
                                                                                                                                gTTAATTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCAt  >  1:796543/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                      
TAAAACTGCGTTAATTACGCGTCTTAAATTACTGATTGCTGAGCAACAGGTAATTACGATTTTAGTTAAGGGTTCACGTAGTGCCGCCATGGAAGAGGTAGTACGCGCTTTACAGGAGAATGGGACATGTTAGTTTGGCTGGCCGAACATTTGGTCAAATATTATTCCGGCTTTAACGTCTTTTCCTATCTGACGTTTCGCGCCATCGTCAGCCTGCTGACCGCGCTGTTCATCTCATTGTGGATGGGCCCGCGTATGATTGCTCAT  >  NC_000913/95876‑96142

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