Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 833,307 G→T 57.6% T22T (ACC→ACA ybiA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb833,3070GT57.6% ‑3.5 / 37.1 26T22T (ACC→ACAybiAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (5/6);  new base T (9/6);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

GGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATT  >  W3110S.gb/833239‑833373
                                                                    |                                                                  
ggCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTTGt                                                                  >  1:528074/1‑71 (MQ=255)
ggCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGt                                                                  <  1:672916/71‑1 (MQ=255)
               tCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTTGTGCTGTAGAAATTTat                                                   <  1:1043117/71‑1 (MQ=255)
               tCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTTGTGCTGTAGAAATTTat                                                   >  1:1639737/1‑71 (MQ=255)
               tCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTTGTGCTGTAGAAATTTat                                                   >  1:585896/1‑71 (MQ=255)
                    cTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTTGTGCTGTAGAAATTTATGATgg                                              >  1:634991/1‑71 (MQ=255)
                     ttGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGgt                                             >  1:790506/1‑71 (MQ=255)
                           ggCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTTGTGCTGTAGAAATTTATGATGGTGTc                                          >  1:528461/1‑68 (MQ=255)
                             acAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTa                                                            >  1:825377/2‑48 (MQ=255)
                              cAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTTGTGCTGTAGAAATTTATGATGGTGTCTTGCAt                                    <  1:8253/71‑1 (MQ=255)
                                      taaGTTAGAAAAATCACCGTAGTCATCGCTTGTGCTGTAGaaa                                                        <  1:758027/42‑1 (MQ=255)
                                       aaGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGttt                           <  1:402023/71‑2 (MQ=255)
                                            aGAAAAATCACCGTAGTCATCGCTTGTGCTGTAGAAATTTATGATGGTGTCTTGc                                      >  1:137727/1‑55 (MQ=255)
                                            aGAAAAATCACCGTAGTCATCGCTTGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATc                       <  1:251934/70‑1 (MQ=255)
                                             gAAAAATCACCGTAGTCATCGCTTGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCt                     >  1:761935/1‑71 (MQ=255)
                                             gAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCt                     <  1:669726/71‑1 (MQ=255)
                                                aaaTCACCGTAGTCATCGCTGGTGCTGTAGAAAttt                                                     >  1:972719/1‑36 (MQ=255)
                                                aaaTCACCGTAGTCATCGCTGGTGCTGTAGAAAttt                                                     >  1:329150/1‑36 (MQ=255)
                                                      ccGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCg            <  1:156116/71‑1 (MQ=255)
                                                         tAGTCATCGCTTGTGCTGTAGAAATTTATGATGGTGTCTTGCAt                                    <  1:15913/44‑1 (MQ=255)
                                                            tCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACgt                               <  1:117947/46‑1 (MQ=255)
                                                              aTCGCTTGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTg                           >  1:1115466/1‑48 (MQ=255)
                                                              aTCGCTTGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCa    <  1:29324/71‑1 (MQ=255)
                                                               tCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCtt                    <  1:832093/54‑1 (MQ=255)
                                                                 gCTTGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCAtt  >  1:1631100/1‑70 (MQ=255)
                                                                 gCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTtg                  >  1:753136/1‑54 (MQ=255)
                                                                    |                                                                  
GGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATT  >  W3110S.gb/833239‑833373

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: