Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,793,663 G→A 17.0% V24V (GTC→GTT ydiV ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,793,6630GA17.0% 131.7 / 18.8 53V24V (GTC→GTTydiVconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (13/31);  new base A (4/5);  total (17/36)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

ATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTT  >  W3110S.gb/1793621‑1793731
                                          |                                                                    
aTCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGa                                          <  1:910437/71‑1 (MQ=255)
aTCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGGTGTCTCTTATCGGCAGGa                                          <  1:1776917/71‑1 (MQ=255)
 tCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGa                                          <  1:9990/70‑1 (MQ=255)
 tCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGa                                          <  1:869157/70‑1 (MQ=255)
 tCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGa                                          <  1:1917015/70‑1 (MQ=255)
 tCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGa                                          <  1:1773339/70‑1 (MQ=255)
 tCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGa                                          <  1:344047/70‑1 (MQ=255)
 tCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGa                                          <  1:1707096/70‑1 (MQ=255)
 tCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGa                                          <  1:1574730/70‑1 (MQ=255)
 tCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGa                                          <  1:1529404/70‑1 (MQ=255)
 tCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGa                                          <  1:1351231/70‑1 (MQ=255)
   ttCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGaaa                                        >  1:204330/1‑70 (MQ=255)
   ttCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGaaa                                        >  1:1614570/1‑70 (MQ=255)
   ttCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGaaa                                        >  1:1746803/1‑70 (MQ=255)
     cACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGAACTTGCTGGTTGTCTCTTATCGGCAGGAAATAg                                     >  1:311349/1‑71 (MQ=255)
     cACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGAACTTGCTGGTTGTCTCTTATCGGCAGGAAATAg                                     >  1:1659544/1‑71 (MQ=255)
     cACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGAACTTGCTGGTTGTCTCTTATCGGCAGGAAATAg                                     >  1:335641/1‑71 (MQ=255)
     cACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGAACTTGCTGGTTGTCTCTTATCGGCAGGAAATAg                                     >  1:1150323/1‑71 (MQ=255)
       cTGGAGAAATGGGTAATTAACTCAACGCCGATAAGAACTTGCTg                                                              <  1:1847770/44‑1 (MQ=38)
       cTGGAGAAATGGGTAATTAACTCAACGCCGACAAGAACTTGCTg                                                              <  1:476601/44‑1 (MQ=255)
       cTGGAGAAATGGGTAATTAACTCAACGCCGACAAGAACTTGCTg                                                              <  1:757755/44‑1 (MQ=255)
       cTGGAGAAATGGGTAATTAACTCAACGCCGACAAGAACTTGCTg                                                              <  1:1638643/44‑1 (MQ=255)
       cTGGAGAAATGGGTAATTAACTCAACGCCGACAAGAACTTGCTg                                                              <  1:1432676/44‑1 (MQ=255)
               aTGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAAt                           >  1:854735/1‑71 (MQ=255)
               aTGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAAt                           >  1:602124/1‑71 (MQ=255)
               aTGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAAt                           >  1:1823115/1‑71 (MQ=255)
               aTGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAAt                           >  1:319542/1‑71 (MQ=255)
               aTGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAAt                           >  1:1589115/1‑71 (MQ=255)
               aTGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAAt                           >  1:310475/1‑71 (MQ=255)
               aTGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAAt                           >  1:1117525/1‑71 (MQ=255)
               aTGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAAt                           >  1:908203/1‑71 (MQ=255)
                      aTTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATaaa                     <  1:476102/70‑1 (MQ=255)
                      aTTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATaaa                     <  1:849297/70‑1 (MQ=255)
                      aTTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATaaa                     <  1:1038727/70‑1 (MQ=255)
                      aTTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATaaa                     <  1:1749853/70‑1 (MQ=255)
                      aTTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATaaa                     <  1:1265028/70‑1 (MQ=255)
                       ttAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGa                   <  1:1887340/71‑1 (MQ=255)
                       ttAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGa                   <  1:1876665/71‑1 (MQ=255)
                       ttAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGa                   <  1:1783297/71‑1 (MQ=255)
                       ttAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGa                   <  1:348845/71‑1 (MQ=255)
                       ttAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGa                   <  1:351176/71‑1 (MQ=255)
                       ttAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGa                   <  1:1420532/71‑1 (MQ=255)
                       ttAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGa                   <  1:566925/71‑1 (MQ=255)
                         aaCTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGAt                  >  1:1727320/1‑70 (MQ=255)
                         aaCTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGAt                  >  1:1332784/1‑70 (MQ=255)
                            tCAACGCCGACAAGGACTTGCTGGGTGTCTCTTATCGGCAGGa                                          <  1:935892/43‑1 (MQ=255)
                                        aGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCtt  <  1:1861681/71‑1 (MQ=255)
                                        aGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCtt  <  1:732872/71‑1 (MQ=255)
                                        aGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCtt  <  1:277752/71‑1 (MQ=255)
                                        aGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCtt  <  1:264310/71‑1 (MQ=255)
                                        aGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCtt  <  1:138714/71‑1 (MQ=255)
                                        aGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCtt  <  1:1127187/71‑1 (MQ=255)
                                        aGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCtt  <  1:1063053/71‑1 (MQ=255)
                                          |                                                                    
ATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTT  >  W3110S.gb/1793621‑1793731

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: