Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	581,227	0	A	G	66.7%	10.2 / 10.5	18	I132V (ATC→GTC)	nohD	putative DNA‑packaging protein NohD
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base G (5/7);  total (9/9)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGCAGGCCAGCGAGACAGATCTCCAGCCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGA‑TCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGG  >  NC_000913/581040‑581425                                                                                                                                                                                            |                                                                                                                                                                                                        gctgcAGGCCAGCGAGACAGATCTCCAGCCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGA‑TCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCATGAAAc                                                                                                                                         >  1:122829/1‑251 (MQ=17)                    gATCTCCAGCCAGGGATTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGAGAAACTGCCAGTATTCTCGACGGGA‑TCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAACCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCAt                                                                                                                      >  1:175336/1‑251 (MQ=9)                                             aCGAACGCCATTGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGA‑TCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAg                                                                                                                          <  1:3542/222‑1 (MQ=14)                                                                                    aggaactaaaaaaTGCAAAAGATTCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGAGGATGGCCGGAGACGTAGCCAGCATTCTTGATGGAG‑TTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTc                                                       >  1:77666/9‑250 (MQ=255)                                                                                               aatgcaaaagattCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTTATGGAG‑TTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACAGCATTAAAGCCATGAACAAATCAGCTGCGCTGGATGAAATAATACCGTGGTTGCTGAGTGAATATATCGAACAGTCAGGTTaaggggc                                           <  1:120193/239‑6 (MQ=255)                                                                                                     aaagattCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAG‑TTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGcagc                                     <  1:41001/245‑4 (MQ=2)                                                                                                         attCCGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAG‑TTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTaaggggcttcagcacnc                                 <  1:136550/249‑16 (MQ=2)                                                                                                             cGCTGAGGTGGTGGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAG‑TTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcactctgcg                             <  1:72937/251‑8 (MQ=2)                                                                                                                 naaGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGA‑TCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGACACCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGaaa                                                           >  1:22387/2‑216 (MQ=17)                                                                                                                    ggggtggAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAG‑TTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcactctgcgcgcgcgg                      >  1:159252/3‑237 (MQ=255)                                                                                                                     tggtggAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAG‑TTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcactctgcgcgcgcggg                     >  1:13835/1‑236 (MQ=255)                                                                                                                       gtggAAACCGCATTCTGCACGTTCGTGCTGTCGCTGATGGCCGGAGAAGAAGCCAGCATTCTTGATGGAG‑TTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGCCATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcactctgcgcgcgcgggtc                   >  1:39174/1‑234 (MQ=2)                                                                                                                         ggAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCAGAGAAGTAGCCAACATTCTTGATGGAG‑TTCCTCTGTCGGTTCAGCGGAGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCCGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcactctgcgcgcgcgggtct                  >  1:98379/1‑232 (MQ=255)                                                                                                                              ccGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAG‑TTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcactctgcgcgcgcgggtctactttc            <  1:50260/251‑25 (MQ=255)                                                                                                                              ccGCAGTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAG‑TTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATAGTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTGGCGGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCGGCactctgcgctctcgggtcttctttc            <  1:135619/251‑25 (MQ=255)                                                                                                                               cGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTG‑CTGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAACACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAACTAATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcac                                   >  1:94972/1‑226 (MQ=2)                                                                                                                                  aTTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAG‑TTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATCTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCTGCGCTGGATGAAATAATACCGGGGTTGCGGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTgcagcactctgcgcgcgcggg                     <  1:121768/238‑16 (MQ=255)                                                                                                                                                                              gTATTCTCGACGGGA‑TCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAgg  <  1:50281/213‑1 (MQ=255)                                                                                                                                                                                            |                                                                                                                                                                                                        GCTGCAGGCCAGCGAGACAGATCTCCAGCCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGA‑TCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGG  >  NC_000913/581040‑581425

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