Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,883,198 A→G 22.5% F429L (TTC→CTC)  rumA ← 23S rRNA (uracil‑5)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,883,1980AG22.5% 74.5 / 20.9 40F429L (TTC→CTC) rumA23S rRNA (uracil‑5)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (13/18);  new base G (6/3);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

AACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTA  >  minE/1883128‑1883265
                                                                      |                                                                   
aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                     <  1:192118/71‑1 (MQ=255)
aaCAATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                     <  1:1478940/71‑1 (MQ=255)
  ccATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                     >  1:1587042/1‑69 (MQ=255)
        tCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCTAGAGAAg                                                                  >  1:2167994/1‑66 (MQ=255)
        tCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAg                                                                  >  1:1976482/1‑66 (MQ=255)
          ctctCCTTTAGGTACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATc                                                           <  1:1535057/71‑1 (MQ=255)
          ctctCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATc                                                           <  1:336605/71‑1 (MQ=255)
          ctctCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATc                                                           <  1:753933/71‑1 (MQ=255)
              ccTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaaaa                                                                   >  1:832375/1‑59 (MQ=255)
              ccTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaaaa                                                                   >  1:1323113/1‑59 (MQ=255)
               cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGa                                                         <  1:1442828/68‑1 (MQ=255)
               cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTc                                                      >  1:1886353/1‑71 (MQ=255)
                tttAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGt                                                                 <  1:729129/59‑1 (MQ=255)
                 ttAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                     <  1:1926004/54‑1 (MQ=255)
                   aGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCAta                                                           >  1:750280/1‑61 (MQ=255)
                   aGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCAt                                                            >  1:929322/1‑61 (MQ=255)
                   aGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATc                                                           >  1:1739600/1‑62 (MQ=255)
                       aCCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                     <  1:727212/48‑1 (MQ=255)
                       aCCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                     <  1:2051491/48‑1 (MQ=255)
                       aCCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCa                                                    <  1:1095673/65‑1 (MQ=255)
                       aCCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCa                                                    <  1:147322/65‑1 (MQ=255)
                         cAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                     <  1:1812781/46‑1 (MQ=255)
                         cAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCa                                                    <  1:607400/63‑1 (MQ=255)
                         cAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCa                                                    <  1:1287487/63‑1 (MQ=255)
                          aGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGTAAAGTACCATCGATTc                                                      <  1:796254/60‑1 (MQ=255)
                             ccTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgtg                                         <  1:1660006/70‑1 (MQ=255)
                              cTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCgtgtg                                       >  1:1419669/1‑71 (MQ=255)
                              cTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCgtgtg                                       >  1:1930120/1‑71 (MQ=255)
                                  cGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGa                                                  >  1:1731058/1‑56 (MQ=255)
                                   gAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgtgtg                                       >  1:661644/1‑66 (MQ=255)
                                   gAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTggg                                   >  1:1292422/1‑70 (MQ=255)
                                     aaTCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTAc                                                               <  1:1213410/40‑1 (MQ=255)
                                      aTCGGCAAATCGCAACTATTTAACGCGCGAGaaaa                                                                   >  1:1556475/1‑35 (MQ=255)
                                             aaTCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTc                                                      <  1:1784565/41‑1 (MQ=255)
                                               tCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCa                                                    <  1:120541/41‑1 (MQ=255)
                                               tCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCa                                                    <  1:2176973/41‑1 (MQ=255)
                                                cGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCAt                     >  1:1811235/1‑71 (MQ=255)
                                                   aaCTATTTAACGCGCGAGAAAAGTACCATCGATTCCa                                                    <  1:1281355/37‑1 (MQ=255)
                                                    aCTATTTAACGCGCGAGAAAAGAACCATCGATTCCAGATGTCCCgtgtgt                                      >  1:325842/1‑50 (MQ=255)
                                                       aTTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGc                       >  1:1867361/1‑62 (MQ=255)
                                                       aTTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGt              >  1:370666/1‑71 (MQ=255)
                                                       aTTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGt              >  1:1447195/1‑71 (MQ=255)
                                                          tAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgc           >  1:2120895/1‑71 (MQ=255)
                                                          tAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgc           >  1:1188574/1‑71 (MQ=255)
                                                            aCGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCgcgc         <  1:2106602/71‑1 (MQ=255)
                                                            aCGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCgcgc         <  1:106860/71‑1 (MQ=255)
                                                               cgcgAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAAt      <  1:101776/71‑1 (MQ=255)
                                                                gcgAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCTCGCAATg     <  1:148450/71‑1 (MQ=255)
                                                                  gagaAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGt   <  1:1929028/71‑1 (MQ=255)
                                                                    gaAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGta  >  1:20072/1‑70 (MQ=255)
                                                                      |                                                                   
AACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTA  >  minE/1883128‑1883265

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: