Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,332,723 A→G 22.2% A84A (GCT→GCC rpmA ← 50S ribosomal subunit protein L27

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,332,7230AG22.2% 86.0 / 22.7 45A84A (GCT→GCCrpmA50S ribosomal subunit protein L27
Reads supporting (aligned to +/- strand):  ref base A (17/18);  new base G (7/3);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAG  >  W3110S.gb/3332656‑3332793
                                                                   |                                                                      
gAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCt                                                                     >  1:604647/1‑71 (MQ=255)
  aTGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTc                                                                   <  1:3727992/71‑1 (MQ=255)
     tAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCGGCTTCGAt                                                                >  1:4296272/1‑71 (MQ=255)
     tAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCGGCTTCGAt                                                                >  1:4301697/1‑71 (MQ=255)
           gCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGAt                                                          <  1:405992/71‑1 (MQ=255)
              ccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGc                                                                      <  1:1831564/56‑1 (MQ=255)
              ccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATaaa                                                       >  1:1676833/1‑71 (MQ=255)
               cGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAAt                                                      >  1:1485348/1‑71 (MQ=255)
                gCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAAtt                                                     <  1:1622716/71‑1 (MQ=255)
                 cAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGAt                                                                >  1:627777/1‑59 (MQ=255)
                    cGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCGGCTTCGa                                                                 >  1:3114532/1‑55 (MQ=255)
                    cGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCGGCTTCGa                                                                 >  1:3551592/1‑55 (MQ=255)
                        ttGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTc                                                                   <  1:7663/49‑1 (MQ=255)
                        ttGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATaaa                                                       <  1:241965/61‑1 (MQ=255)
                        ttGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGttt                                             <  1:1678597/71‑1 (MQ=255)
                          gCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACgg                                                >  1:2871815/1‑66 (MQ=255)
                           cGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGttttt                                           <  1:1211454/70‑1 (MQ=255)
                           cGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTc                                          <  1:1412522/71‑1 (MQ=255)
                             gggCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTc                                                                   <  1:1565689/44‑1 (MQ=255)
                                cTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGATTCGATGCTGa                                                           >  1:3287550/1‑49 (MQ=255)
                                 tttCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTa                                                   >  1:4394239/1‑56 (MQ=255)
                                  ttCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTc                                                                   <  1:2818303/39‑1 (MQ=255)
                                  ttCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTc                                                                   <  1:619328/39‑1 (MQ=255)
                                   tCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCttt                                  >  1:1349629/1‑71 (MQ=255)
                                    cATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGa                                                                 >  1:3464036/1‑39 (MQ=255)
                                    cATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCttt                                  >  1:3043283/1‑70 (MQ=255)
                                       ccGTTACCGGGACGCGAAAAACTTATTCAGCTTCGa                                                                 >  1:1535593/1‑36 (MQ=255)
                                       ccGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACg                                                 >  1:2158847/1‑52 (MQ=255)
                                             ccGGGACGCGAAAAACTTATTCGGCTTCGATGCTGATaa                                                        <  1:896515/39‑1 (MQ=255)
                                                ggACGCGAAAAACTTATTCGGCTTCGATGCTGATAAATTTACg                                                 <  1:1376232/43‑1 (MQ=255)
                                                       aaaaaCTTATTCGGCTTCGATGCTGATAAATTTACGGttttt                                           >  1:292751/1‑42 (MQ=255)
                                                       aaaaaCTTATTCGGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGa                          <  1:4157505/59‑1 (MQ=255)
                                                       aaaaaCTTATTCAGCTTCGATGCTGATAAATTTACGGttt                                             <  1:4217132/40‑1 (MQ=255)
                                                       aaaaaCTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTc                                          <  1:3894809/43‑1 (MQ=255)
                                                       aaaaaCTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCg                           <  1:3403187/58‑1 (MQ=255)
                                                       aaaaaCTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTa               <  1:1831396/70‑1 (MQ=255)
                                                         aaaCTTATTCAGCTTCGATGCTGATAAATTTACgg                                                >  1:2197119/1‑35 (MQ=255)
                                                         aaaCTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGc                                      >  1:4053679/1‑45 (MQ=255)
                                                         aaaCTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGc                                      >  1:3420049/1‑45 (MQ=255)
                                                         aaaCTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCa                    >  1:4413652/1‑63 (MQ=255)
                                                           aCTTATTCGGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTc          >  1:1265792/1‑71 (MQ=255)
                                                           aCTTATTCGGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTc          >  1:464354/1‑71 (MQ=255)
                                                                ttCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTcttctt     <  1:3534122/71‑1 (MQ=255)
                                                                ttCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCtt     <  1:58421/71‑1 (MQ=255)
                                                                  cAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTAc              >  1:2773041/1‑60 (MQ=255)
                                                                   aGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCtt     >  1:315919/1‑68 (MQ=255)
                                                                   aGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAg  >  1:3790464/1‑71 (MQ=255)
                                                                   |                                                                      
GAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAG  >  W3110S.gb/3332656‑3332793

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: