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

CGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTT  >  NC_000913/2386724‑2387010                                                                                                                                            |                                                                                                                                                     cGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCCCAAAAAACGACCTGccca                                                                                                                                           >  1:80207/1‑150 (MQ=255)         aCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACattttatttt                                                                                                                                  <  2:231066/151‑1 (MQ=255)          cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GGTCCGAAAAAAAAAAAAGCACattttattttt                                                                                                                                 >  2:325788/1‑150 (MQ=255)                            gATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGGCATGCACATTTTATTTTCCCTCATCTGTTTTAAGAt                                                                                                               <  1:145239/151‑1 (MQ=255)                            gATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa                                                                                                                <  1:289284/150‑1 (MQ=255)                                 cGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTTCAGAAAAAAAGACAAGCAAAATTTATTTTTCCTCCACTTTTTTTAAATggg                                                                                                            >  2:104064/1‑147 (MQ=255)                                   cagacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAAc                                                                                                        <  1:459593/151‑1 (MQ=255)                                    agacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAAAAAAAACGACATGCAAATTTTATTTTTCCTCATCTGTTTTTAAATGCGaaa                                                                                                         >  1:261664/1‑149 (MQ=255)                                       cagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAAAGACAAGCACATTTTATTTTTCCTCATCTTTTTTTAAAAGCGaaaaa                                                                                                       >  1:329729/1‑146 (MQ=255)                                         gaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAAAGACATGCACATTTTTTTTTTCCTCCTCTGTTTTTAGATGCGAAACATGaa                                                                                                   >  1:609491/1‑150 (MQ=255)                                         gaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAAAGACATGCACATTTTATTTTTTCTCATCTATTTTTAGAATCGAAAAAAGaa                                                                                                   >  1:174509/1‑150 (MQ=255)                                            tGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAAtt                                                                                                <  2:125385/150‑1 (MQ=255)                                               aTTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAAc                                                                                            <  1:212907/151‑1 (MQ=255)                                               aTTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAAc                                                                                            <  2:297675/151‑1 (MQ=255)                                                       aaaaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGaaa                                                                                    <  1:208861/151‑1 (MQ=255)                                                                     aGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GGTCCGAAAAAAAGACAAGCACATTTTTTTTTTTCTCCTCTGTTTTTAAATGAGGAAAAAGAAAATTAAAGGAGGGAAAAAAAATagg                                                                          >  2:537971/1‑147 (MQ=255)                                                                            ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc                                                                <  1:71732/150‑1 (MQ=255)                                                                            ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGAGAATGATTGGCt                                                               <  1:325795/151‑1 (MQ=255)                                                                                ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAAAGACAAGCACATTTTATTTTTTCTCCTCTTTTTTTAGATGCGAAAAATGAAAATTAAAGGAGGAAAAATAAATGGGGATGATTTGCTAttt                                                           >  2:143174/1‑151 (MQ=255)                                                                                 tataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGTGTCC‑CAAAAAAAGACAAGAACAATTTTTTTTTCCTCATCTGTTTTTATATGCGGAAAATTAAAATTTACCGAGGAAAAAAAAATGGGGATGATTTGGTAtttt                                                          >  1:579169/1‑151 (MQ=255)                                                                                  ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTAtttt                                                          <  2:149834/150‑1 (MQ=255)                                                                                  ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAAAGACAAGCACATTTTATTTTTCCTCATCTGTTTTTAGATGCGAAAAAAGAAAATTAACGGAGGAAAAAAAAATTGGAAAGATTTGCTAttttt                                                         >  2:539300/1‑150 (MQ=255)                                                                                    aaCTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGc                                                        <  2:518097/150‑1 (MQ=255)                                                                                        ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GGCCCGAAAAAAACAGAAACACATTTTTTTTTTCCTCCTCTGTTTTTAAATTAGAAACAAGAAAATTAAAGGGGGGGAAAAAAAAGGGGATGATTTGGTTTTTTTCAGaa                                                    >  2:36853/1‑150 (MQ=255)                                                                                              taCAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCCGAAAAAAAGGCATGCACATTTTTTTTTTCCTCCTCTGTTTTTAAATGCGAAAAAAGAAAATCAAACGGGGAAAAAAAAAAGGGGATGATTGGGTATTTTTTAGaaaat                                                 >  1:54418/1‑145 (MQ=255)                                                                                                    aGTAAATAATAAAAATTTCTCGAGGCGG‑GGCCAGAAAAAAAGACATGAACATTTTATTTTTCCTCATCTGTTTTTAAAAGTGAAAAAAGAAAATTAAAGGAGGGAAAAAAAATAGGAATGATTGGCTATTTTTCAGAAATTGACTCAGaa                                        >  2:251252/1‑150 (MQ=255)                                                                                                               aaaaTTTCTCGAGGCGG‑GTTCAGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTGTTTTTAAATGCGAAAAAAGAAAATTAAAGGAGGAAAAAAAAAAAGGGAAGATTGGCTTTTTTTTAGAAAATTAATCAGAAAAAAATaa                                >  2:289435/1‑147 (MQ=255)                                                                                                                aaaTTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAAc                           <  2:428928/151‑1 (MQ=255)                                                                                                                         gAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGt                  <  1:289442/151‑1 (MQ=255)                                                                                                                           ggCGG‑GTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCt                <  2:26394/151‑1 (MQ=255)                                                                                                                                      gAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAc      >  2:60161/1‑151 (MQ=255)                                                                                                                                          aaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCttt  <  2:609474/151‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                     CGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTT  >  NC_000913/2386724‑2387010

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