Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_pgi	3,183,350	0	A	.	23.5%	41.0 / 10.4	17	intergenic (+27/+463)	zupT/ribB	zinc transporter/3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (6/7);  new base . (3/1);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCAAATGCA  >  NC_000913_3_pae_pgi/3183202‑3183493                                                                                                                                                      |                                                                                                                                                atgaATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTCGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACAACAGAAGTTGTCCA‑GCGCac                                                                                                                                                >  1:55620/1‑149 (MQ=255)    aaTTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGCGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGa                                                                                                                                             <  1:157024/150‑1 (MQ=255)             cGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTg                                                                                                                                    <  1:130434/150‑1 (MQ=255)                     aaaGAGATTGACCCTAATAATCACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACTTCt                                                                                                                                     <  4:41138/141‑1 (MQ=255)                     aaaGAGATTGACCCTAATAATCACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCATGCGC‑TGGACAACTTCt                                                                                                                                     >  3:41138/1‑141 (MQ=255)                         agaTTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGCAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGa                                                                                                                        >  1:77423/1‑150 (MQ=255)                                   taataataaCCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTAtt                                                                                                              >  1:67602/1‑150 (MQ=255)                                           aCCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAACACAACAACAGAAGTTATCCA‑TCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCttt                                                                                                      >  3:117127/1‑150 (MQ=255)                                                       gCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGtctct                                                                                          >  2:22777/1‑150 (MQ=255)                                                                       gTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTa                                                                          <  4:105508/150‑1 (MQ=255)                                                                                               gctgctACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCa                                                  >  2:55250/1‑150 (MQ=255)                                                                                                                                  aTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCg                                                                                   >  4:70440/1‑82 (MQ=255)                                                                                                                                  aTCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCg                                                                                   <  3:70440/82‑1 (MQ=255)                                                                                                                                         gTTGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCa        >  1:89594/1‑150 (MQ=255)                                                                                                                                          ttGTCCATGCGC‑TGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCaa       >  3:3837/1‑150 (MQ=255)                                                                                                                                           tGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCaaa      <  1:41670/150‑1 (MQ=255)                                                                                                                                               cA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGGGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCAAATGCa  <  1:149991/150‑1 (MQ=255)                                                                                                                                               cA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCAAATGCa  <  2:67602/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                ATGAATTAATGCCGCTCGCCAAAGAGATTGACCCTAATAATAACCCCAGCTATGGCGTACTATGTGGAATGTCAGTGATGGGATTCAGTTTAGTGCTGCTACAAACGGCGGGAATTGGTTAAAAACAACATCAGAAGTTGTCCA‑GCGCATGGACAACTTCTGCCAGCTTCACGAAAGCCTTATTCTACCTTTTCGGCTGTCTCTTCTCTCGTACTGTTTAAGTATTTGTTCCGCGCTTCAAGCAGTGCGTCTCTGTCCTGGCGAAACGCCCGACAATAACGTTTCAAATGCA  >  NC_000913_3_pae_pgi/3183202‑3183493

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