Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913566,2450GA53.8% ‑5.1 / 12.3 13intergenic (‑267/‑439)intD/renDputative integrase/protein RenD
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base A (2/5);  total (6/7)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTCACCTGGACAGTATGCAGCGGAGACTGAAGCGCCGCAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGT  >  NC_000913/566009‑566451
                                                                                                                                                                                                                                            |                                                                                                                                                                                                              
cTCACCTGGACAGTATGCAGCGGAGACTGAAGCGCCGCAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCCTTGCTCCCCAAATACAAAACCAATTTCAGCCCGTGCCTCGTCCATTTTTTTGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGtt                                                                                                                                                                                                  >  1:319264/1‑251 (MQ=255)
                      nagaCTGAAGCGCCGCAATGCAGGCTTGCCGTGTTGTGACGTACCCGGATATTATCGTGTGGATGCGTCATCGACATTGCTCCCCAACTACAAAACCAATTTCAGCCAGTGCCTCGCCCATTTTTTCGATGAAATCCGGACCCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACAAGGAATCtttt                                                                                                                                                                            >  1:31084/2‑251 (MQ=255)
                                               cgaaagtgattgcgcctACCCGGATATTATCGTGAGGATGCGTCATCGCAATTGCTCCACAAATACAAAACCAATTTAAGCAAGTGCCTCGTCCATTTTTTAGATGAACTCCGGCACGATCTAGTAAAAACTCGACATGTACTTTTAATACCGCTCAATCACGCCATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGtt                                                                                                                                                                                                  >  1:242095/17‑204 (MQ=255)
                                                                                  ggATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGGTGAACTCCGGCACTATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAACGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTCGCAAAGTACCAGGCATGTTTTCTCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTAtggc                                                                                                                       <  1:290849/244‑4 (MQ=255)
                                                                                       cGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGGGCCTCGTCCATGTGTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTGTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTgn                                                                                                                                                                                                <  1:35079/166‑2 (MQ=255)
                                                                                           aTCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCAATTTTTTCGATGAACTACGGCACGATCTCGTCAAAACTCGACATGTACTTTTCATCCCGCTCAATCCCGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGc                                                                                                       >  1:300101/1‑251 (MQ=255)
                                                                                                                     ccAATTTCAGCCAGTGCCTCGTCCATTTTATCGATGAACTCAGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTAAATCACGGCATAATGCAGGCCTTCACGCTTCGTGCGCGGGTCATAGATGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTGCTGCACCTGGGACATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAACTCCCGGGAGGt                                                                             >  1:100358/1‑251 (MQ=255)
                                                                                                                       aaTTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCCTGCAATTTTCATTCCGCTCAATCACGCCATAATGCAGGCATTCACGCTTCATGCGCGGGTCATAGTTGGCAACGTCCCCGGCATCTTTTCGCGTCACCCACCTGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATg                                                                                           >  1:186575/1‑235 (MQ=255)
                                                                                                                                    gCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTAtgg                                                                                                                        <  1:103940/193‑3 (MQ=255)
                                                                                                                                                    cGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTGTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATc                                               <  1:132252/250‑1 (MQ=255)
                                                                                                                                                      aTGAACTCCGGCACGATCTCGTCAAAACTCGCCATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGtn                                                                                                                                                                      <  1:85010/129‑2 (MQ=255)
                                                                                                                                                                                                 cATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtggt  <  1:25037/250‑1 (MQ=255)
                                                                                                                                                                                                          cAATCTCGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTcg                                                                                                                                                                          <  1:186822/73‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                              
CTCACCTGGACAGTATGCAGCGGAGACTGAAGCGCCGCAATGCAGGCTTGCCGTGTTGTGAGGTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAATTTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACCATCTCGTCAAAACTCGCCATGTACTTTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGT  >  NC_000913/566009‑566451

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