Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,606,089 G→A 31.4% N37N (AAC→AAT dnaT ← DNA biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,606,0890GA31.4% 45.0 / 27.1 35N37N (AAC→AATdnaTDNA biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (10/14);  new base A (8/3);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

TTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCA  >  W3110S.gb/4606025‑4606153
                                                                |                                                                
tttCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGc                                                             >  1:621312/1‑70 (MQ=255)
    ttcCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTATTAGCAAATa                                                        >  1:3100613/1‑71 (MQ=255)
      ccAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATa                                                        <  1:68287/69‑1 (MQ=255)
         gcgcCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGc                                                    >  1:4369005/1‑70 (MQ=255)
                   ttCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTATTAGCAAATACGGCAACcac                                              <  1:4000847/66‑1 (MQ=255)
                     cAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCt                                       >  1:3138275/1‑71 (MQ=255)
                      aGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGt                                                                 <  1:1497249/44‑1 (MQ=255)
                         caggcgtgcaggcgtgACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCa                                    <  1:2916381/70‑1 (MQ=255)
                            gcgtgcaggcgtgACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCtt                                >  1:3306844/1‑71 (MQ=255)
                                 caggcgtgACGGCATAAAACGCCGGGGCATTATTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTg                             >  1:4537224/1‑69 (MQ=255)
                                  aggcgtgACGGCATAAAACGCCGGGGCATTATTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGc                            >  1:3188832/1‑69 (MQ=255)
                                   ggcgtgACGGCATAAAACGCCGGGGCATTATTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAg                         <  1:4648567/71‑1 (MQ=255)
                                     cgtgACGGCATAAAACGCCGGGGCATTATTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGaa                       <  1:1707522/71‑1 (MQ=255)
                                         aCGGCATAAAACGCCGGGGCATTATTAGCAAATACGGCAACca                                               >  1:2233019/1‑43 (MQ=255)
                                           ggCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGttt                 <  1:4767357/71‑1 (MQ=255)
                                              aTAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTgg               >  1:4285415/1‑70 (MQ=255)
                                              aTAAAACGCCGGGGCATTATTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGt              >  1:352960/1‑71 (MQ=255)
                                               tAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCtttt                              <  1:1732629/54‑1 (MQ=255)
                                               tAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGt              >  1:3391847/1‑70 (MQ=255)
                                               tAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTg             <  1:3941174/71‑1 (MQ=255)
                                                 aaaCGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAg                                   >  1:2416625/1‑47 (MQ=255)
                                                 aaaCGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGa            <  1:3595993/70‑1 (MQ=255)
                                                 aaaCGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGa            <  1:2525791/70‑1 (MQ=255)
                                                 aaaCGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGAt           >  1:340589/1‑71 (MQ=255)
                                                 aaaCGCCGGGGCATTGTTAGCAAATACGGCAACAACACCGCCTTCAg                                   >  1:299416/1‑47 (MQ=255)
                                                 aaaCGCCGGGGCATTATTAGCAAATACGGCAACcaca                                             >  1:3062823/1‑37 (MQ=255)
                                                 aaaCGCCGGGGCATTATTAGCAAATACGGCAACca                                               >  1:574671/1‑35 (MQ=255)
                                                 aaaCGCCGGGGCATTATTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCGGAACGGTTTGGTGa            >  1:4674596/1‑70 (MQ=255)
                                                   aCGCCGGGGCATTGTTAGCAAATACGGCAACCacac                                            >  1:5085359/1‑36 (MQ=255)
                                                   aCGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCg         <  1:4054864/71‑1 (MQ=255)
                                                   aCGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCg         <  1:3818375/71‑1 (MQ=255)
                                                      ccGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTAAGCTTTTGCCAGAACGGTTTGGTGATCgtgt      <  1:2941646/71‑1 (MQ=255)
                                                          ggCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCa  <  1:4777973/71‑1 (MQ=255)
                                                              ttGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAg                         <  1:4050938/44‑1 (MQ=255)
                                                              ttGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAg                         <  1:1342836/44‑1 (MQ=255)
                                                                |                                                                
TTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTACCA  >  W3110S.gb/4606025‑4606153

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: