Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,386,862	0	G	A	42.9%	28.2 / 18.4	28	intergenic (+31/‑72)	yfbL/yfbM	putative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/16);  new base A (12/0);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 3.29e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAA  >  NC_000913/2386718‑2387005                                                                                                                                                 |                                                                                                                                                  ttttATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACa                                                                                                                                              >  2:12270/1‑150 (MQ=255)                    aCCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAGGACATGCACatttttttttttc                                                                                                                            >  2:255238/1‑146 (MQ=255)                           gCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGCAAAAAAGAAAAGCACATTTTTTTTTTCCTCATCTGtt                                                                                                                   >  1:354369/1‑151 (MQ=255)                                   aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAACGACATGGACAATTTTTTTTTCCTCCTCTGtttt                                                                                                                 >  1:14135/1‑145 (MQ=255)                                            acagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGACATGCACATTTTATTTTTTCTCATCTGTTTTTAA‑ATGCGaaaatcg                                                                                                 >  1:127459/1‑147 (MQ=255)                                            acagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGa                                                                                                       <  2:127496/145‑1 (MQ=255)                                            acagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATGGCCACTTTATTTTTCCTCATTTGTTTTTAA‑ATGCGaa                                                                                                      >  2:242453/1‑146 (MQ=255)                                            acagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATGCACATTTTTTTTTTCCTCATCTTTTTTTAA‑AAGCGaaaaa                                                                                                   >  2:432584/1‑147 (MQ=255)                                              agaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCACAAAAAAAGACATGCACATTTTTTTTTTTCTCATCTGTTTTTAA‑ATGCGaaaaaaga                                                                                                >  2:489119/1‑145 (MQ=255)                                                aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGAAACATGAAAt                                                                                             <  2:186841/151‑1 (MQ=255)                                                    aaTTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGAAACATGAAATTCaa                                                                                         <  1:302780/151‑1 (MQ=255)                                                        aTCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGAAACATGAAATTCAACgg                                                                                      <  2:265686/150‑1 (MQ=255)                                                                aTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGAAACATGAAATTCAACGGA‑GGAAaaat                                                                            <  1:175086/151‑1 (MQ=255)                                                                  ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGAAACATGAAATTCAACGGA‑GGAAaaataa                                                                          <  2:94721/151‑1 (MQ=255)                                                                   gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGAAACATGAAATTCAACGGA‑GGAAaaataaa                                                                         <  1:310082/151‑1 (MQ=255)                                                                      cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATgg                                                                      <  2:14137/151‑1 (MQ=255)                                                                       aGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCACAAAAAAAGACAAGCAAATTTTTTTTTTTCTCATCTGTTTTTAT‑ATGCGGAACATGAAAATTAAAGGA‑GGAAAAATAAATgg                                                                      >  1:8509/1‑150 (MQ=255)                                                                            gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGCAATGAACATTTTTTTTTTCCTCATCTGTTTTTAA‑ATGGGGAAAAAAAAATTTAAAGGA‑GGGAAAAAAAATggg                                                                     >  2:221465/1‑146 (MQ=255)                                                                                        ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTAtttt                                                     <  2:398145/150‑1 (MQ=255)                                                                                            cTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGGAAAAAAGAAATGCCAAATTTTTTTTTCCTCCTCTTTTTTTAG‑ATGAGAAAAAAGAAATTAAAAGGAGGGAAAAAAAAAAGGGAAGGA‑TGGGTAtttt                                                     >  2:350514/1‑146 (MQ=255)                                                                                              ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGa                                                <  2:383711/149‑1 (MQ=255)                                                                                                gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAtt                                            <  2:409840/151‑1 (MQ=255)                                                                                                       aacaGTAAATAATAAAAATTTCTCGAGGCGGGGTCCGAAAAAAAGACATGAAAATTTTATTTTTCCCCCTTTGTTTTTAGAATGGAAAAAATAAAATT‑AAAGGA‑GGGAAAAAAAATTGGAATTATTGGGTAttttttaaaaaattactc                                       >  1:524632/1‑135 (MQ=255)                                                                                                             aaaTAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAg                                <  1:102799/150‑1 (MQ=255)                                                                                                                              cGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTg                <  2:27480/149‑1 (MQ=255)                                                                                                                              cGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGa               <  1:161505/150‑1 (MQ=255)                                                                                                                                  gCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa          <  2:451021/151‑1 (MQ=255)                                                                                                                                          aGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaaaa  <  2:210986/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                  TTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAG‑ATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAA  >  NC_000913/2386718‑2387005

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