Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 669,029 C→T 51.4% E37K (GAA→AAA)  ybeA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb669,0290CT51.4% ‑5.5 / 50.6 35E37K (GAA→AAA) ybeAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (6/11);  new base T (14/4);  total (20/15)
Fisher's exact test for biased strand distribution p-value = 1.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

TGCTCACCCTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACC  >  W3110S.gb/668962‑669093
                                                                   |                                                                
tGCTCACCCTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAAt                                                               >  1:1491555/1‑71 (MQ=255)
tGCTCACCCTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAAt                                                               >  1:1628488/1‑71 (MQ=255)
       ccTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCg                                                        <  1:101457/71‑1 (MQ=255)
       ccTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCg                                                        <  1:2103125/71‑1 (MQ=255)
          ctTTGTCGAGTATGCGCTTGATGTCCGTATTCTTGCCGCGTTTTCCGGCCGGAATTTTAAt                                                               >  1:2708556/1‑61 (MQ=255)
          ctTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTTAAt                                                               >  1:3537624/1‑61 (MQ=255)
           tttGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGAAgg                                                    >  1:189434/1‑71 (MQ=255)
             tGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGAAggg                                                   <  1:3609992/70‑1 (MQ=255)
                        cgcTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGAAGGGCatat                                              >  1:379922/1‑64 (MQ=255)
                        cgcTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGAAGGGCatat                                              >  1:804419/1‑64 (MQ=255)
                        cgcTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGAAGGGCATATCTTTCg                                        >  1:439215/1‑70 (MQ=255)
                        cgcTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGAAGGGCATATCTTTCg                                        >  1:814838/1‑70 (MQ=255)
                          cTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGaa                                     <  1:658108/71‑1 (MQ=255)
                                 tCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGa                                                       <  1:909884/46‑1 (MQ=255)
                                  ccGCATTCTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGa                                                       <  1:2250011/45‑1 (MQ=255)
                                   cGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGc                             <  1:245961/70‑1 (MQ=255)
                                   cGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGc                             <  1:30530/70‑1 (MQ=255)
                                   cGCATACTTGCCGCGTTTTCCGGCCGGAATTTTAATCa                                                             >  1:3065228/1‑38 (MQ=38)
                                    gCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAg                           <  1:239296/71‑1 (MQ=255)
                                         cTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGAAGGGCATATCTTTCGGAAAacga                                >  1:1085450/1‑61 (MQ=255)
                                         cTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGAAGGGCATATCTTTCGGAAAacga                                >  1:1410659/1‑61 (MQ=255)
                                         cTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGAAGGGCATATCTTTCGGAAAacga                                >  1:2979464/1‑61 (MQ=255)
                                         cTTGCCGCGTTTTCCGGCCGGAATTTTAATCAGCTCGAAGGGCATATCTTTCGGAAAacg                                 >  1:1832339/1‑60 (MQ=255)
                                         cTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACt                      >  1:194432/1‑71 (MQ=255)
                                         cTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACt                      >  1:1411127/1‑71 (MQ=255)
                                                  ttttCCGGCCGGAATTTTAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAAc             >  1:1733954/1‑71 (MQ=255)
                                                    ttCCGGCCGGAGTTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCg           >  1:256864/1‑71 (MQ=255)
                                                    ttCCGGCCGGAATTTTAATCAGCTCGAAGGGCATATCTTTCGGaaaa                                   >  1:1583159/1‑47 (MQ=255)
                                                              aaTTTTAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTAcc  <  1:443093/70‑1 (MQ=255)
                                                               aTTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTAcc  >  1:1685692/1‑69 (MQ=255)
                                                                tttCAATCAGCTCGAAGGGCATATCTTTCGGAAAacg                                 <  1:2944144/37‑1 (MQ=255)
                                                                tttCAATCAGCTCGAAGGGCATATCTTTCGGAAAacg                                 <  1:3488815/37‑1 (MQ=255)
                                                                tttCAATCAGCTCGAAGGGCATATCTTTCGGAAAacg                                 <  1:2896430/37‑1 (MQ=255)
                                                                tttCAATCAGCTCGAAGGGCATATCTTTCGGAAAacg                                 <  1:1713901/37‑1 (MQ=255)
                                                                tttCAATCAGCTCGAAGGGCATATCTTTCGGAAAacg                                 <  1:20115/37‑1 (MQ=255)
                                                                   |                                                                
TGCTCACCCTCTTTGTCGAGTATGCGCTTGATGTCCGCATTCTTGCCGCGTTTTCCGGCCGGAATTTCAATCAGCTCGAAGGGCATATCTTTCGGAAAACGACGCAGGTACTCGGTAAAACCGGTTTGTACC  >  W3110S.gb/668962‑669093

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