Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,633,9840CA55.6% ‑5.4 / 11.2 18pseudogene (408/528 nt)nohQpseudogene, Qin prophage; predicted packaging protein;Phage or Prophage Related; DNA packaging protein NU1 homolog from lambdoid prophage Qin
Reads supporting (aligned to +/- strand):  ref base C (0/8);  new base A (8/2);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 1.05e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCC  >  NC_000913/1633893‑1634065
                                                                                           |                                                                                 
cacccAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAgg                                                                          <  2:4442134/98‑1 (MQ=1)
                       attgCTTTGATTATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGc                                                   <  2:6690353/98‑1 (MQ=0)
                               tattatatCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCAc                                           <  1:8494260/97‑1 (MQ=1)
                                 tgctatCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCt                                         <  2:5477537/98‑1 (MQ=1)
                                      tcacGTTTCATGAAATCAACATGTCGGTATTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGAt                                    <  2:3584923/98‑1 (MQ=0)
                                           cttgaGGAAATCAATATGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCcgcg                               >  1:12761905/5‑101 (MQ=0)
                                           cttgaGGAAATCAATATGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCcgcg                               >  1:2986271/5‑101 (MQ=0)
                                           cttgaGGAAATCAATATGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCcgcg                               >  1:4952798/5‑101 (MQ=0)
                                           cttgaGGAAATCAATATGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCcgcg                               >  1:5658642/5‑101 (MQ=0)
                                           cttgaGGAAATCAATATGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGCCAATCcgcg                               >  1:12842701/5‑101 (MQ=0)
                                            ttaaGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACACAGCTGCACCCACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGa                              <  2:4507280/98‑1 (MQ=1)
                                                       aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGaacg                   >  1:12358707/4‑99 (MQ=0)
                                                       aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTATTCGGCCATCCGCGACAGCACGAAAg                   >  1:13918625/4‑101 (MQ=0)
                                                       aataTGTCGGTTTGCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGaacg                   <  1:5003389/98‑3 (MQ=0)
                                                             tctgTTTTCCAGTTCCGTAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGa             <  2:6894158/98‑1 (MQ=0)
                                                              cgattttCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGaa            <  1:6045292/98‑1 (MQ=1)
                                                                       agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc   >  1:13605162/4‑101 (MQ=0)
                                                                       agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc   <  1:11380328/98‑1 (MQ=0)
                                                                        gtgccGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTcc  <  2:5370333/98‑1 (MQ=1)
                                                                                           |                                                                                 
CATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCC  >  NC_000913/1633893‑1634065

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