Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,478,240	1	.	C	23.5%	29.7 / 13.9	17	coding (258/591 nt)	slyD	FKBP‑type peptidyl‑prolyl cis‑trans isomerase SlyD
Reads supporting (aligned to +/- strand):  ref base . (10/3);  new base C (2/2);  total (12/5)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCACGACGTGATCGTCTTCAACCGCAGTGATTTCAACCGGTACCGGACCCTGGTCGGTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACG‑CCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGG  >  NC_000913/3478103‑3478379                                                                                                                                           |                                                                                                                                            aaTTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCACGACGTGATCGTCTTCAACCGCAGTGATTTCAACCGGTACCGGACCCTGGTCGGTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACG‑cc                                                                                                                                           <  1:435223/140‑1 (MQ=255)                             ggTTACCATCAACCACGACGTGATCGTCTTCAACCGCAGTGATTTCAACCGGTACCGGACCCTGGTCGGTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCCATAAATACGTCTTTAGGAACACg                                                                                                                   <  2:552467/137‑1 (MQ=255)                             ggTTACCATCAACCACGACGTGATCGTCTTCAACCGCAGTGATTTCAACCGGTACCGGACCCTGGTCGGTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCCATAAATACGTCTTTAGGAACACg                                                                                                                   <  2:542992/137‑1 (MQ=255)                             ggTTACCATCAACCACGACGTGATCGTCTTCAACCGCAGTGATTTCAACCGGTACCGGACCCTGGTCGGTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCCATAAATACGTCTTTAGGAACACg                                                                                                                   >  1:542992/1‑137 (MQ=255)                             ggTTACCATCAACCACGACGTGATCGTCTTCAACCGCAGTGATTTCAACCGGTACCGGACCCTGGTCGGTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCCATAAATACGTCTTTAGGAACACg                                                                                                                   >  1:552467/1‑137 (MQ=255)                              gTTACCATCAACCACGACGTGATCGTCTTCAACCGCAGTGATTTCAACCGGTACCGGACCCTGGTCGGTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACG‑CCCATAAATACGTCTTTAGGAACACGTTGCa                                                                                                              >  1:505345/1‑140 (MQ=255)                                             gacgTGATCGTCTTCAACCGCAGTGATTTCAACCGGTACCGGACCCTGGTCGGTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACG‑CCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTcgtcgt                                                                                               >  2:418333/1‑140 (MQ=255)                                                                  aGTGATTTCAACCGGTACCGGACCCTGGTCGGTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACG‑CCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTcg                                                                          <  2:356667/140‑1 (MQ=255)                                                                                gTACCGGACCCTGGTCGGTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACG‑CCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACa                                                            >  1:510844/1‑140 (MQ=255)                                                                                           tGGTCGGTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACG‑CCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTc                                                 >  1:494836/1‑140 (MQ=255)                                                                                                            gAAACGCATACCTACCTGCAGTTCATCAACG‑CCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCtt                                >  2:196350/1‑140 (MQ=255)                                                                                                                           ctGCAGTTCATCAACG‑CCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCa                 >  1:26275/1‑140 (MQ=255)                                                                                                                               aGTTCATCAACG‑CCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGcc             >  2:14776/1‑140 (MQ=255)                                                                                                                                gTTCATCAACG‑CCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCa            <  2:429606/140‑1 (MQ=255)                                                                                                                                 ttCATCAACG‑CCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCag           >  1:408135/1‑140 (MQ=255)                                                                                                                                  tcatcaACG‑CCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCATTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCATCGCCGTTTCCAGGCCaga          <  2:505345/140‑1 (MQ=255)                                                                                                                                    atcaACG‑CCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCaga          >  2:151061/1‑138 (MQ=255)                                                                                                                                    atcaACG‑CCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCaga          <  1:151061/138‑1 (MQ=255)                                                                                                                                         cG‑CCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAgg   >  2:232810/1‑140 (MQ=255)                                                                                                                                          g‑CCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAggg  <  2:510844/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                            AATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCACGACGTGATCGTCTTCAACCGCAGTGATTTCAACCGGTACCGGACCCTGGTCGGTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACG‑CCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGG  >  NC_000913/3478103‑3478379

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