Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,099,547 C→T 61.9% A143V (GCC→GTC)  ycdX → predicted zinc‑binding hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,099,5470CT61.9% 9.2 / 20.3 21A143V (GCC→GTC) ycdXpredicted zinc‑binding hydrolase
Reads supporting (aligned to +/- strand):  ref base C (5/3);  new base T (5/8);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

ACTATCGCCAGCGGCAATGTGCATATAATAAGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTC  >  W3110S.gb/1099477‑1099613
                                                                      |                                                                  
aCTATCGCCAGCGGCAATGTGCATATAATAAGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGc                                                                    >  1:1186217/1‑71 (MQ=255)
    tCGCCAGCGGCAATGTGCATATAATAAGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGTCGtt                                                                >  1:1056169/1‑71 (MQ=255)
        cAGCGGCAATGTGCATATAATAAGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGTCGTTGCTg                                                            <  1:1333778/71‑1 (MQ=255)
                        ataataAGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGTCGTTGCTGAAGCAGCCGCGa                                                <  1:1116346/67‑1 (MQ=255)
                        ataataAGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAagcag                                                      >  1:329719/1‑61 (MQ=255)
                         taataaGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGcc                                                                   >  1:1547567/1‑47 (MQ=255)
                           ataaGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGTCGTTGCTGAAGCAGCCGCGAAACATc                                          <  1:511375/70‑1 (MQ=255)
                            taaGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACATc                                          >  1:517399/1‑69 (MQ=255)
                            taaGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGAGAAACATc                                          >  1:306841/1‑69 (MQ=255)
                                     ccGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTg                               >  1:777148/1‑71 (MQ=255)
                                       ggAAATCCCAAATATGAAATAGATGTGAAAGTCGTTGCTGAAGCAGCCGCGa                                                <  1:809185/52‑1 (MQ=255)
                                       ggAAATCCCAAATATGAAATAGATGTGAAAGTCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCt                                <  1:722728/68‑1 (MQ=255)
                                       ggAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGa                                                <  1:1143516/52‑1 (MQ=255)
                                                     tGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGaa                                               <  1:1307874/39‑1 (MQ=255)
                                                       aaaTAGATGTGAAAGTCGTTGCTGAAGCAGCCGCGAAACATCAGGTGgcgc                                 >  1:311069/1‑51 (MQ=255)
                                                           aGATGTGAAAGTCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCAtttt         >  1:690468/1‑71 (MQ=255)
                                                            gATGTGAAAGTCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCAtttt         <  1:809993/70‑1 (MQ=255)
                                                            gATGTGAAAGTCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCAtttt         <  1:832385/70‑1 (MQ=255)
                                                             aTGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACATCTGGTGGCGCTGGAAATCaat                      <  1:871975/56‑1 (MQ=255)
                                                               gtgAAAGTCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCt                                >  1:1683046/1‑44 (MQ=255)
                                                               gtgAAAGTCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCt                                >  1:1681036/1‑44 (MQ=255)
                                                                tgAAAGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTaca     >  1:961754/1‑70 (MQ=255)
                                                                  aaaGTCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTc  <  1:156371/71‑1 (MQ=255)
                                                                      |                                                                  
ACTATCGCCAGCGGCAATGTGCATATAATAAGCCATCCCGGAAATCCCAAATATGAAATAGATGTGAAAGCCGTTGCTGAAGCAGCCGCGAAACATCAGGTGGCGCTGGAAATCAATAATTCCTCATTTTTACACTC  >  W3110S.gb/1099477‑1099613

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: