Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,414,112 G→A 18.9% L52L (TTG→TTA rlmB → 23S rRNA (Gm2251)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,414,1120GA18.9% 83.8 / 13.0 37L52L (TTG→TTArlmB23S rRNA (Gm2251)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (20/10);  new base A (3/4);  total (23/14)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCG  >  W3110S.gb/4414044‑4414170
                                                                    |                                                          
ggCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTgg                                                           >  1:1504421/1‑70 (MQ=255)
        aGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGcc                                                  <  1:1027938/71‑1 (MQ=255)
          aTAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCaa                                                >  1:393884/1‑71 (MQ=255)
          aTAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCaa                                                >  1:6592/1‑71 (MQ=255)
          aTAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCaa                                                >  1:896048/1‑71 (MQ=255)
          aTAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCaa                                                >  1:1538127/1‑71 (MQ=255)
          aTAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCaa                                                >  1:935358/1‑71 (MQ=255)
          aTAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCaa                                                >  1:1198543/1‑71 (MQ=255)
          aTAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCaa                                                >  1:975995/1‑71 (MQ=255)
          aTAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCaa                                                >  1:1069893/1‑71 (MQ=255)
                   tGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcgac                                       >  1:1263303/1‑71 (MQ=255)
                   tGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcgac                                       >  1:967683/1‑71 (MQ=255)
                   tGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcgac                                       >  1:1591543/1‑71 (MQ=255)
                   tGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcgac                                       >  1:705715/1‑71 (MQ=255)
                   tGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTcgac                                       >  1:137069/1‑71 (MQ=255)
                                          gAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTg                  <  1:1753479/69‑1 (MQ=255)
                                             tctcAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCAt               >  1:932690/1‑69 (MQ=255)
                                             tctcAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCAt               >  1:702642/1‑69 (MQ=255)
                                             tctcAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCAt               >  1:243470/1‑69 (MQ=255)
                                             tctcAGGGCGTGGTTATCCAGTTAGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGt                        <  1:1715096/60‑1 (MQ=255)
                                             tctcAGGGCGTGGTTATCCAGTTAGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGt                        <  1:513298/60‑1 (MQ=255)
                                             tctcAGGGCGTGGTTATCCAGTTAGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGt                        <  1:534035/60‑1 (MQ=255)
                                             tctcAGGGCGTGGTTATCCAGTTAGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGt                        <  1:584144/60‑1 (MQ=255)
                                               tcAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAgg           >  1:1028876/1‑71 (MQ=255)
                                               tcAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAgg           >  1:381358/1‑71 (MQ=255)
                                               tcAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAgg           >  1:1780481/1‑71 (MQ=255)
                                               tcAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAgg           >  1:1565067/1‑71 (MQ=255)
                                               tcAGGGCGTGGTTATCCAGTTAGCAAACCGCCAATGTCTCGACGAGAAAAGCGACGGTGc                      >  1:1627680/1‑60 (MQ=255)
                                               tcAGGGCGTGGTTATCCAGTTAGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGc                      >  1:887565/1‑60 (MQ=255)
                                               tcAGGGCGTGGTTATCCAGTTAGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGc                      >  1:1241152/1‑60 (MQ=255)
                                                 aGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAggg          <  1:486136/70‑1 (MQ=255)
                                                 aGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAggg          <  1:25794/70‑1 (MQ=255)
                                                   ggCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAggg          <  1:17402/68‑1 (MQ=255)
                                                         gTTATCCTGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCg  <  1:1917331/70‑1 (MQ=255)
                                                         gTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCg  <  1:469667/70‑1 (MQ=255)
                                                         gTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCg  <  1:1440839/70‑1 (MQ=255)
                                                         gTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCg  <  1:117511/70‑1 (MQ=255)
                                                               cAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAggg          <  1:15112/56‑1 (MQ=255)
                                                                    |                                                          
GGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTCAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCG  >  W3110S.gb/4414044‑4414170

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: