Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,634,5970CA57.9% ‑5.1 / 17.0 19V97F (GTT→TTT) tfaQputative tail fiber assembly protein TfaQ
Reads supporting (aligned to +/- strand):  ref base C (2/6);  new base A (7/4);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCTTCCAGCAACGAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTT  >  NC_000913/1634389‑1634826
                                                                                                                                                                                                                |                                                                                                                                                                                                                                     
gCTTCCAGCAACGAGATTTCTTCCTCCGTTGCAATTTCCAGCTCTGCAGCATCCTGAAGCGGCGCAGTATGCTACCTGGCTAACTGCATAAGGTTGTTTTTTGTTTCTTCCGCCTCCAGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTTCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACAGAGTTCAGAAAAAAATAACGCGTc                                                                                                                                                                                             >  1:277571/1‑251 (MQ=11)
                 ttcgtcctccGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGGATCCTTCACCCACGCTTTGCCGTGCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGc                                                                                                                                                                                           <  1:222883/232‑1 (MQ=11)
                               cAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATagacggttaccccg                                                                                                                                                              >  1:142025/1‑238 (MQ=11)
                                            ngcagcaTCCTGAAGCGGCGCAATATGCTCACTGGCTAACTGCATCAGGCTTTTTTTTGTTTCTTCAGCATCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACACAGGCTGTGCCGTTCCCATTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGc                                                                                                                                                                                           >  1:63739/2‑209 (MQ=11)
                                                             gcgcAATCTGCTCACTGGCTCCCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCCGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGaaccc                                                                                                                                >  1:272420/1‑247 (MQ=11)
                                                              cgcAATATGCTCACTGGCTACCTGCATCAGGCTTTTTTTTGATTCTTCCGCCTCCCGGATCCGGACTAGTTTTTCTGCTTCCGTATCCTTCACCCCGGCTGCGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGCAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATgcct                                                                                                                               >  1:134495/1‑248 (MQ=11)
                                                                      gCTCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACt                                                                                                                        <  1:283827/250‑1 (MQ=11)
                                                                        tCACTGGCTACCTGCATCAGGCTTTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAAGTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGc                                                                                                                                         <  1:183297/231‑1 (MQ=11)
                                                                                               tttttttGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTGTCCTTCACCCAGGCTGTGCCGGTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAAGTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGTAATCCACGTCATAGACGGGTTTGCGCCGATGGTCTTCAACTAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAagccagcc                                                                                               <  1:312865/249‑1 (MQ=11)
                                                                                                                     cGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCACTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATc                                                                        >  1:279277/1‑251 (MQ=11)
                                                                                                                                   tttttCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCTGTAACGTACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAAcag                                                                                                                                                                      <  1:312597/143‑3 (MQ=12)
                                                                                                                                       tCTGCTTCCGTATCCTTCACCCAGTCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTTCCCCGATGATCTTCATCGAGATGCGTCGAGGACTCAGCACTGTTGTAATCAACTACTAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCa                                                       <  1:307472/250‑1 (MQ=11)
                                                                                                                                                              ggCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAatatatn                               <  1:79650/251‑2 (MQ=14)
                                                                                                                                                                  gtgCCGTTCCTCTTCTGAAACTCCCCTTCCGGTGTTAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAAGAAATAACGCTTCGCCGGAAGCCACGTCAGAAGCCGTTTTTCCCCTATGATCCTCAACGAGATGCCACGATGACTCATCACTGTTGAGAACAGCCACGAAGCCAGCCCCAATATCGGTCGTTGCAATATCGGTACGGTTTGCTGGCAGACCTGTATGAGGCGGAatatat                                <  1:167149/246‑1 (MQ=11)
                                                                                                                                                                                          cctcccGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCCCGTCATAGACGGTTTTAcccc                                                                                                                                                             >  1:120024/5‑97 (MQ=11)
                                                                                                                                                                                             tCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATTACTCATCTCTGTTGAAAACAGCCGCGAAGCCAGCCGGAATATCTGGCGGTGTAGTATCGGTATTGTTAGCAAGCATACCTGTATGAGGCTGAATATATGCATCACCTTCGCCAATAAATTCATTAGtt  <  1:162797/249‑1 (MQ=11)
                                                                                                                                                                                             nccGGTGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGAAAACAGCCACAAag                                                                                                     >  1:63017/2‑150 (MQ=11)
                                                                                                                                                                                                      aaCCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCtn                                                          <  1:80517/184‑2 (MQ=17)
                                                                                                                                                                                                         cAGGTAACATTTTCCCCTACCGGACCGAGTTCTGAAATAACTAACGCGTCGCCGGAAGCCACGTc                                                                                                                                                                              >  1:116770/1‑65 (MQ=255)
                                                                                                                                                                                                                |                                                                                                                                                                                                                                     
GCTTCCAGCAACGAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTT  >  NC_000913/1634389‑1634826

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑