Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,901,6280GA59.5% 90.1 / 151.3 126R471R (CGC→CGTrecBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base G (48/3);  new base A (75/0);  total (123/3)
Fisher's exact test for biased strand distribution p-value = 6.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 130 total aligned reads displayed.

TTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTCACGCTGTTCACCATTCCTGGTGCGGAAC  >  minE/1901558‑1901698
                                                                      |                                                                      
tttCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTATTTcg                                                                        >  1:1041938/1‑71 (MQ=255)
tttCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTATTTcg                                                                        >  1:3076548/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGt                                                 >  1:2672116/1‑70 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:3108289/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:3208311/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:2998294/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:2950145/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:2934310/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:2910087/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:2870708/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:2787920/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:101086/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:2622392/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:2444015/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:2427119/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:2387899/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:2225044/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:2191351/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:583472/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:961602/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:690493/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:676561/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:663287/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:3165589/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:3235677/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:2651548/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:1707321/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:1016640/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:1193037/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:1219705/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:123869/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:1329430/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:1347456/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:1388430/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:1542115/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:1641963/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:1656889/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:1710464/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:1796973/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:1818474/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:1826499/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:1873522/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:1882213/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:2122144/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:2063658/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:2062919/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:1993919/1‑71 (MQ=255)
                        aCGCCTGATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCCTCATCAGt                                                 >  1:2195654/1‑70 (MQ=255)
                        aCGCCTAATTTTTGCCGGCAGATTTCACAGGAATAAACGGTATTTCACGAAACATGAACGCGTCATCAGTc                                                >  1:417378/1‑71 (MQ=255)
                                         gCTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:379493/1‑71 (MQ=255)
                                         gCTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:823335/1‑71 (MQ=255)
                                          cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:2741646/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:1280698/1‑70 (MQ=255)
                                          cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:384838/1‑70 (MQ=255)
                                          cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:723564/1‑70 (MQ=255)
                                          cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:732856/1‑70 (MQ=255)
                                          cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:3239803/1‑70 (MQ=255)
                                          cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:2694305/1‑70 (MQ=255)
                                          cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:3197843/1‑70 (MQ=255)
                                          cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:1691566/1‑70 (MQ=255)
                                          cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:2122866/1‑70 (MQ=255)
                                          cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:2052255/1‑70 (MQ=255)
                                          cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:265262/1‑70 (MQ=255)
                                          cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:1780431/1‑70 (MQ=255)
                                          cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:2526524/1‑70 (MQ=255)
                                          cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:2592959/1‑70 (MQ=255)
                                          cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:2701097/1‑70 (MQ=255)
                                          cTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTAtt                               >  1:1842206/1‑70 (MQ=255)
                                                tCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTc                              >  1:2767363/1‑65 (MQ=255)
                                                tCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTCa                             >  1:1202341/1‑66 (MQ=255)
                                                tCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTCa                             >  1:772046/1‑66 (MQ=255)
                                                tCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTCa                             >  1:1126628/1‑66 (MQ=255)
                                                tCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTCa                             >  1:275641/1‑66 (MQ=255)
                                                tCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTCa                             >  1:1389784/1‑66 (MQ=255)
                                                tCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTCa                             >  1:426741/1‑66 (MQ=255)
                                                tCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTCa                             >  1:3274251/1‑66 (MQ=255)
                                                tCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTCa                             >  1:154337/1‑66 (MQ=255)
                                                tCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTCa                             >  1:2427618/1‑66 (MQ=255)
                                                tCACTGGAATAAACGGAATTTCGCGAAACATGAAcc                                                           >  1:2856480/1‑35 (MQ=39)
                                                              ggTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCt                                   <  1:2673397/46‑1 (MQ=255)
                                                              ggTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCt                                   <  1:2984548/46‑1 (MQ=255)
                                                                tATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCt                                   <  1:2336098/44‑1 (MQ=255)
                                                                tATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTCACGCTGTTCACCATTCCTGGTg        >  1:1473524/1‑71 (MQ=255)
                                                                tATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTCACGCTGTTCACCATTCCTGGTg        >  1:732489/1‑71 (MQ=255)
                                                                tATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTCACGCTGTTCACCATTCCTGGTg        >  1:3175096/1‑71 (MQ=255)
                                                                tATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTCACGCTGTTCACCATTCCTGGTg        >  1:1070273/1‑71 (MQ=255)
                                                                tATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTCACGATGTTCACCATTCCTGGTg        >  1:2082053/1‑71 (MQ=255)
                                                                      gcgAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTCACGCTGTTCACCATTCCTGGTGCGGAAc  <  1:2157882/71‑1 (MQ=255)
                                                                      |                                                                      
TTTCACCTTTAAATACAAAACGTAACGCCTGATTTTTCCCGGCTGATTTCACTGGAATAAACGGTATTTCGCGAAACATGAACGCGTCATCAGTCTGGCTGAAAAGCTTATTCACGCTGTTCACCATTCCTGGTGCGGAAC  >  minE/1901558‑1901698

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: