Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	336,865	0	G	.	63.6%	‑0.8 / 13.8	22	coding (88/1548 nt)	yahF	putative acyl‑CoA synthetase YahF
Reads supporting (aligned to +/- strand):  ref base G (2/6);  new base . (14/0);  total (16/6)
Fisher's exact test for biased strand distribution p-value = 3.75e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AAAATACAAAGAGGAATTTTATCTCGGTCTGCAACGCGGCAAAAATAATACCACATTATTAAGTGCCAT‑AACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCG‑TTTTATTCACAACAT‑TATTTATGACATCCCTGGG‑AACGCAACTCAGGCAATAGAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTT‑AAAATAG‑TCA‑TTAAACCGAATACCTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AAAC‑TC‑GACGGCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCC‑GGAGCTGGAGCAGGCG‑AAAAA‑CGGCG‑ACCTGATG‑ATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGC  >  NC_000913/336515‑337082                                                                                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                 aaaaTACAAAGAGGAA‑‑‑TATCTCGGTCTGC‑ACGCGGCAAAAATAATACCACATATTAAAGTGCCATAAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTTATTCACAACATATTTTTATGACATCCCTGGG‑AACGCAACTCAGGCAATAGAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCC‑AACCTACGGAGGGAATTATGTCAGTT‑AAAATAG‑TCA‑TTAAACCGAATACCTATTTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AAAC‑TC‑GACGGCGTCGAGCAGGCATTTTGTGGCGATGGCGACCGaa                                                                                                                                                                                              <  1:35914/386‑1 (MQ=255)                       ctcGGTCTGCAACGCGGCAAAAATAATACCACATTATTAAGTGCCAT‑AACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCG‑TTTTATTCACAACA‑‑TATTTATGACATCCCTGGG‑AACGCAACTCAGGCAATAGAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTT‑AAAATAG‑TCA‑TTAAACCGAATA‑CTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AAAC‑TC‑GACGGCGTCGa                                                                                                                                                                                                                           <  1:35915/334‑1 (MQ=255)                                                                                       ttACAACAACGCTGCCGGGAAAATATTCATCG‑TTTTATTCACAACA‑‑TATTTATGACATCCCTGGGAAACGCAACTCAGGCAATAGAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTT‑AAAATAG‑TCA‑TTAAACCGAATACCTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AAAC‑TC‑GACGGCGTCGAGCAGGCA‑tt                                                                                                                                                                                                                 <  1:754170/282‑1 (MQ=255)                                                                                                                                               gACATCCCTGGG‑AACGCAACTCAGGCAATAGAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTT‑AAAATAG‑TCA‑TTAAACCGAATACCTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GTAGACAAT‑AAAC‑TC‑GACGGCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAgg                                                                                                                                                                <  1:35930/277‑1 (MQ=255)                                                                                                                                                                     cAGGCAATAGAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTT‑AAAATAG‑TCA‑TTAAACCGAATACCTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AAAC‑TC‑GAC‑GCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCC‑GGAGCTGGAGCAgg                                                                                                                                     >  1:754183/1‑280 (MQ=255)                                                                                                                                                                                                  cTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTT‑AAAATAG‑‑CA‑TTAAACCGAATACCTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGCAGT‑AGCAAT‑AAACGTC‑GACGGCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCC‑GGAGCTGGAGCAGGCG‑aaa                                                                                                                               <  1:754184/258‑1 (MQ=255)                                                                                                                                                                                                                             tggcatggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTT‑AAAATAG‑TCA‑TTAAACCGAATACCTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AAAC‑TC‑GAC‑GCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCC‑GGAGCTGGAGCAGGCG‑AAAAA‑CGGCG‑ACCTGATG‑ATTGTCATCAATGGt                                                                                              >  1:754188/1‑259 (MQ=255)                                                                                                                                                                                                                             tggcatggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTT‑AAAATAG‑TCA‑TTAAACCGAATACCTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AAAC‑TC‑GAC‑GCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCC‑GGAGCTGGAGCAGGCG‑AAAAA‑CGGCG‑ACCTGATG‑ATTGTCATCAATGGTAAATCGGGTg                                                                                    >  1:754187/1‑269 (MQ=255)                                                                                                                                                                                                                                          ggTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTT‑AAAATAG‑TCA‑TTAAACCGAATACCTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AAAC‑TC‑GAC‑GCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCC‑GGAGCTGGAGCAGGCG‑AAAAA‑CGGCG‑ACCTGATG‑ATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGAtt                                                       >  1:754189/1‑285 (MQ=255)                                                                                                                                                                                                                                                               aaCGGAGGGAATTATGTCAGTTAAAAATAGCTCATTTAAACCGAATACCTA‑TTTTGATTCTGTCTCGCTGATGCTCTATCTCCACGC‑GT‑GCAAATAAAAC‑TCGGACGGCGTCGAGCAGGCA‑‑TTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCC‑GGAGCTGGAGCAGGCG‑AAAAA‑CGGCG‑ACCTGATG‑ATTGtcatc                                                                                                    <  1:35935/224‑1 (MQ=255)                                                                                                                                                                                                                                                                aCGGAGGGAATTATGTCAGTTAAAAATAG‑TCA‑TTAAACCGAATACCTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AAAC‑TC‑GACGGCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCC‑gg                                                                                                                                                 >  1:35936/1‑179 (MQ=255)                                                                                                                                                                                                                                                                aCGGAGGGAATTATGTCAGTT‑AAAATAG‑TCA‑TTAAACCGAATACCTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AAAC‑TC‑GAC‑GCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGGAGCTGGAGCAGGCG‑AAAAA‑CGGCG‑ACCTGATG‑ATTGTCATCAATGGTAAATCg                                                                                        >  1:754190/1‑231 (MQ=255)                                                                                                                                                                                                                                                                aCGGAGGGAATTATGTCAGTT‑AAAATAG‑TCA‑TTAAACCGAATACCTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AAAC‑TC‑GAC‑GCGTCGAGCAGGCA‑TTTGTGGCGATGGCG‑CCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGGAGCTGGAGCAGGCGAAAAAA‑CGGCG‑ACCTGATG‑ATTGTCATCAATGGTAAATCg                                                                                        >  1:754191/1‑231 (MQ=255)                                                                                                                                                                                                                                                                                          tAG‑TCA‑TTAAACCGAATACCTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AAAC‑TC‑GAC‑GCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCC‑GGAGCTGGAGCAgg                                                                                                                                     >  1:754192/1‑164 (MQ=255)                                                                                                                                                                                                                                                                                               cattaaaaCCGAATACCTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAATAAAAC‑TC‑GAC‑GCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCC‑GGAGCTGGAGCAGGCGAAAAAA‑ACGCGCACTCGATG‑ATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGtt                                                                    >  1:35937/6‑225 (MQ=255)                                                                                                                                                                                                                                                                                                a‑TTAAACCGAATACCTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AAAC‑TC‑GACGGCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCC‑GGAGCTGGAGCAGGCGAAAAAA‑CGGCG‑ACCTGATGAATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACg   >  1:754193/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                                     aaCCGAATACCTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AAAC‑TC‑GAC‑GCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCC‑GGAGCTGGAGCAGGCG‑AAAAA‑CGGCG‑ACCTGATG‑ATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTAc    >  1:754194/1‑280 (MQ=255)                                                                                                                                                                                                                                                                                                         gAATACCTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AAAC‑TC‑GAC‑GCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCC‑GGAGCTGGAGCAGGCG‑AAAAA‑CGGCg                                                                                                                       >  1:754196/1‑163 (MQ=255)                                                                                                                                                                                                                                                                                                              ccTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AAAC‑TC‑GAC‑GCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCC‑GGAGCTGGAGCAGGCG‑AAAAG‑CGGCG‑‑‑C‑GATG‑ATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGc  >  1:35938/1‑270 (MQ=255)                                                                                                                                                                                                                                                                                                                     ttGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AACT‑AC‑GAC‑GCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAATAAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCC‑GGAGCTGGAGCAGGCG‑AAAAAGCGGCG‑‑‑CTGATG‑ATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGtt                                                                    >  1:35941/1‑201 (MQ=255)                                                                                                                                                                                                                                                                                                                                       tG‑TCTATCTCCACGC‑GT‑GCAAATAAAAC‑TC‑GAC‑GCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAATAAAAGGCGTGc                                                                                                                                                                             >  1:35942/1‑85 (MQ=255)                                                                                                                                                                                                                                                                                                                                                tcCACGC‑GT‑GCAAAT‑AAAC‑TC‑GAC‑GCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTgaagaa                                                                                                                                                                      >  1:35943/1‑82 (MQ=255)                                                                                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                 AAAATACAAAGAGGAATTTTATCTCGGTCTGCAACGCGGCAAAAATAATACCACATTATTAAGTGCCAT‑AACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCG‑TTTTATTCACAACAT‑TATTTATGACATCCCTGGG‑AACGCAACTCAGGCAATAGAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTT‑AAAATAG‑TCA‑TTAAACCGAATACCTA‑TTTTGATTCTGTCTCGCTGATG‑TCTATCTCCACGC‑GT‑GCAAAT‑AAAC‑TC‑GACGGCGTCGAGCAGGCA‑TTTGTGGCGATGGCGACCGAAATGAAT‑AAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCC‑GGAGCTGGAGCAGGCG‑AAAAA‑CGGCG‑ACCTGATG‑ATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGC  >  NC_000913/336515‑337082

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