Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 344,381 G→A 100% 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,3810GA94.4% 150.2 / 1.0 53intergenic (+166/‑247)yahL/yahMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base A (30/20);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACAT  >  W3110S.gb/344314‑344448
                                                                   |                                                                   
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                  >  1:3531827/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                  >  1:2050509/1‑71 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                  >  1:1583224/1‑71 (MQ=255)
 cATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                  <  1:1889255/70‑1 (MQ=255)
 cATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                  <  1:3088437/70‑1 (MQ=255)
 cATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                  <  1:536142/70‑1 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAg                                                             >  1:2737029/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAg                                                             >  1:2721928/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAg                                                             >  1:1185780/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAg                                                             >  1:491988/1‑71 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAg                                                             >  1:1204740/1‑71 (MQ=255)
         ttttCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCAt                                                         >  1:2874927/1‑71 (MQ=255)
         ttttCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCAt                                                         >  1:1733899/1‑71 (MQ=255)
            tCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGc                                                       <  1:3556457/70‑1 (MQ=255)
            tCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGc                                                       <  1:2807673/70‑1 (MQ=255)
              gCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCATTAGTGATTTCGCAGCCATGCgg                                                     <  1:1758145/70‑1 (MQ=255)
                       ttgttgTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATaaaa                                            <  1:3070586/70‑1 (MQ=255)
                       ttgttgTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATaaaa                                            >  1:1473774/1‑70 (MQ=255)
                       ttgttgTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATaaaa                                            <  1:1771853/70‑1 (MQ=255)
                        tgttgtGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGGAGCCATGCGGTGTATaaaaa                                           <  1:3667125/70‑1 (MQ=255)
                        tgttgtGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATaaaaaa                                          >  1:1618931/1‑71 (MQ=255)
                        tgttgtGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATaaaaa                                           <  1:347167/70‑1 (MQ=255)
                        tgttgtGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATaaaaa                                           <  1:1064234/70‑1 (MQ=255)
                                  ggCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTTTATAAAAAATGATCTCAt                                 <  1:3176510/70‑1 (MQ=255)
                                  ggCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCAt                                 <  1:361600/70‑1 (MQ=255)
                                  ggCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCAt                                 <  1:3628710/70‑1 (MQ=255)
                                    cGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCa                              <  1:3401057/71‑1 (MQ=255)
                                    cGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCa                              <  1:91143/71‑1 (MQ=255)
                                     gACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAg                             <  1:2663945/71‑1 (MQ=255)
                                         cTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATg                                        <  1:722818/56‑1 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAg                             <  1:3140916/65‑1 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                       >  1:3191698/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                       >  1:401398/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                       >  1:2687533/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                       >  1:1139418/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                       >  1:2542739/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                       >  1:1244021/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                       >  1:2356621/1‑71 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGttt                       >  1:1657003/1‑71 (MQ=255)
                                            cGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGAt                                      >  1:1494214/1‑55 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                         >  1:1836613/1‑63 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                         >  1:722448/1‑63 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                         >  1:2567014/1‑63 (MQ=255)
                                                 aaaaaaGCACCAGTAGTGATTTCGCAGCCATGCGCTGTATAAAAAATGATCTCATGCAGATGTTTTGTGaa                 >  1:1300592/1‑71 (MQ=255)
                                                    aaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAAtgt              <  1:470648/71‑1 (MQ=255)
                                                        cACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTgttg          >  1:2969513/1‑71 (MQ=255)
                                                        cACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTgttg          >  1:2084346/1‑71 (MQ=255)
                                                        cACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTgttg          >  1:2478786/1‑71 (MQ=255)
                                                        cACCAGTAGTGATTTCGCAGACATGCGGTGTATAAAAAATGAtctc                                   >  1:1849809/1‑46 (MQ=255)
                                                                gtGATTTCGCAGCCATGCGGTGTATAAAAAATGAt                                      >  1:1399795/1‑35 (MQ=255)
                                                                gtGATTTCGCAGCCATGCGGTGTATAAAAAATGAt                                      >  1:2373016/1‑35 (MQ=255)
                                                                gtGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACAt  >  1:3534386/1‑71 (MQ=255)
                                                                gtGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACAt  >  1:422620/1‑71 (MQ=255)
                                                                   |                                                                   
CCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACAT  >  W3110S.gb/344314‑344448

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