Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,618,668 A→G 6.9% intergenic (+76/‑74) yneK → / → ydeA hypothetical protein/predicted arabinose transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,618,6680AG 6.9% 163.5 / 3.8 58intergenic (+76/‑74)yneK/ydeAhypothetical protein/predicted arabinose transporter
Reads supporting (aligned to +/- strand):  ref base A (31/23);  new base G (2/2);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGAT  >  W3110S.gb/1618602‑1618733
                                                                  |                                                                 
ttACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATg                                                                 <  1:1288823/69‑1 (MQ=255)
ttACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGc                                                                <  1:410032/70‑1 (MQ=255)
ttACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGc                                                                <  1:1968352/70‑1 (MQ=255)
   ccGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTg                                                            <  1:3263673/71‑1 (MQ=255)
   ccGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTg                                                            <  1:31870/71‑1 (MQ=255)
    cGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGt                                                           <  1:259330/71‑1 (MQ=255)
    cGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGt                                                           <  1:1716717/71‑1 (MQ=255)
     gACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTa                                                          <  1:2890081/71‑1 (MQ=255)
            tGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTcac                                                     <  1:2254811/69‑1 (MQ=255)
                     aaaGAGTCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCg                                          >  1:2871853/1‑71 (MQ=255)
                     aaaGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGTTCACAAATGCATTGTATTCACATCATTAACCg                                          >  1:1550063/1‑71 (MQ=255)
                     aaaGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCg                                          >  1:3410619/1‑71 (MQ=255)
                     aaaGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCg                                          >  1:3261226/1‑71 (MQ=255)
                     aaaGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCg                                          >  1:3446625/1‑71 (MQ=255)
                     aaaGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCg                                          >  1:2306154/1‑71 (MQ=255)
                     aaaGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCg                                          >  1:1143875/1‑71 (MQ=255)
                     aaaGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCg                                          >  1:1638023/1‑71 (MQ=255)
                     aaaGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCg                                          >  1:1519073/1‑71 (MQ=255)
                       aGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGtt                                        >  1:169126/1‑71 (MQ=255)
                             aGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAAc                                            >  1:1769013/1‑61 (MQ=255)
                                   agagAAACGGAAAACGCGATCCAGATCACAAGTGCATTGTATTCACATCATTAACCGTTTTAAGATCAttt                            >  1:3011263/1‑71 (MQ=255)
                                   agagAAACGGAAAACGCGATCCAGATCACAAGTGCATTGTATTCACATCATTAACCGTTTTAAGATCAttt                            >  1:2722056/1‑71 (MQ=255)
                                   agagAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCAtt                             <  1:3442091/70‑1 (MQ=255)
                                   agagAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCAtt                             <  1:1520612/70‑1 (MQ=255)
                                   agagAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCAtt                             <  1:3289888/70‑1 (MQ=255)
                                       aaaCGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGAt                                 <  1:2456957/62‑1 (MQ=255)
                                       aaaCGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGAt                                 <  1:2538464/62‑1 (MQ=255)
                                       aaaCGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGAt                                 <  1:3726109/62‑1 (MQ=255)
                                       aaaCGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTtcat                         >  1:3328154/1‑70 (MQ=255)
                                         aCGGAAAACGCGATCCAGATCACAAGTGCATTGTATTCATTCATTAACCGTTTTAAGATCATTTCATCACt                     <  1:2438339/71‑1 (MQ=255)
                                            gAAAACGCGATCCAGATCACAAGTGCATTGTATTCACATCATTAACCGTTTTAAGATCAt                              <  1:1108131/60‑1 (MQ=255)
                                               aaCGCGATCCAGATCACAAATGCATTGTATTCAcatca                                                 >  1:3271474/1‑38 (MQ=255)
                                               aaCGCGATCCAGATCACAAATGCATTGTATTCAcatca                                                 >  1:3256183/1‑38 (MQ=255)
                                               aaCGCGATCCAGATCACAAATGCATTGTATTCAcatca                                                 >  1:3350596/1‑38 (MQ=255)
                                               aaCGCGATCCAGATCACAAATGCATTGTATTCAcatca                                                 >  1:1866759/1‑38 (MQ=255)
                                               aaCGCGATCCAGATCACAAATGCATTGTATTCAcatca                                                 >  1:2398671/1‑38 (MQ=255)
                                               aaCGCGATCCAGATCACAAATGCATTGTATTCAcatca                                                 >  1:2112158/1‑38 (MQ=255)
                                                aCGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCg               >  1:610516/1‑71 (MQ=255)
                                                aCGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCg               >  1:1472304/1‑71 (MQ=255)
                                                aCGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCg               >  1:3577145/1‑71 (MQ=255)
                                                aCGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATATCATCACTTTTTCg               >  1:1614523/1‑71 (MQ=255)
                                                 cgcgATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCg               >  1:1507176/1‑70 (MQ=255)
                                                 cgcgATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGc              <  1:3502501/71‑1 (MQ=255)
                                                 cgcgATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGc              <  1:529167/71‑1 (MQ=255)
                                                 cgcgATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGc              <  1:280953/71‑1 (MQ=255)
                                                 cgcgATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGc              <  1:2902859/71‑1 (MQ=255)
                                                 cgcgATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGc              <  1:1009414/71‑1 (MQ=255)
                                                       ccAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCa                               <  1:3205228/48‑1 (MQ=255)
                                                       ccAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCa                               <  1:1736717/48‑1 (MQ=255)
                                                       ccAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCa                               <  1:729395/48‑1 (MQ=255)
                                                         aGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCAtt                             >  1:2656223/1‑48 (MQ=255)
                                                         aGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCAtt                             >  1:1209890/1‑48 (MQ=255)
                                                            tCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGa   >  1:3597256/1‑71 (MQ=255)
                                                            tCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGa   >  1:487361/1‑71 (MQ=255)
                                                             cacaAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGa   >  1:3539243/1‑70 (MQ=255)
                                                              acaAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGAt  >  1:3026182/1‑70 (MQ=255)
                                                              acaAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGAt  >  1:242114/1‑70 (MQ=255)
                                                              acaAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGAt  >  1:561381/1‑70 (MQ=255)
                                                                  |                                                                 
TTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGAT  >  W3110S.gb/1618602‑1618733

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