Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,901,6010TA60.0% 83.6 / 157.8 125S480S (TCA→TCTrecBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base T (50/0);  new base A (75/0);  total (125/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 = 9.36e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 125 total aligned reads displayed.

GTTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATT  >  minE/1901557‑1901669
                                            |                                                                    
gtTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGa                                                                    >  1:1994370/1‑47 (MQ=255)
gtTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGa                                                                    >  1:1723038/1‑47 (MQ=255)
gtTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGa                                                                    >  1:255439/1‑47 (MQ=255)
gtTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGa                                                                    >  1:867162/1‑47 (MQ=255)
gtTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGa                                                                    >  1:14647/1‑47 (MQ=255)
gtTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTAttt                                             >  1:2791764/1‑70 (MQ=255)
gtTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTAttt                                             >  1:2657565/1‑70 (MQ=255)
gtTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTATTTc                                            >  1:2538746/1‑71 (MQ=255)
gtTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTATTTc                                            >  1:2081800/1‑71 (MQ=255)
gtTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTATTTc                                            >  1:2975130/1‑71 (MQ=255)
gtTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTATTTc                                            >  1:2242257/1‑71 (MQ=255)
gtTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTATTTc                                            >  1:1330003/1‑71 (MQ=255)
gtTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTATTTc                                            >  1:2436985/1‑71 (MQ=255)
gtTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTATTTc                                            >  1:1075603/1‑71 (MQ=255)
gtTTCACATTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTAttt                                             >  1:1057727/1‑70 (MQ=255)
 tttCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTATTTcg                                           >  1:3076548/1‑71 (MQ=255)
 tttCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTATTTcg                                           >  1:3010723/1‑71 (MQ=255)
 tttCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTATTTcg                                           >  1:1041938/1‑71 (MQ=255)
 tttCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTATTTc                                            >  1:2167467/1‑70 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGGACGCGTCATCAGTc                   >  1:1814648/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGt                    >  1:2672116/1‑70 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:2870708/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:2797618/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:2934310/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:2935748/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:2950145/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:2632883/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:2787920/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:2651548/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:3235677/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:583472/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:101086/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:955678/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:93477/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:831120/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:805321/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:709967/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:690493/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:676561/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:641706/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:3108289/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:574723/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:556129/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:496942/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:429955/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:3250604/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:3237284/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:961602/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:3208311/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:3165589/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:1590449/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:2062919/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:204382/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:1882213/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:1873522/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:1796973/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:1707321/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:1706350/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:1656889/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:1826499/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:1542115/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:1388430/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:1329430/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:1285164/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:1219705/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:1193037/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:1016640/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:2122144/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:2086390/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:2387899/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:216443/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:2444015/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:2191351/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:2427119/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:2410649/1‑71 (MQ=255)
                         aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCCTCATCAGt                    >  1:2195654/1‑70 (MQ=255)
                         aCGCCTCATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                   >  1:419518/1‑71 (MQ=255)
                                          gCTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:823335/1‑71 (MQ=255)
                                          gCTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:379493/1‑71 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGtc                          >  1:777197/1‑46 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:265262/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:2592959/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:723564/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:1280698/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:1093859/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:2689309/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:2694305/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:2701097/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:1691566/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:2365655/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:2330239/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:487710/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:1780431/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:384838/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:365796/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:2526524/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:1842206/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:3197843/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:2122866/1‑70 (MQ=255)
                                           cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt  >  1:2052255/1‑70 (MQ=255)
                                            |                                                                    
GTTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATT  >  minE/1901557‑1901669

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: