Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,812,792	0	C	A	45.5%	36.8 / 15.3	33	T21N (ACC→AAC)	dfp	fused 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase, FMN‑binding
Reads supporting (aligned to +/- strand):  ref base C (7/11);  new base A (15/0);  total (22/11)
Fisher's exact test for biased strand distribution p-value = 1.71e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGCGGT‑TTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGCCGAA  >  NC_000913/3812644‑3812940                                                                                                                                                      |                                                                                                                                                      cAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAAacc                                                                                                                                                     >  2:203662/1‑151 (MQ=255)             cAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAACCCCTGAACTg                                                                                                                                             <  2:144974/147‑1 (MQ=255)                      gCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAAACCCTGAACTGGGGCGGCGttt                                                                                                                                  >  1:389386/1‑149 (MQ=255)                        aCGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAAACCCTGAACTGGTGGGTCGttt                                                                                                                                  >  1:170790/1‑147 (MQ=255)                                   aGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTTTTGGCGCCTT‑TAAAAACACCCAACTGGGGCGTCGGTTGCGCGGTCGCgggg                                                                                                                    >  2:561330/1‑150 (MQ=255)                                       ttGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAAACCCCTAAATGGGGCGTCGGTTTGGCGGTAGCgggg                                                                                                                    >  2:423425/1‑146 (MQ=255)                                        tgtgATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTTGTTGCTA‑TAAAAACCCCGAACTAGTGCGTCGTTTGTGCGATCGCGGGGcc                                                                                                                  >  1:272266/1‑147 (MQ=255)                                               aaaTCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGGCGACGTCcgcg                                                                                                        >  1:189307/1‑150 (MQ=255)                                                      cAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTTCTGCCTA‑TAAAAACCCTGAACTGGGGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGGAGcc                                                                                                   >  1:214700/1‑148 (MQ=255)                                                           tCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGT‑AAAAAATCGTTCTCGGCGTTAGCGGCGGGATTGGTGCCTA‑TAAAAACCCTGCAATGGGGCGGCGTTTGCGCGCTCGCGGGGCCGAAGGCCGCGGAGCCCTGACCg                                                                                            >  2:581342/1‑149 (MQ=255)                                                               gtgtCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCg                                                                                       <  2:111484/151‑1 (MQ=255)                                                                    aCCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTTCTGCCTA‑TAAAAACCCTGAACTGGGGCGTCGTTTTCGCGATAGCGGGGGCGACGGCCGCGTAGGCCTGACCCAAGAGGcca                                                                                   >  1:311378/1‑148 (MQ=255)                                                                      ccTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAAACCCTGAACTGGTGCGGCGTTTGCGCGATCGCGGGGCCGACGTCCGCGGAGCCATGACCGAAGCGGCaaaa                                                                                 >  1:543045/1‑150 (MQ=255)                                                                              aaaaaTCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAAACCCTGAACTGGTGCGGCGTTTGCGCGATCGCGGGGCCGACGTCCGCGGAGCCCTGACCGAAGCGGCAAAAGCCttttt                                                                         >  1:205284/1‑148 (MQ=255)                                                                                    catcatGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCa                                                                  <  1:129906/151‑1 (MQ=255)                                                                                            gCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCtt                                                          <  1:687243/151‑1 (MQ=255)                                                                                             ccTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCCGCCTA‑TAAAAACCCTGAACCGGGGGGTCGTTTGCGCGGTCGCGGGGCCGCCGTCCCCGTATCCCAGGCCGCAGCGGCCAAAACCTTTTTCACCCCCCCTAGCt                                                           >  2:86195/1‑149 (MQ=255)                                                                                               tGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAAACCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGGAGCCCTGACCGAAGCGGGAAAAGCCTTTTTCACCCCACCTAGCtt                                                          >  1:276456/1‑148 (MQ=255)                                                                                                ggCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGc                                                        <  2:99731/149‑1 (MQ=255)                                                                                                  ccGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGc                                                    <  2:345387/151‑1 (MQ=255)                                                                                                       aaaaaaaTCGTTCTCGGCGTTAGCGGCGGTATTTCTGCCCA‑TAAAAACCCCTAACTGCGGCGGCGGTTGCGCGAACGCGGGGCCGGCGTCCGCGCAGCCCTGACCCAAGCGGCCAAAGACCTTTTCACCCCCCCTAACTTGCAGGCGGg‑tt                                              >  2:128688/1‑151 (MQ=255)                                                                                                       aaaaaaaTCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAAACCCCGAACTGGTGCGGCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAACCTTTTTCACCCCACTTAGCTTGCAGGCGGg‑tt                                              >  1:587644/1‑151 (MQ=255)                                                                                                           aaaTCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGCGGT‑TTCTgg                                          <  2:264502/151‑1 (MQ=255)                                                                                                              tCGTTCTCGGCGTTAGCGGCGGCATTGTCTCCTATTAAAACCCCCG‑ACTGGGGCGTCGTTTTCTGGATCGCGGGGCCGACGTCCGCGGAGGCCTGACCCAAGCGGGCAAAGACCTTTTCACCCCACCTCGCTTGCCGGCGGTGTTTTGGtt                                        >  2:613160/1‑151 (MQ=255)                                                                                                                   ctcGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAAACCCTGAAATGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGGAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCCCCTTGCTTGCAGGCGGg‑tttttgtgtttcc                                   >  1:367490/1‑139 (MQ=255)                                                                                                                       gCGTTAGCGGCGGTATTGCTGCCTA‑TAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGCGGT‑TTCTGGTTATCCCGTTTc                              <  2:64728/151‑1 (MQ=255)                                                                                                                                       tgctgcCTA‑TAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGCGGT‑TTCTGGTTATCCCGTTTCCGACAGTCTGCTGGa               >  1:42520/1‑150 (MQ=255)                                                                                                                                         ctgcCTA‑TAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGCGGT‑TTCTGGTTTTCCCGTTTCCGACAGTCTGCTGGAccc            >  2:365425/1‑151 (MQ=255)                                                                                                                                           gcCTA‑TAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGCGGT‑TTCTGGTTATCCCGTTTCCGACAGTCTGCTGGAccc            <  2:84211/149‑1 (MQ=255)                                                                                                                                           gcCTA‑TAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGCGGT‑TTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCg           <  2:497506/150‑1 (MQ=255)                                                                                                                                            ccTA‑TAAAACCCCTGAACTGGCGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGCGGT‑TTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCgg          >  2:420992/1‑150 (MQ=255)                                                                                                                                                 taAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGCGGT‑TTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGcc     >  2:712274/1‑151 (MQ=255)                                                                                                                                                    aaCCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGCGGT‑TTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCGGCCGa   >  2:235278/1‑150 (MQ=255)                                                                                                                                                    aaCCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGCGGT‑TTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGCCGaa  <  1:469850/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTA‑TAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGCGGT‑TTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGCCGAA  >  NC_000913/3812644‑3812940

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