Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	4,808,173	0	G	.	45.5%	15.5 / 15.8	11	intergenic (‑282/‑273)	PP_4234/dsbD‑II	Aldehyde oxidase and xanthine dehydrogenase/thiol/disulfide interchange protein 2
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base . (3/2);  total (6/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGACATGCTTGATCCTCTGGGCAGTTGGGGCCTGCCCTTTCTTTTGATTCTTGGCCTGGCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTC  >  AE015451/4807886‑4808390                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                          ggACATGCTTGATCCTCTGGGCAGTTGGGGCCTGCCCTTTCTTTTGATTCTTGGCCTGGCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCcgcgcggcgct                                                                                                                                                                                                              <  1:828162/301‑8 (MQ=255)                     gCAGTTGGGGCCTGCCCTTTCTTTTGATTCTTGGCCTGGCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCcgcgcggcgcttcgcgggtaaccccg                                                                                                                                                                                               >  1:397310/1‑274 (MQ=255)                                    ccTTTCTTTTGATTCTTGGCCTGGCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCcgcgcggcgctcgatct                                                                                                                                                                                                        <  1:224382/272‑14 (MQ=255)                                    ccTTTCTTTTGATTCTTGGCCTGGCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCcgcgcggcgctcgatct                                                                                                                                                                                                        >  2:224377/1‑259 (MQ=255)                                                                 gaAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGcgccgcgcgggcggcgctcgatctcaca                                                                                                                                                                                                    <  2:641313/247‑23 (MQ=255)                                                                                 cGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGcgccgcgcgggcggcgctcgatc                                                                                                                                                                                                         >  1:308476/1‑209 (MQ=255)                                                                                                 aTTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGcgccgcgcgggcggcgctcgatc                                                                                                                                                                                                         <  2:308468/210‑18 (MQ=255)                                                                                                                  aaaGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGccgccccct                                                                                                                                                                                                                     <  1:11788/181‑4 (MQ=255)                                                                                                                  aaaGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGccgccccct                                                                                                                                                                                                                     >  2:11786/1‑178 (MQ=255)                                                                                                                    aGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTcc                                                                                                                                       <  1:227982/243‑1 (MQ=255)                                                                                                                    aGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTcc                                                                                                                                       >  2:227977/1‑243 (MQ=255)                                                                                                                                   gtgtCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCAcagctcgga                                                                                                                             >  2:689625/1‑230 (MQ=255)                                                                                                                                   gtgtCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCAc                                                                                                                                     <  1:689650/230‑1 (MQ=255)                                                                                                                                                                                                                                    aTTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTc                                                                                                                                   <  1:920811/135‑1 (MQ=21)                                                                                                                                                                                                                                    aTTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTc                                                                                                                                   >  2:920785/1‑135 (MQ=21)                                                                                                                                                                                                                                             gcCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTc  >  1:472587/1‑255 (MQ=255)                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                          GGACATGCTTGATCCTCTGGGCAGTTGGGGCCTGCCCTTTCTTTTGATTCTTGGCCTGGCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCGTTAACCTTC  >  AE015451/4807886‑4808390

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