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

CGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAT‑GCA‑CATTTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATG‑AAATTC‑AAC‑GGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCA‑GAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAA  >  NC_000913/2386724‑2387003                                                                                                                                           |                                                                                                                                                     cGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAT‑GCA‑Cat                                                                                                                                        <  2:492897/151‑1 (MQ=255)       aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGGAAAAACGACAA‑GCA‑Caattattt                                                                                                                                  >  2:49222/1‑147 (MQ=255)       aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAGGACAA‑GGA‑Caatttttt                                                                                                                                  >  2:81371/1‑151 (MQ=255)         aCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGACAT‑GGA‑Caatttttttt                                                                                                                                >  2:22930/1‑151 (MQ=255)            aTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGCCAC‑GCG‑Cattttatttttc                                                                                                                              >  2:672424/1‑148 (MQ=255)                          gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCACAAAAAAAGACAT‑GAA‑CATTTTATTTTTTCTCATT‑TGtttttaa                                                                                                              >  2:275307/1‑148 (MQ=255)                           tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAT‑GCA‑CATTTTATTTTCCCTCATC‑TGTTTTAAGa                                                                                                             <  2:678457/151‑1 (MQ=255)                               aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCTGAAAAAAGAAAT‑GGA‑CATTTTATTTTTTCTTATT‑TGttttttaatgtg                                                                                                         >  1:216985/1‑143 (MQ=255)                               aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACAT‑GCA‑CCTTTTATTTTTCCTCATC‑TGTTTTTAAATGCg                                                                                                         >  2:127908/1‑151 (MQ=255)                                                      gAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAT‑GCA‑CATTTTTTTTTTCCTCCTC‑TTTTTTTAGATGCGGAAAAAG‑AAAATT‑AAA‑GGGGGaa                                                                               >  2:283458/1‑151 (MQ=255)                                                            ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAA‑GCA‑CATTTTTTTTTTCCCCCTC‑TTTTTTTAAAAGCGGAAAATGAAAATTC‑ACC‑GGGGGaaaaaaa                                                                          >  2:775495/1‑149 (MQ=255)                                                                cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCGAAAAAAAGACAT‑GCA‑CAATTTTTTTTTCCTCATC‑TGTTTTTTGAAGCGAAAAAAG‑AAAATC‑AAA‑GGAGGGAAAATAAATgg                                                                     >  1:483993/1‑151 (MQ=255)                                                                  ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCACAAAAAAAGAAAA‑GAA‑CATTTTTTTTTTCCTCATC‑TTTTTTTAGAAGCGAAAAAAT‑AAAATC‑AAA‑GGAGGAaaaaaaaaaagga                                                                   >  2:596523/1‑146 (MQ=255)                                                                  ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGAAAA‑GCA‑CATTTTTTTTTTTCTCCTCTTTTTTTAAAATGGGAAA‑AAG‑AAATTT‑AAC‑GGGGGGAAAAAAAATgggg                                                                   >  2:287623/1‑150 (MQ=255)                                                                      gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCCAAAAAAACGCATGGAA‑CTTTTTTTTTCCCCTCCT‑‑TGTTTTTATAGGCGCAACATA‑AAATTCAAAC‑AGGGGAAAAA‑AAATggggaggt                                                               >  2:800736/1‑146 (MQ=255)                                                                        ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGACAA‑GAA‑AAATTTATTTTTTCTCATC‑TGTTTTTAGAAGAGAAAAAAG‑AAAATT‑AAA‑AGAGGGaaaaaaaaatagaaagatt                                                             >  1:472563/1‑140 (MQ=255)                                                                        ggATGGCGTTATAACTTTGTTATACAACAGTAAATAATAAAAATCTCTCGAGGCGGGTTCACTAAAAGAGACAT‑GCA‑CACTTTATTTTTACTCATC‑TGTTTTAAGAAAGGAAACACA‑CAAACC‑CAC‑GGAGGAAAAATAAATTGGAATGAtt                                                             >  1:796598/1‑151 (MQ=255)                                                                          aTGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAAT‑ACA‑AAATTTATTTTTTTTCATC‑TGTTTTTAAAATCGAAACAAG‑AAAAAC‑AACAGGGGGAAAAAAAAATGGGatgagttg                                                            >  1:711992/1‑144 (MQ=255)                                                                             gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGTTCCAGCAAAAAAGACAT‑AAA‑AATTTTATTTTTTCTCCTC‑TTTTTTATAGTGCGAAAAAAG‑AAATTC‑AAA‑GGAGGGAaaaaaaattgtgtagatttgcta                                                        >  1:756028/1‑136 (MQ=255)                                                                             gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACAT‑GCA‑CATTTTTTTTTTCCTCCTT‑TTTTTTTAAATGCGAAAAAAG‑AAAATT‑AAA‑AGAGGGAAAATAAATGGGGATGGTTGGCTa                                                        >  2:435452/1‑151 (MQ=255)                                                                             gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAT‑GGACCATTTTTTTTTCCC‑CCTC‑TTTTTTTAGATGGGGAAAAAT‑AAAATC‑CAC‑GGGGGGAAAAAAAATAGGGAAGAtttggtt                                                        >  2:638380/1‑146 (MQ=255)                                                                             gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAA‑GCA‑CCTTTTTTTTTTCCTCCTC‑TTTTTTTAAATGCGGAAAAAG‑AAAATT‑AAA‑GGAGGAAAAAAAAATGGGGATGAtttggtt                                                        >  2:170652/1‑146 (MQ=255)                                                                             gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAA‑GCA‑CATTTTTTTTTTCCTCCTC‑TTTTTTTAAAAGCGAAAAAAGAAAATTA‑AAG‑GGGGGAAAAAAAAAAGGGaagatttggtt                                                         >  2:145284/1‑142 (MQ=255)                                                                             gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAAA‑GAA‑AAATTTTTTTTTCCTCCTCTTGTTTTAA‑AAGAGAAAAAAGAAAATTA‑AAG‑GGGGGAAAAAAAAATGGGaagatttggtt                                                         >  1:304567/1‑142 (MQ=255)                                                                              tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAT‑GCA‑CATTTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATG‑AAATTC‑AAC‑GGAGGAAAAATAAATGGGAATGATTGGCTa                                                        <  2:715704/150‑1 (MQ=255)                                                                                ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAT‑GCA‑CATTTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATG‑AAATTC‑AAC‑GGAGGAAAAATAAATGGGAATGATTGGCTAttt                                                     <  2:680158/151‑1 (MQ=255)                                                                                          gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGGAAAAAAGACAT‑TGA‑CATTTTTTTTTTTCTCCTC‑TGTTTTTAGAATCGAAAAAAT‑AAAATC‑AAC‑GGTGGGAAAAAAAATGGGAATGATTTGGTATTTTTTA‑Gaaaat                                          >  1:632718/1‑149 (MQ=255)                                                                                           ttATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAT‑GCA‑CATTTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATG‑AAATTC‑AAC‑GGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCA‑GAAAtt                                          <  2:739001/150‑1 (MQ=255)                                                                                           ttATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAT‑GCA‑CATTTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATG‑AAATTC‑AAC‑GGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCA‑GAAATTg                                         <  1:583767/151‑1 (MQ=255)                                                                                             ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACAA‑GGA‑CATTTTTTTTTTCCTCCTC‑TTTTTTAAAAAGGGGAAAAAG‑AAAATT‑AAA‑GGAGGGAAAAAAAATGGGGATGATTTGGTATTTTTCA‑GAAAATGa                                        >  2:489581/1‑150 (MQ=255)                                                                                                caacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAAT‑GGA‑CACTTTTTTTTTTCTCATC‑TGTTTTTTGATGTGAAAAAAG‑AAAATC‑AAA‑GGAGGGAAAATAAATAGGAAAAAATGGGTATTTTTCAGGAAAttaatca                                     >  1:713571/1‑146 (MQ=255)                                                                                                         ataataAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAT‑GCA‑CATTTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATG‑AAATTC‑AAC‑GGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCA‑GAAATTGACTCAGAAAAGAtt                           <  1:325716/151‑1 (MQ=255)                                                                                                            ataaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAT‑GCA‑CATTTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATG‑AAATTC‑AAC‑GGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCA‑GAAATTGACTCAGAAAAGATTAAc                        <  1:340023/151‑1 (MQ=255)                                                                                                              aaaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACAT‑GCA‑CATTTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATG‑AAATTC‑AAC‑GGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCA‑GAAATTGACTCAGAAAAGATTAACCa                      <  2:339246/151‑1 (MQ=255)                                                                                                                          aGGCGGGTCCAGAAAAAGCGACAT‑GCA‑CATTTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATG‑AAATTC‑AAC‑GGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCA‑GAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGt           <  1:757941/150‑1 (MQ=255)                                                                                                                            gCGGGTCCAGAAAAAGCGACAT‑GCA‑CATTTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATG‑AAATTC‑AAC‑GGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCA‑GAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc          <  1:283346/149‑1 (MQ=255)                                                                                                                               ggTCCAGAAAAAGCGACAT‑GCA‑CATTTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATG‑AAATTC‑AAC‑GGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCA‑GAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa        >  1:729604/1‑148 (MQ=255)                                                                                                                               ggTCCAGAAAAAGCGACAT‑GCA‑CATTTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATG‑AAATTC‑AAC‑GGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCA‑GAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACt      <  2:346951/150‑1 (MQ=255)                                                                                                                                  ccAGAAAAAGCGACAT‑GCA‑CATTTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATG‑AAATTC‑AAC‑GGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCA‑GAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaa  <  2:786702/151‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                     CGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAT‑GCA‑CATTTTATTTTCCCTCATC‑TGTTTTAAGATGCGAAACATG‑AAATTC‑AAC‑GGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCA‑GAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAA  >  NC_000913/2386724‑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: ‑