Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,068,568 T→A 78.0% G22G (GGT→GGA ymdF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,068,5680TA78.0% 53.4 / 10.1 27G22G (GGT→GGAymdFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base A (5/16);  total (8/19)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTG  >  W3110S.gb/1068503‑1068626
                                                                 |                                                          
aTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaa                                                        >  1:756851/1‑70 (MQ=255)
aTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaa                                                        >  1:1917523/1‑70 (MQ=255)
 tGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                      <  1:2102792/71‑1 (MQ=255)
 tGGCAAACCATCGAGGCGGTTCCGGCAATTTTGC‑‑AAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAgg                                                    <  1:1246321/71‑1 (MQ=255)
   gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAgg                                                    >  1:465668/1‑71 (MQ=255)
          aTCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGacag                                             <  1:934004/71‑1 (MQ=255)
                gCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGacagcacagc                                       <  1:2315518/71‑1 (MQ=255)
                     tCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCggggg                                  >  1:2540071/1‑71 (MQ=255)
                     tCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCggggg                                  >  1:931914/1‑71 (MQ=255)
                     tCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAg                                                     >  1:428369/1‑52 (MQ=255)
                      ccGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaa                                                        <  1:1451959/48‑1 (MQ=255)
                        ggCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGacag                                             <  1:1545136/57‑1 (MQ=255)
                         gCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaa                                                        <  1:1490229/45‑1 (MQ=255)
                         gCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAAt                              >  1:896796/1‑71 (MQ=255)
                             ttttGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGacagcac                                          <  1:695422/55‑1 (MQ=255)
                             ttttGCAGAAGACCGCGAAAGAGCATCAGAAGCAGTaaaaa                                                       <  1:967493/41‑1 (MQ=38)
                                tGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAg                        <  1:1302307/70‑1 (MQ=255)
                                tGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAg                        <  1:1117081/70‑1 (MQ=255)
                                  cAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGatt                     <  1:663540/71‑3 (MQ=255)
                                  cAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAAt                     <  1:460794/71‑1 (MQ=255)
                                    gaagaCCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGa                   <  1:2129463/71‑1 (MQ=255)
                                                      tCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAg                        <  1:427889/48‑1 (MQ=255)
                                                      tCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGaa                      <  1:2566478/50‑1 (MQ=255)
                                                      tCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGaa                      <  1:1161177/50‑1 (MQ=255)
                                                      tCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCa      <  1:2724279/66‑1 (MQ=255)
                                                      tCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTg  >  1:500979/1‑70 (MQ=255)
                                                             cAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAgcgcg        <  1:1616121/57‑1 (MQ=255)
                                                                 |                                                          
ATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTG  >  W3110S.gb/1068503‑1068626

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: