Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,386,862	0	G	A	63.2%	12.8 / 34.7	38	intergenic (+31/‑72)	yfbL/yfbM	putative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/14);  new base A (24/0);  total (24/14)
Fisher's exact test for biased strand distribution p-value = 1.03e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-05

CTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAGCGACATGCA‑CA‑TTTTA‑TTTTCCCTCATCTG‑TTTTAAG‑ATGC‑G‑AAACATGAAATTC‑AACGGA‑GGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAA  >  NC_000913/2386716‑2387003                                                                                                                                                    |                                                                                                                                                        cTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCATACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCC‑AAAAAACgga                                                                                                                                                     >  2:568300/1‑149 (MQ=255)   ttttATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAGCGACAt                                                                                                                                                   <  1:466500/151‑1 (MQ=255)            aataaACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAGCGACATGCA‑Ca‑tttt                                                                                                                                        <  2:12535/151‑1 (MQ=255)                aaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAGCGACATGCA‑Ca‑tttta‑ttt                                                                                                                                   <  2:113393/151‑1 (MQ=255)                  cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCG‑AAAAAACGACAAGCA‑Ca‑ttttt‑ttttt                                                                                                                                 >  2:42833/1‑150 (MQ=255)                       ccACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGCTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCG‑AAAAAAAGACAAGGC‑Ca‑attta‑ttttttctc                                                                                                                             >  2:375559/1‑145 (MQ=255)                         ccTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACTACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAGCGACATGCA‑CA‑TTTTA‑TTTTCCCTCATc                                                                                                                          <  1:487615/150‑1 (MQ=255)                                cGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAG‑AAAAAAAGACAAGCA‑CA‑TTTTA‑TTTTTCCTCCTCTG‑ttttt                                                                                                                  >  2:119857/1‑150 (MQ=255)                                cGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAACCACATTCA‑CA‑ATTTT‑TTTTTCCCCACCtt‑ttttt                                                                                                                  >  2:508712/1‑150 (MQ=255)                                  gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAG‑AAAAAAAGAAAAGCA‑AA‑ATTTT‑TTTTTTCTCAttcctttttaa                                                                                                                 >  2:156959/1‑141 (MQ=255)                                   tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAGCGACATGCA‑CA‑TTTTA‑TTTTCCCTCATCTG‑TTTTAAG‑a                                                                                                              <  2:223787/151‑1 (MQ=255)                                   tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAGCGACATGCA‑CA‑TTTTA‑TTTTCCCTCATCTG‑TTTTAAG‑a                                                                                                              <  2:478043/151‑1 (MQ=255)                                       aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCG‑AAAAAAAGAAAAGCC‑AA‑ATTTT‑TTTTTTCTCTTCTGATTTTAAA‑ATgg                                                                                                           >  2:74459/1‑150 (MQ=255)                                              acagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAGCGACATGCA‑CA‑TTTTA‑TTTTCCCTCATCTG‑TTTTAAG‑ATGC‑G‑AAAc                                                                                                    <  1:47619/148‑1 (MQ=255)                                                              gAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAG‑AAAAAAAGAAAAGGA‑CA‑TTTTT‑TTTTTCCTCCTCTGTTTTTAAA‑ATG‑‑G‑AAAAAAAAAATTT‑AACGGG‑GGaa                                                                               >  2:130600/1‑151 (MQ=255)                                                                   tGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAGCGACATGCA‑CA‑TTTTA‑TTTTCCCTCATCTG‑TTTTAAG‑ATGC‑G‑AAACATGAAATTC‑AACGGA‑GGAAaaat‑a                                                                         <  1:297305/151‑1 (MQ=255)                                                                   tGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAGCGACATGCA‑CA‑TTTTA‑TTTTCCCTCATCTG‑TTTTAAG‑ATGC‑G‑AAACATGAAATTC‑AACGGA‑GGAAaaat‑a                                                                         <  2:81919/151‑1 (MQ=255)                                                                    ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCG‑AAAAAACCACATGGA‑CA‑TTTTT‑TTTTTCCTCCTCTG‑TTTTTTA‑ATGC‑G‑AAAAAATAAAATT‑AAAAGA‑GGGAaaat‑aa                                                                        >  1:452251/1‑151 (MQ=255)                                                                          ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAAAGGAAAGGCC‑CA‑‑TTTA‑TTTTTTCTCCTCTG‑TTTTTAA‑ATGCGG‑AAAAAAAAAATTA‑AACGGG‑GGAAAAAA‑AAATggg                                                                   >  2:605037/1‑151 (MQ=255)                                                                          ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAG‑AAAAAAAGACAAGGA‑CA‑TTTTT‑TTTTTCCTCATCTG‑TTTTTAG‑ATGG‑G‑AAACAAAAAAATT‑AACAGG‑GGAAAAAA‑AAAAgggg                                                                  >  2:374552/1‑150 (MQ=255)                                                                              gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCTG‑AAAAAACGACAACAA‑CA‑ATTTA‑TTTTTCCTCATCTG‑TTTTTAG‑AAGC‑G‑AAAAAAGAAAATT‑AACAGG‑GGGAAAAT‑TAATTGGGATGa                                                              >  2:443370/1‑151 (MQ=255)                                                                                ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCG‑AAAAAACGACAAGCA‑AA‑ATTTA‑TTTTTCCTCCTCTG‑TTTTTAG‑AAGG‑G‑AAAAAAGAAAATC‑AACGGA‑GGAAAAAA‑AAATTGGGAAGAtt                                                            >  2:118914/1‑151 (MQ=255)                                                                                     gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAT‑AAAAAAAGAAATGCA‑CA‑TTTTT‑TTTTTTCTCCTCTG‑TTTTTAA‑ATGC‑G‑GAAAATGAAAATT‑AAACGA‑GGGAAAAT‑AAATGGGGATGATTgggtt                                                       >  1:174597/1‑148 (MQ=255)                                                                                     gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAAAGACAAGGA‑CA‑TTTTT‑TTTTTCCTCTTTTG‑TTTTTAGAAGGT‑G‑AAAAATAAAATTT‑AAAGGG‑GGAAAAAA‑AAATTGG‑ATGAtttggta                                                       >  2:246896/1‑146 (MQ=255)                                                                                     gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAAAGACAAACA‑CA‑TTTTT‑TTTTTCCTCTTTTG‑TTTTTAT‑AAGC‑GAAAAAAATAAATTA‑AAAGGG‑GGAAAATA‑AAATggaatggtttggtt                                                        >  2:306852/1‑139 (MQ=255)                                                                                     gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCG‑AAAAAAAGACAAGGA‑AA‑TTTTT‑TTTTTCCTCATCTG‑TTTTTAG‑ATAC‑G‑AAAAATGAAATTT‑AAATGA‑GGGAAAAA‑AAAATGGAATGATTgggtt                                                       >  1:44061/1‑148 (MQ=255)                                                                                       gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCTG‑AAAAAAAGAAATAGA‑AA‑ATTTA‑TTTTCTCTTTTCAT‑TTTTAAA‑ATGCGG‑AAAAATAAAATTT‑AACGG‑‑GGAAAAAA‑AAATTGGGGAGATTGGGTAtt                                                     >  2:186374/1‑151 (MQ=255)                                                                                                ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCG‑AAAAAAAGAAATAAA‑CA‑TTTTA‑TTTTTCCTCATCTT‑TTTTAAG‑ATGC‑G‑AAAAAAAAAAATT‑AACGGAGGGAAAAAA‑AAAGGGGAAT‑ATTTGGTAttttttataaa                                            >  2:568349/1‑144 (MQ=255)                                                                                                ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAAAGACAAGCA‑AA‑TTTTT‑TTTTTCCTCATCTT‑TTTTTAG‑ATGC‑G‑AAACAAGAAATTC‑AAAGGA‑GGGAAAAT‑AAATTGGGATGATTGGCTATTTTTCAGaaa                                            >  2:414605/1‑151 (MQ=255)                                                                                                  gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAACGACATGCA‑CA‑TTTTT‑TTTTTCCTCATCTG‑TTTTTAA‑AAGC‑G‑AAACAAGAAAATC‑AACGGA‑GGAAAAAA‑AAATTGGAATGATTGGCTAtttttttgaaaat                                          >  2:459912/1‑142 (MQ=255)                                                                                                  gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAACGACAAGAACCA‑TTTTT‑TTTTCCC‑CACCTT‑TTTTTAA‑ATGC‑G‑GAACAAGAAAATCAAACGGG‑GGAAAAA‑‑AAATGGGGATGATTTGCTAttttttcgaaaat                                          >  2:543104/1‑142 (MQ=255)                                                                                                     ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAGCGACATGCA‑CA‑TTTTA‑TTTTCCCTCATCTG‑TTTTAAG‑ATGC‑G‑AAACATGAAATTC‑AACGGA‑GGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGAc                                       <  1:267949/151‑1 (MQ=255)                                                                                                          acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCCG‑AAAAAACGAACAGTA‑CACTTTTT‑TTTTCCCCCTTTTT‑TTTTAAA‑AGCG‑G‑AAAAATAAAATTA‑AAAGAG‑GGAAAAAA‑AAAGGGGATTAATTGG‑TATTTTTTAGaaaattaaacaaa                                  >  2:140607/1‑141 (MQ=255)                                                                                                             gTAAATAATAAAAATTTCTCGAGGCGGGGCCAG‑AAAAAAAGGCAAGCA‑CA‑ATTTATTTTTCCCCCACTTT‑TTTTAAA‑AT‑C‑G‑AAAAAAGAAAATT‑AAAAGA‑GGAAAAATAAAATGGGAA‑AATTGGGTTTTTTTTAGAAAATTACACCGaaaa                               >  2:460866/1‑151 (MQ=255)                                                                                                                aataatAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAGCGACATGCA‑CA‑TTTTA‑TTTTCCCTCATCTG‑TTTTAAG‑ATGC‑G‑AAACATGAAATTC‑AACGGA‑GGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGAt                            <  1:5093/151‑1 (MQ=255)                                                                                                                           tttCTCGAGGCGGGTCCAG‑AAAAAGCGACATGCA‑CA‑TTTTA‑TTTTCCCTCATCTG‑TTTTAAG‑ATGC‑G‑AAACATGAAATTC‑AACGGA‑GGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGc                 <  2:523061/151‑1 (MQ=255)                                                                                                                                cGAGGCGGGTCCAG‑AAAAAGCGACATGCA‑CA‑TTTTA‑TTTTCCCTCATCTG‑TTTTAAG‑ATGC‑G‑AAACATGAAATTC‑AACGGA‑GGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAg            <  2:384862/151‑1 (MQ=255)                                                                                                                                          ccAG‑AAAAAGCGACATGCA‑CA‑TTTTA‑TTTTCCCTCATCTG‑TTTTAAG‑ATGC‑G‑AAACATGAAATTC‑AACGGA‑GGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaa  <  1:106817/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                        CTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAG‑AAAAAGCGACATGCA‑CA‑TTTTA‑TTTTCCCTCATCTG‑TTTTAAG‑ATGC‑G‑AAACATGAAATTC‑AACGGA‑GGAAAAAT‑AAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAA  >  NC_000913/2386716‑2387003

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