Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,157,289	0	A	G	43.8%	56.3 / 19.1	49	G171G (GGA→GGG)	ycfH	putative DNase
Reads supporting (aligned to +/- strand):  ref base A (6/21);  new base G (21/0);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 1.66e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGAAAAGAGAATCAAC  >  NC_000913/1157144‑1157437                                                                                                                                                   |                                                                                                                                                      tgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAAT‑GGGTCTCCGGttttt                                                                                                                                                  >  1:396979/1‑150 (MQ=255)          cGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATATT‑TGGTCTCTGAGTTTATATTTCCtt                                                                                                                                         >  1:684714/1‑151 (MQ=255)          cGGTTATAGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTCAATAATT‑TGATCTCTGCttttta                                                                                                                                                 >  2:744391/1‑141 (MQ=255)           ggTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTAAT‑TGGTCTCGGAttttttatttccctt                                                                                                                                        >  1:733691/1‑140 (MQ=255)            gTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCtttt                                                                                                                                       <  1:896220/151‑1 (MQ=255)                 cGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGc                                                                                                                                  <  1:416486/151‑1 (MQ=255)                 cGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACT‑GGGTCTCGGGGTTTAAATCTCCTTTTTCGGc                                                                                                                                  >  2:399978/1‑151 (MQ=255)                    ccATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCAtt                                                                                                                               <  2:192016/151‑1 (MQ=255)                      aTACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATtgt                                                                                                                             <  2:557741/151‑1 (MQ=255)                       tACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACT‑GGGTCTCGGGTTTTAAATCTCCTTTTTCCGCAtt                                                                                                                               >  2:216755/1‑148 (MQ=255)                        aCCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATtgtg                                                                                                                            <  1:827247/150‑1 (MQ=255)                                   ccGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAAt                                                                                                                <  2:287969/151‑1 (MQ=255)                                     gTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACT‑GGGTCTCGGGTTTTAAATCTCCTTTTTCGGCATTGTGACCTTCCGTaa                                                                                                                 >  2:194462/1‑148 (MQ=255)                                              cGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACt                                                                                                     <  2:658564/151‑1 (MQ=255)                                                  ggCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAATTAAT‑GGGTCTCGGGTTTTAAATCTTCTTTTTCCGCCTTGTGGCCTTCCCTAATGTGGGGGAAATgcg                                                                                                  >  1:890179/1‑150 (MQ=255)                                                      aTTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACT‑GGTTTTCGGGTTTTACATCTCCCTTTCCCGCCTTGTGGCCTTCCGTCATGCGGAGCAACAGCGCGCtg                                                                                             >  2:1068446/1‑151 (MQ=255)                                                              cgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGt                                                                                    <  2:983488/151‑1 (MQ=255)                                                               gaagaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAAT‑GGGTTTCGGGTTTTAAATCTCCTTTTTCCGCATTGTTGACTTCCCTAATGCGGGGGAACTGcgcggt                                                                                              >  1:998286/1‑139 (MQ=255)                                                                aagaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACT‑GGGACTCGGGTTTTAAATTTCCTTTTTCCGCATTTTTACCTTCCCTAAAGCGGGGGAACTGCGCGATGCTGcg‑cggt                                                                                   >  2:399440/1‑148 (MQ=255)                                                                 agaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTATT‑GGGTCTCGGGTTTTACATCTCCTTTTTCCGCCTTTTGGCCCTCCGTAATGCGGGGCAACAGCGCGAGGCTGCG‑CGGTAt                                                                                 >  2:465668/1‑151 (MQ=255)                                                                  gaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATTTCCGGTTTTTCAACTCCTTTTTCCGGATTGTGGACCTCCGTAAAGAGGGGCAAACGCGCGGTGCTGgg‑cgcttt                                                                                 >  2:1106352/1‑146 (MQ=255)                                                                   aaaaaGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTTCT‑GGGACTCCGGTTTTACATCTCCTTTTCCCGCATTGTGACCTTCCCTAATGCGGAGAAACTGCGCGATGCTGTG‑CGTTAt                                                                                 >  1:681929/1‑149 (MQ=255)                                                                     aaaGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACT‑GGATCTCGGGTTTTAAATCTCCTTTTTCGGCCTTGTGACCTTCCGTAATGCGGAGCAAAAGCGCGGTGCTGCG‑CGGTATTTcc                                                                             >  2:451667/1‑151 (MQ=255)                                                                          gACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAAATAAT‑AGATCTCCGGTTTTATATCTTCTTTTTCGCCATTGTGACATTCCCTAACTCGGAGCGACTGCGCGCTGCTCCG‑CGTCTTGTcaacc                                                                          >  1:687822/1‑145 (MQ=255)                                                                             ggATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGCAAATACT‑AGGTCTCGGAGTTTAAATTTCCTTTTCCCGCATTGTGAGCTTCCCTCAAGAGGCGCAACCGCGCGCTGCTGCG‑CGCTTTGTGCCCCCggg                                                                      >  2:288654/1‑149 (MQ=255)                                                                                  gCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAAT‑TTTGCTTTGATCTTAAATCTCCCTTTTTCGCCTGGGGTACCTCCCTAATTCGGTGGAACTGCGCGCTGCTTCGCCGTATTGTcccccgcgttgggt                                                                 >  1:145896/1‑142 (MQ=255)                                                                                   cGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTGAAATACC‑GGGACTCCGGTTTTACAACTCCTTTTTCACTAATGGGACCTTCCCGAATGCCGGGCCCCTGGGCGCTGCTGGC‑CGGTTTGTCCCCCtcgtacggt                                                                 >  1:179414/1‑141 (MQ=255)                                                                                   cGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTa                                                               <  1:940827/151‑1 (MQ=255)                                                                                    ggTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTAc                                                              <  2:647063/151‑1 (MQ=255)                                                                                       ggCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAAATACT‑GGATATCGGGTTTTAAATTTTCTTTTCCCGCGTTGTGAACTTCCGGCAAGTGGGGCAACCGCGCCGTGCGGCG‑CGGTTTGTCCCCCCCGATCgggttc                                                              >  2:665707/1‑144 (MQ=255)                                                                                        gCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACtggt                                                          <  2:245406/151‑1 (MQ=255)                                                                                           tactacACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATATT‑TGATCTCGGGTTTTTAATTTCCTTTTTCGGCCATTTGACCTTCCTTCATTCGGGGCAACAGCGCGATGTTGCG‑CGGTATTTCCCCCCCGGTCGGTTTCtggtgg                                                        >  2:859806/1‑150 (MQ=255)                                                                                            actacACTGTTTTACAGAGGACAGAGAAACGGCGGGAAGATTACTCGGAT‑TGGGGGTTTTCAACCCCCTTTTCGGCTATTTGGACTTCCGTGAAGCGGGGCAAAAGCGCGGGGCAGAG‑CGGTATTTGCCCCTCGacggggttttggtgg                                                        >  2:78392/1‑135 (MQ=255)                                                                                             ctacACTGTTTTACAGAGGACAGAGAAACGGCGGGCAAATTAAT‑TGATCTCCGGTTTTTTATATTCTTTTTCCGCATTTTGGTCTTCCTGAATGCGGAGGCAATGCGGGGTGGTCGG‑CGTTATTTTCCCCCGggttggttcatggggg                                                        >  1:760982/1‑140 (MQ=255)                                                                                              tacACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAAATAAT‑TGATCTCGTGTTTTAAATCTCCCTTTCCCGCCTTGTGGTCCTTCGTAATACTGGGGAACTGCGCGGTGGTGCG‑CGGTTTTTCCCCCCGGGTCCGGTACTGGTGGgaa                                                     >  1:228055/1‑150 (MQ=255)                                                                                               acacTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACt                                                   <  2:823529/151‑1 (MQ=255)                                                                                                cacTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTg                                                  <  1:1104951/151‑1 (MQ=255)                                                                                                    gTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGAc                                                <  2:747628/149‑1 (MQ=255)                                                                                                      tttACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCa                                             <  2:110102/150‑1 (MQ=255)                                                                                                               gACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTg                                   <  2:1018874/151‑1 (MQ=255)                                                                                                                  agagAAACGGCGGGGAAATTTAT‑TGATTTCTGGTTTTAAATCTCCTTTTTCGGCCGTGTGACCTTCCCTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTTTGTTCCCCCGGGTCGGGTACTGGTGGGGAATGACTCaccctacccttccc                                >  1:153232/1‑145 (MQ=255)                                                                                                                    agaAACGGCGGGGAAAATACT‑GGATCTCGGGTTTTAAATCTCCTTTTCCGGCATTGTGACCTTCCGTAATTCGGGGCAACTGCGCGGTGCTGCG‑CGTTTTGTTCCCCTGGGTCGGTTACTGGTGGGAACTGACTCACCCTACCTTTCGcc                               >  1:261794/1‑150 (MQ=255)                                                                                                                         cggcggGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTAcc                         <  1:1091705/151‑1 (MQ=255)                                                                                                                                 aaaTTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAg                 <  2:127103/151‑1 (MQ=255)                                                                                                                                   aTTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGa               <  2:397043/151‑1 (MQ=255)                                                                                                                                      aCT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGaaa             >  2:148575/1‑150 (MQ=255)                                                                                                                                               tcGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGAAAAGAGAATc     >  2:866741/1‑150 (MQ=255)                                                                                                                                                cGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGAAAAgaga        <  2:1107512/146‑1 (MQ=255)                                                                                                                                                 ggATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGAAAAGAGAATCAAc  <  2:173666/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      TGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACT‑GGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGAAAAGAGAATCAAC  >  NC_000913/1157144‑1157437

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