Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb36,3850AG73.1% 77.4 / 30.4 52F485F (TTT→TTCcaiCpredicted crotonobetaine CoA ligase:carnitine CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (14/0);  new base G (38/0);  total (52/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTG  >  W3110S.gb/36324‑36436
                                                             |                                                   
ttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaaa                                                     >  1:1076845/1‑62 (MQ=255)
ttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaaa                                                     >  1:82851/1‑62 (MQ=255)
ttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaaa                                                     >  1:779898/1‑62 (MQ=255)
ttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaaa                                                     >  1:512646/1‑62 (MQ=255)
ttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaaa                                                     >  1:322912/1‑62 (MQ=255)
ttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaaa                                                     >  1:2897304/1‑62 (MQ=255)
ttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaaa                                                     >  1:2777424/1‑62 (MQ=255)
ttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaaa                                                     >  1:2767212/1‑62 (MQ=255)
ttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaaa                                                     >  1:2387148/1‑62 (MQ=255)
ttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaaa                                                     >  1:2090216/1‑62 (MQ=255)
ttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaaa                                                     >  1:1994454/1‑62 (MQ=255)
ttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaaa                                                     >  1:1888159/1‑62 (MQ=255)
ttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaaa                                                     >  1:1850547/1‑62 (MQ=255)
ttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaaa                                                     >  1:1768569/1‑62 (MQ=255)
ttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaaa                                                     >  1:1230921/1‑62 (MQ=255)
ttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaaa                                                     >  1:1574467/1‑62 (MQ=255)
ttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGAGAAGCGGaaaa                                                    >  1:2712599/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAAtt                                                 >  1:1683625/1‑51 (MQ=37)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCtt                                           >  1:2822151/1‑57 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:2663633/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:576343/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:2693807/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:2704337/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:2725915/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:1444201/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:1388225/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:2802156/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:2804155/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:753160/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:710857/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:133231/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:707867/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:384785/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:697895/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:432426/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:4673/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:697402/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:75740/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:1710758/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:1680754/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:1670311/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:1990423/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:823014/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:206821/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:2077663/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:790526/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:1129408/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:2237087/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:2323084/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:2330280/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:1526521/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:2504066/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                      >  1:2613506/1‑62 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTACACTTCGCt                                       >  1:1529667/1‑61 (MQ=255)
               aaGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATACCTCTTCGCTc                                      >  1:2047356/1‑62 (MQ=255)
                            aaaTTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                         >  1:2104418/1‑62 (MQ=255)
                                                  cAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCttt   >  1:802930/1‑62 (MQ=255)
                                                  cAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCttt   >  1:1074235/1‑62 (MQ=255)
                                                  cAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCttt   >  1:501008/1‑62 (MQ=255)
                                                  cAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCttt   >  1:1289977/1‑62 (MQ=255)
                                                  cAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCttt   >  1:2853405/1‑62 (MQ=255)
                                                  cAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCttt   >  1:2841133/1‑62 (MQ=255)
                                                  cAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCttt   >  1:2554089/1‑62 (MQ=255)
                                                  cAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCttt   >  1:2210863/1‑62 (MQ=255)
                                                  cAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCttt   >  1:2090086/1‑62 (MQ=255)
                                                  cAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCttt   >  1:1998806/1‑62 (MQ=255)
                                                  cAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCttt   >  1:1557126/1‑62 (MQ=255)
                                                  cAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCttt   >  1:1694404/1‑62 (MQ=255)
                                                  cAGAAGCGGAAAAATTCCTCTTCGCTAATGTTTCACCTTCATTCAGCACCACAAATGCTTTg  >  1:409777/1‑62 (MQ=255)
                                                             |                                                   
TTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTG  >  W3110S.gb/36324‑36436

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