Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,386,862	0	G	A	54.8%	0.1 / 52.1	42	intergenic (+31/‑72)	yfbL/yfbM	putative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (3/16);  new base A (23/0);  total (26/16)
Fisher's exact test for biased strand distribution p-value = 5.82e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-07

GCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGC‑G‑AAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAC  >  NC_000913/2386715‑2387006                                                                                                                                                    |                                                                                                                                                    gCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGa                                                                                                                                                  >  1:297994/1‑151 (MQ=255)        aTCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGcac                                                                                                                                           <  2:94743/151‑1 (MQ=255)           gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACaagaaaatt                                                                                                                                        >  1:478246/1‑143 (MQ=255)                aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCCGAAAAAACGACATGCACattttttt                                                                                                                                   >  1:478127/1‑151 (MQ=255)                               cccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGttt                                                                                                                    <  1:478399/151‑1 (MQ=255)                                ccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGtttt                                                                                                                   <  1:627884/151‑1 (MQ=255)                                   gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACATTTTTTTTTTTCTCCTCTGtttttaa                                                                                                                >  1:507655/1‑148 (MQ=255)                                   gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGACACGCACATTTTTTTTTTCCCCaagttgttttag                                                                                                                >  2:167872/1‑140 (MQ=255)                                    tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa                                                                                                               <  2:668756/151‑1 (MQ=255)                                    tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa                                                                                                               <  2:189585/151‑1 (MQ=255)                                                   aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATTGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAAGTTCTCGAGGCGGGTCCAGAAAACGCGACATGCACATTTTATTTTCCCTCAGCTGTTTTAAGATGC‑G‑AAACATG‑AAAt                                                                                             <  1:400189/151‑1 (MQ=255)                                                    ttGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAGCGACATGCACATTTTATTTTCCCACATCTGTTTTAAGAGGC‑G‑AAACATG‑AAAtt                                                                                            <  1:68941/151‑1 (MQ=255)                                                          tATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGC‑G‑AAACATG‑AAATTCAACgg                                                                                      <  1:400928/151‑1 (MQ=255)                                                           aTCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGGAAAAAAGAAATGCACATTTTTTTTTTCCTCCTCTGTTTTTAGAAGC‑G‑GAAAAAG‑aaaattaaaggg                                                                                     >  1:382278/1‑142 (MQ=255)                                                                     ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAACACATGCACATTTTATTTTTTCTCATCTGTTTTTAGATGG‑GAAAAAAGG‑AAATTAACTGGGGGaaaaaaa                                                                            >  1:229488/1‑149 (MQ=255)                                                                        tCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGC‑G‑AAACATG‑AAATTCAACGGAGGAAAAATAAATg                                                                        <  1:21867/151‑1 (MQ=255)                                                                           ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGACAAGAACATTTTTTTTTTTCTCCTCTGTTTTTAAAAGC‑G‑AAAAAAG‑AAATTTTACGGAGGGAAAAAAAATgggg                                                                     >  1:54052/1‑150 (MQ=255)                                                                           ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACATTTTTTTTTTCCTCATCTTTTTTTAAATGC‑G‑AAACAAG‑AAAATCAACGGAGGGAAAAAAAATgggg                                                                     >  2:61951/1‑150 (MQ=255)                                                                           ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGAAAAGCAAATTTTATTTTTTCTCATCTCGTTTTAGATGA‑G‑AAACAAG‑AAAATTAAAGGGGGGAaaataaaagtgg                                                                     >  1:173093/1‑146 (MQ=255)                                                                           ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGAAAAGCACATTTTTTTTTTCCTCATCTGTTTTTAGATGG‑G‑GAAAATG‑AAAATTAACGGAGGAAAAAAAAATgggg                                                                     >  2:591741/1‑150 (MQ=255)                                                                           ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGACGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGC‑G‑AAACATG‑AAATTCAACGGAGGAAAAATAAATGGGa                                                                     >  2:214665/1‑151 (MQ=255)                                                                            gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCAAAAAAACGACATGCACAATTTATTTTTCCTCTTCTGTTTTTAGAAGC‑G‑GAAAATG‑AAAATTAACAGGGGGAAAAAAAATgggta                                                                    >  2:115832/1‑149 (MQ=255)                                                                               gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGACATGCACATTTTTTTTTTCCTCATCTGTTTTTAGATGG‑GAAAAAAGG‑AAATTCAAAGAGGGAAAAAAAAAGGGGatga                                                                  >  1:456759/1‑148 (MQ=255)                                                                               gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACATGGACCTTTTTTTTTTCCTCATCTGTTTTTAAAAGC‑G‑AAAAAAG‑AAAATCAACGGAGGAAAAAAAAAAGGGAATGa                                                                 >  2:431044/1‑151 (MQ=255)                                                                                 ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGC‑G‑AAACATG‑AAAt                                                                                             <  1:344470/121‑1 (MQ=255)                                                                                      gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGACAAACACATTTTTTTTTTTCTCATCTTTTTTAAAAAGC‑G‑AAACATT‑AAAATTAAAAGAGGGAAAAAAAATGGGAATGAtttgttt                                                          >  1:40044/1‑146 (MQ=255)                                                                                      gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAACACAAGAACATTTTATTTTTTCTCATCTGTTTTTAAATGA‑G‑AAAAAAG‑AAAATTAAAGGAGGAAAAATAAATAGGGATGATTGGGTa                                                          >  1:42286/1‑151 (MQ=255)                                                                                      gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCTGCAAAAAAAGAAAATCACAATTTTTTTTTTCTCATTTTTTTTTAAAAGCGG‑AAAAAGA‑AAATTCAA‑TGAGGGAAAAAAAATTGGGATGAATGGCtt                                                          >  1:576800/1‑150 (MQ=255)                                                                                                   gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGAACATTTTATTTTTCCTCATCTGTTTTTAGATGG‑G‑AAAAAAGAAAATTCAA‑GGAGGAAAAAAAAATGGGGATGATTTGGCAtttttttgaaaat                                             >  1:262923/1‑142 (MQ=255)                                                                                                   gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACATTTTTTTTTTTCTCCTCTGTTTTTAGAAGA‑G‑AAAAATG‑AAAATTAAAGGAGGGAAAAAAAATGGGAATGATTGGCTATTTTTCAGaaaat                                             >  1:662881/1‑149 (MQ=255)                                                                                                    ttATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGGAAAAACGACAAGCACATTTTATTTTTCCTCCTCTGTTTTTAGATTG‑G‑AAAAATG‑AAATTTAAAGGAGGAAAAAAAAATTGGGATTATTGGGTATTTTTGAGAAAATg                                            >  1:522770/1‑151 (MQ=255)                                                                                                           acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCATCTGTTTTTAGATGC‑G‑AAACATG‑AAATTTAACGGGGGAAAAAAAAATGGGGATGATTGGCTATTTTTCAGAAATTTAATCAGa                                     >  2:333910/1‑151 (MQ=255)                                                                                                           acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGAAAATTTTATTTTTTCTCATCTGGTTTTAGATGG‑G‑GAAAATG‑GAATTTAAAGGAGGGAAAAAAAAATGGAATGATTTGATTTTTTTTAGAAAATGAATCAGa                                     >  1:592802/1‑151 (MQ=255)                                                                                                              gTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa                                                                                                               <  1:239102/77‑1 (MQ=255)                                                                                                              gTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTTTTTTTCCTTCTCTGTTTTTAGa                                                                                                               >  2:239101/1‑77 (MQ=255)                                                                                                                     ataaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGC‑G‑AAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAAc                           <  2:416222/151‑1 (MQ=255)                                                                                                                      taaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGC‑G‑AAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAAcc                          <  1:312350/151‑1 (MQ=255)                                                                                                                                     gCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGC‑G‑AAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc             <  2:9626/149‑1 (MQ=255)                                                                                                                                        ggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGC‑G‑AAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa           >  2:364933/1‑148 (MQ=255)                                                                                                                                        ggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGC‑G‑AAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa           >  2:130656/1‑148 (MQ=255)                                                                                                                                         gTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGC‑G‑AAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTg        <  2:192761/150‑1 (MQ=255)                                                                                                                                               aaaaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGC‑G‑AAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAc  <  1:61952/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGC‑G‑AAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAC  >  NC_000913/2386715‑2387006

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