Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,634,7780TG60.0% ‑6.4 / 23.1 25T36T (ACA→ACC) tfaQputative tail fiber assembly protein TfaQ
Reads supporting (aligned to +/- strand):  ref base T (8/2);  new base G (9/6);  total (17/8)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTCTGAATGCCATT  >  NC_000913/1634535‑1634886
                                                                                                                                                                                                                                                    |                                                                                                            
aTCCTTCACCCAGGCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATg                                                                                                         >  2:108261/1‑249 (MQ=11)
            ggCTGTGCCGTTCCACTTCTGATATTCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATg                                                                                           >  2:73111/1‑251 (MQ=11)
               tgtgCCGTTCCACTTCTCAAACTCCCCTTCCGGCGATAACCCGGAAACATTTTCCGGTAACGGACCGAGTTCAGAAATCAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCAGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGTTGCAATATCGGTCCTGCTTGCTGGCAGACCTGTATGAGGCGGAATATATGCAt                                                                                        >  1:351066/1‑251 (MQ=11)
                        ccACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCa                                                                                <  2:71250/250‑1 (MQ=31)
                             tctgatattCCCCTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGTCGGAATATATTCGTCGCCTTCACCAAt                                                                           <  1:82552/242‑1 (MQ=11)
                             tCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAATCCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAAt                                                                           >  1:373253/1‑250 (MQ=14)
                                  aaCTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGCAATAAATAACGCGTCGACGGAAGCCACGTCCTAAACCGTTTTAACCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGCGCAATATCGGTACTGTTTGCTGGCAGAACTGTATTCGGCGGAATATATGCCTCACCTTCACCAATAAATTc                                                                     >  2:146645/1‑251 (MQ=11)
                                     tacccTCCCGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAATCCACGTCATAGACGGTTTTACACCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTGTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCAt                                                                   <  2:76618/248‑1 (MQ=11)
                                      ccccTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCGCTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGATGCGGAATAGATGCATCACCTGCACCAATAAATTCAt                                                                   <  2:233735/249‑1 (MQ=9)
                                        cctcccGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGt                                                               >  2:201437/5‑251 (MQ=11)
                                        cctcccGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGt                                                               >  2:70946/5‑251 (MQ=11)
                                         ctcccGGCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGGATGAGGGGGAATATATGCGTCACCTGCACCAATAAATTCATTAGtt                                                              >  2:294115/4‑251 (MQ=11)
                                              ggCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGc                                                         >  2:64722/1‑251 (MQ=31)
                                               gCGATAACCAGGTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGTGACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGcc                                                        <  2:240020/251‑1 (MQ=11)
                                                        aGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGAtt                                                >  2:196474/1‑250 (MQ=31)
                                                         ggTAAAATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATtat                                              <  2:117849/251‑1 (MQ=11)
                                                              aaTTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTCTGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAAtt                                         <  2:364233/250‑1 (MQ=11)
                                                                      ggTAACGGACCGAGTTCAGAAATAAATAACGAGTCGCCGGAAGCAACGACATAGACGGTATTACCCCGATGGTCTTCCACGAGCTGCCCAGCAGCCTCATCACTGTTGACAAAC‑GCCACAAAGCCCGCAGGAAAATCTGGCGGTGCAATATAGGTACTGTTTGCAGGCAGACCGGTCTGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGt                                 >  1:221240/1‑251 (MQ=2)
                                                                        tAACGGAACGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCCTTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTcc                               >  2:217669/1‑251 (MQ=11)
                                                                           cGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCAACGTCATAAACAGTTTTACCCCGATGATCTTCAACGAGATGCAACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTCCTGTTTGCTGGAAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCCTTAGTTCCAGCCAGCAGATTATAAATTTTTATGGTCCGTg                            >  1:59684/1‑251 (MQ=11)
                                                                                  aGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTc                     <  2:320875/251‑1 (MQ=11)
                                                                                        gAAATAAATAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCAt               >  1:22930/1‑251 (MQ=11)
                                                                                         aaataaatAACGCGTCGCCGGAAGCCACGTCATAGACGGTTTTACCCCGATGGTCTTCAACGAGATGCCACGATGCCTCATCACTGTTGA‑AAACAGCCACAAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCAGGCAGACCGGTATGAGGCGGAATATATGCGTCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCAtt              >  1:79928/1‑251 (MQ=11)
                                                                                          aataaataACGCGTCGCCGGAAGCCACGTCATACACCGTTTTACACCGATGATCTTCAACGAGATGCCACGATAACTCATCCCTGTTGA‑AAACAGCCACGAAGCAAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTc             >  1:210167/1‑251 (MQ=9)
                                                                                                      gTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTATTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTCTGAATGCCAtt  >  2:256949/1‑250 (MQ=14)
                                                                                                                                                                                                                                                    |                                                                                                            
ATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGA‑AAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCTGTATGAGGCGGAATATATGCATCACCTTCACCAATAAATTCATTAGTTCCGGCCAGCAGATTATAAATTTTTATGGTCCGTGGTTGTTCACTCATTCTGAATGCCATT  >  NC_000913/1634535‑1634886

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑