Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,633,9750AG52.9% ‑5.1 / 12.4 17pseudogene (417/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 A (0/8);  new base G (8/1);  total (8/9)
Fisher's exact test for biased strand distribution p-value = 4.11e-04
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.

CAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCC  >  NC_000913/1633876‑1634065
                                                                                                   |                                                                                          
cacccccGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAAc                                                                                           <  1:3797682/98‑1 (MQ=1)
       gggaTCGGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTCTTCCAGTTCCGGAAAACGCCGCTg                                                                                    <  2:3493986/98‑1 (MQ=0)
       gggaTCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGATTCAGGAAATCAACATGTCGGGTTTCCAGTTCCGGGAAACGCCGCtc                                                                                    >  1:3368044/4‑100 (MQ=1)
                 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)
                                                                                                   |                                                                                          
CAACCCCGGTATCAGTTCATCCAGCGCGGCTGCTTTGTTCATGGCTTTGATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCC  >  NC_000913/1633876‑1634065

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