Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,027,436	0	T	G	40.9%	12.7 / 18.7	22	intergenic (+208/+173)	fre/fadA	NAD(P)H‑flavin reductase/3‑ketoacyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base T (7/6);  new base G (9/0);  total (16/6)
Fisher's exact test for biased strand distribution p-value = 4.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGT  >  NC_000913/4027193‑4027682                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         gATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATTCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGGACGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGAGGCCTGATGCGACGCTGGCGCGTCt                                                                                                                                                                                                                                                   >  1:201438/1‑251 (MQ=255)         gTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGTGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATcaggc                                                                                                                                                                                                                                           >  1:106558/1‑247 (MQ=255)                     gTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTCAACTGTTAATCAGTGGCATTTAGCTCTATGACGTATCTTGCAAACCATGCCCGAGGCGACGCTGTCGCGTCTTATCGTGCCTACCAATAGTCCGAACCGGAGGCCGGATAAGGCGTTTACGCCGCCTCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATcaggccttcattggtgc                                                                                                                                                                                                                               >  1:319633/1‑235 (MQ=21)                                    agaTATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACCCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACattggtgccggatcggtagaccg                                                                                                                                                                                                                >  1:71816/1‑229 (MQ=18)                                                                    aagaaCGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAT‑TGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCTCATCCGGCaagagg                                                                                                                                                                               >  2:55047/1‑247 (MQ=11)                                                                    aagaaCGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAT‑TGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTgg                                                                                                                                                                               >  1:305515/1‑251 (MQ=11)                                                                    aagaaCGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAT‑TGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTgg                                                                                                                                                                               >  2:455701/1‑251 (MQ=11)                                                                     agaaCGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAT‑TGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTggg                                                                                                                                                                              >  1:176919/1‑250 (MQ=11)                                                                       aaCGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAT‑TGGTGCCGGATCGGTAGACCGGATAAGGCGATCACGCCGCATCCGGCAAGTGGtt                                                                                                                                                                             >  1:397517/1‑249 (MQ=11)                                                                                                                                                                                                         agatctatgacgtatctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGAc                                         >  1:58029/26‑251 (MQ=18)                                                                                                                                                                                                          gatctatgacgtatctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGAcc                                        >  1:475718/25‑251 (MQ=18)                                                                                                                                                                                                             ctatgacgtatctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGAccc                                       >  2:206459/22‑249 (MQ=18)                                                                                                                                                                                                              tatgacgtatctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGaacc                                       >  1:202749/21‑245 (MQ=18)                                                                                                                                                                                                                 gacgtttctggcaacccaggcccgATGCGACGCTGTCGCGACTTATCGTGCCTACAAATAGT‑CCG‑AACCGGAGGCCGGATAAGGCGGTTACGCCGCATCCGGCAATTGGTGCATGATGCCGGATGCGACGCTGGCGCGTCTTAGCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCAGCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGAcc                                  <  2:164906/227‑1 (MQ=21)                                                                                                                                                                                                                   cgtatctggcagaccctgCCCGATGCGACGCTGTCGCGTCTTAGCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCTGGATAAGGCGTTTACGCCGCTTCCGGCAATTGGTGCATGATGCCTTATGCGTCGCTGGCGCGTCTTATCAGGCGTCCATTGGTGCCGGATCGGTAGACCGGGTAAGGCGGTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTGAAACACCGTCGCAATACCCTGACCCAGACCGAt                               <  1:131843/235‑1 (MQ=255)                                                                                                                                                                                                                    gtatctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATa                              <  1:310256/237‑1 (MQ=255)                                                                                                                                                                                                                             aaaccatgCCCGATGCGGCGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACAtcgtc                     >  1:106010/6‑251 (MQ=255)                                                                                                                                                                                                                                  atgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGcc                  >  1:412508/1‑249 (MQ=255)                                                                                                                                                                                                                                    gcccgATGCGCCGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGAc              <  2:334302/247‑1 (MQ=255)                                                                                                                                                                                                                                         aTGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAAc         >  1:446678/1‑251 (MQ=255)                                                                                                                                                                                                                                           gCGACGCGGTCGCGTCTTATCGTGCCTACAAATAGC‑CCG‑ACCCGTAGGCCGGATAAGTCGTTTACGCCGCATCCGGGAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTg       <  2:202749/251‑1 (MQ=255)                                                                                                                                                                                                                                                gCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGACAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGt  <  1:344537/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         GATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CCG‑AACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGT  >  NC_000913/4027193‑4027682

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