Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	670,328	0	G	A	60.0%	2.4 / 17.0	20	P82L (CCC→CTC)	nadD	nicotinate‑nucleotide adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (5/3);  new base A (6/6);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGGTGAAGATCTTCCGGGTTATGTGTCAAATGATCTTCCAGCCATTGCTGGTATTGCGGTTGCGCCATTTCAAGTGGGTAACCTGGACGCCGACAGACGATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAGCATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTT  >  NC_000913/670085‑670513                                                                                                                                                                                                                                                    |                                                                                                                                                                                          aGGTGAAGATCTTCCGGGTTATGTGTCAAATGATCTTCCAGCCATTGCTGGTATTGCGGTTGCGCCATTTCAAGTGGGTAACCTGGACGCCGACAGACGATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCAtt                                                                                                                                                                                    >  1:188764/1‑251 (MQ=255)                         gtCAAATGATCTTCCAGCCATTGCTGGTTTTGCGGTTGCGCCATTTCAAGTGGGTAACCTGGACGCCGACAGACGAACAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGAGGGCATTGCGCTTTAGCTCGCGTTCATCAag                                                                                                                                                            <  1:275207/251‑1 (MQ=255)                                gATCTTCCAGCCATTGCTGGTATTGCGGTTGCGCCATTTCAAGTGGGTAACCTGGACGCCGACAGACGATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTGTCAGTGTTTGCGCAGTGTAAGAGAGGGCATTg                                                                                                                                                                                   <  1:114023/221‑1 (MQ=255)                                                                                          ccGACAGACGATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAGCATGTGTTTACGCTGCACGCTGTTCGCtt                                                                                           >  1:2092/1‑251 (MQ=255)                                                                                                                                cttattcgtaCCAGGTCGGAAAGGCCAGCAGTGAATCATGACCAATAATAAACGACCCCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTAAGTGTTTGCGCAGTGTAAGAGAGGGCATTGCGCTTTAGCTCGCGTTCCTCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGCTCAAGCATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACCATGCGGAGGAACATTATTAGGgatg                                                     >  1:308698/3‑251 (MQ=255)                                                                                                                                         aCCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAGCATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGAc                                             <  1:27636/250‑1 (MQ=255)                                                                                                                                                 ggAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCCGTGTAAGAGAGGGCATTGCGCTTTAGCTCGCGTTCATCCAGAGTAAATAATGGCTTGTCGGCAAt                                                                                                                                    >  1:162108/1‑155 (MQ=255)                                                                                                                                                              gTGAATCCTGACCAATAATAAACGCCAGAGGCACGTCCGGTACTTGTTCTTGCCCCCCCTCTTTCAGTGATTGAGCAGTGTGAGAGGGGGCATTGCGCTTTCGCTCGCGATCATCAAGAGAAAATAATGACTACTCGGCCATCGCCAGTTCAAGCATGGGTTTACGCTGGACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGAGGATTGTGACCAGCGTCAGACCAATCAAAttt                       >  1:165107/1‑250 (MQ=255)                                                                                                                                                                naaTCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTACTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATACTGGCTTg                                                                                                                                            >  1:39521/2‑132 (MQ=255)                                                                                                                                                                    ccTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGAGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGt                                                                                                                                           >  1:164809/1‑129 (MQ=255)                                                                                                                                                                         ncAATAATAACCGCCAGCGGCACGTCAGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAACTAATGGCTTGTCGGACATCGCCAGTTCAAGCATGTGTTTACGCTGCACGCTGTTCGc                                                                                             >  1:84806/2‑170 (MQ=255)                                                                                                                                                                               aaaaaCGCCAGCGGCACGTCCGGTCGTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGAGTGCAGTGCGCTTTAGCTCGCGTTCATCAAGAGAAAATAATGGCTTGTCGGCAATCGCCAGTTCAAGCATGGGTTTACGCTGCACGCTGt                                                                                                 <  1:299501/158‑1 (MQ=255)                                                                                                                                                                                 nnaCNNNNNCGGNACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAGCATGTGTTTAc                                                                                                             <  2:197715/144‑1 (MQ=255)                                                                                                                                                                                    cGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCAACTCCTTCAGTGTTTGCGCAGTGTAAGAGAGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAGCATGTGTTTATTTTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTTTTAGGGATGATTGTGAACCGCGTCCGACCAATCAAATTCGCCAGCGTTTCCACGGGtttt  >  1:110354/1‑250 (MQ=255)                                                                                                                                                                                     gCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCAGCAAGAGTAAATAATGGCTTGTCg                                                                                                                                         <  1:260563/114‑1 (MQ=255)                                                                                                                                                                                            gCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGAGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAGCATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAAAATTATTAGGGATGATTGTGACCCGCGTCAGAcc                                 >  1:264095/1‑211 (MQ=255)                                                                                                                                                                                                gTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGAGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAGCATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGAcn                                 <  1:193971/207‑2 (MQ=255)                                                                                                                                                                                                gTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGAGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAGCATGTGTTTACGCNGCANNNNNttnn                                                                                              >  2:193971/1‑144 (MQ=255)                                                                                                                                                                                                                      cACTCTTTCAGTGTTTGCGCAGTGTAAGAGAGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAGCATGTGTTTACGCTGCACGCTGTTCGCTTCCGGc                                                                                      <  1:288338/132‑1 (MQ=255)                                                                                                                                                                                                                                       cgcAGTGTAAGAGAGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAGCATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTc                                      <  1:303673/163‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                          AGGTGAAGATCTTCCGGGTTATGTGTCAAATGATCTTCCAGCCATTGCTGGTATTGCGGTTGCGCCATTTCAAGTGGGTAACCTGGACGCCGACAGACGATCAAATGTGCATTGTCGAGTATCGTTTCGTATTCGTACCAGGTCGGAAAGGTCAGCAGTGAATCCTGACCAATAATAAACGCCAGCGGCACGTCCGGTCCTTGTTCCTGCCGCCACTCTTTCAGTGTTTGCGCAGTGTAAGAGGGGGCATTGCGCTTTAGCTCGCGTTCATCAAGAGTAAATAATGGCTTGTCGGCAATCGCCAGTTCAAGCATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTT  >  NC_000913/670085‑670513

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