Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 344,381 G→A 32.0% intergenic (+166/‑247) yahL → / → yahM hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb344,3810GA32.0% 65.5 / 42.1 50intergenic (+166/‑247)yahL/yahMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (15/19);  new base A (12/4);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 6.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGA  >  W3110S.gb/344314‑344445
                                                                   |                                                                
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGttt                                                               >  1:317264/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGttt                                                               >  1:1745596/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGttt                                                               >  1:1664966/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGttt                                                               >  1:1623443/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGttt                                                               >  1:1778448/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGttt                                                               >  1:180121/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGttt                                                               >  1:412569/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGttt                                                               >  1:682428/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGttt                                                               >  1:116913/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGttt                                                               >  1:740524/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGttt                                                               >  1:1108178/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                               >  1:1914765/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                               >  1:1041384/1‑71 (MQ=255)
                          ttgtGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATg                                     <  1:1654053/71‑1 (MQ=255)
                          ttgtGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATg                                     <  1:1754224/71‑1 (MQ=255)
                               aTCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGAtct                                 <  1:136254/70‑1 (MQ=255)
                               aTCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGAtct                                 <  1:1382098/70‑1 (MQ=255)
                               aTCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGAtct                                 <  1:652242/70‑1 (MQ=255)
                               aTCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGAtct                                 <  1:1842109/70‑1 (MQ=255)
                               aTCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGAtct                                 <  1:1550241/70‑1 (MQ=255)
                                 cGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATg                                     <  1:54322/64‑1 (MQ=255)
                                     gACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAg                          <  1:1794840/71‑1 (MQ=255)
                                     gACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAg                          <  1:154015/71‑1 (MQ=255)
                                     gACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAg                          <  1:1435510/71‑1 (MQ=255)
                                         cTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                      <  1:1854746/71‑1 (MQ=255)
                                         cTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                      <  1:1902589/71‑1 (MQ=255)
                                         cTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                      <  1:1501711/71‑1 (MQ=255)
                                         cTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                      <  1:1939281/71‑1 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                    >  1:97103/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                    >  1:892864/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                    >  1:1208123/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                    >  1:332542/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                    >  1:570827/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                    >  1:168979/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                    >  1:743015/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                    >  1:1162897/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                    >  1:293634/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                    >  1:508368/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                    >  1:42825/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                    >  1:1538573/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                    >  1:1535920/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATCCAGATGtt                     >  1:757878/1‑70 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                <  1:1272674/69‑1 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                <  1:120130/69‑1 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                <  1:350767/69‑1 (MQ=255)
                                                    aaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAAtgt           <  1:787063/71‑1 (MQ=255)
                                                             gtagtGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGa  <  1:784693/71‑1 (MQ=255)
                                                             gtagtGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGa  <  1:1728463/71‑1 (MQ=255)
                                                              tagtGGTTTCTCAGCCATGCGTTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGa  <  1:1002474/70‑1 (MQ=255)
                                                              tagtGGTTTCGCAGCCATGCGGTGTATAAAAAATGAtct                                 <  1:1493408/39‑1 (MQ=255)
                                                                   |                                                                
CCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGA  >  W3110S.gb/344314‑344445

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: