Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,245,234	0	G	T	44.9%	57.9 / 23.1	49	L25F (TTG→TTT)	ycgY	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (8/19);  new base T (22/0);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 1.28e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTG  >  NC_000913/1245086‑1245381                                                                                                                                                      |                                                                                                                                                     tCGATATCCGATAAGGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGTGGCGGATAGTGAGGC‑GTTTTTgg                                                                                                                                                     <  1:76868/150‑1 (MQ=255) tCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTg                                                                                                                                                      <  2:529741/149‑1 (MQ=255)     tatCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGGTgcg                                                                                                                                                >  2:595795/1‑151 (MQ=255)         cGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGGTGCGGtt                                                                                                                                             >  2:338922/1‑150 (MQ=255)           aTAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGCGGGTTTTGGcgcgcctttt                                                                                                                                           >  2:161249/1‑142 (MQ=255)             aaaGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGGGGCGCTtttt                                                                                                                                          >  1:589220/1‑149 (MQ=255)              aaGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGaa                                                                                                                                       <  2:794415/151‑1 (MQ=255)                gTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACa                                                                                                                                     <  1:209058/151‑1 (MQ=255)                       cctggTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGGTGCGCTTTTTTAACAtt‑t                                                                                                                                 >  1:423523/4‑147 (MQ=255)                        ccGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGCTGCGCTTATTTAACATT‑TAGaa                                                                                                                             >  1:495814/1‑150 (MQ=255)                        ccGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGaa                                                                                                                             <  1:605881/150‑1 (MQ=255)                         cGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAAtt                                                                                                                           <  1:251521/151‑1 (MQ=255)                                 acaacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATc                                                                                                                   <  2:611336/151‑1 (MQ=255)                                 acaacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATc                                                                                                                   <  1:486526/151‑1 (MQ=255)                                   aacaGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGCTGCGCTTTTTTAAAAtt‑ttggaattattaa                                                                                                                     >  1:250474/1‑135 (MQ=255)                                     caGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATCCTTc                                                                                                               <  1:392081/151‑1 (MQ=255)                                      aGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGCTGCGCTTATTTTAAATT‑TTGAATTGATCAAAccct                                                                                                                >  1:508720/1‑147 (MQ=255)                                        aaTCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGCTGCGCGTATTGGAAATT‑TAGAATTGATTAAACCTTCTgg                                                                                                            >  1:730707/1‑151 (MQ=255)                                         aTCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGGTTCGCTTTTTTAACATT‑TAGAATTTATTCAACCTTCTg                                                                                                             >  1:337259/1‑149 (MQ=255)                                         aTCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTTTTGAACATT‑TAGAATTGATTAATCCTTCttg                                                                                                            >  1:9953/1‑148 (MQ=255)                                          tCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGGTGCGCTTATTGAACATT‑TTTGATTGATTAATCCtt                                                                                                                >  2:511957/1‑145 (MQ=255)                                          tCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGCTGCGCTTATTGAACATT‑TAGAATTGATTAATCCTTCt                                                                                                              >  2:533458/1‑147 (MQ=255)                                                   aaTTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATCCTTCTGGCAGATTAAGc                                                                                                  <  1:595755/150‑1 (MQ=255)                                                    attattGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATCCTTCTTGCAGATTAAGcc                                                                                                 >  2:641819/1‑150 (MQ=255)                                                            ttCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATCCTTCTGGCAGATTAAGCC‑AAATAccc                                                                                        <  1:788192/150‑1 (MQ=255)                                                               gttgttTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGCTGCGCTTATTTAACATT‑TAGAATTTATCAATCCTTCTTGCAGATTAAGCCAAAATACCCGtt                                                                                     >  2:462343/1‑151 (MQ=255)                                                               gttgttTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGCTGCGCTTATTTAACATT‑TAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATAACCGtt                                                                                     >  2:183665/1‑151 (MQ=255)                                                               gttgttTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGtt                                                                                     <  1:435265/151‑1 (MQ=255)                                                                  gttTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGCTGCGCTTATTGAACATT‑TTGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCgg                                                                                  >  2:285670/1‑151 (MQ=255)                                                                    ttttGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCGTTTTTAAAATT‑TAGAATTGATCAAATCTTCTGGTAGATTAAGTCAAAAAACCCCTTCgg                                                                                  >  2:134261/1‑149 (MQ=255)                                                                       tGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAAt                                                                              <  2:221950/150‑1 (MQ=255)                                                                           gtgATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTc                                                                         <  1:641763/151‑1 (MQ=255)                                                                             gATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCa                                                                        <  2:190339/150‑1 (MQ=255)                                                                             gATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAg                                                                       <  2:8110/151‑1 (MQ=255)                                                                              aTTCAGAAAATAAATGTTTTATACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGTGTGGCATTTTGGACATTGTGGAATTG‑TATATTCTTCTGGGCGTTTAAGCCGAAATACCCGTTCgggagagggtcagt                                                                      >  1:571337/1‑139 (MQ=255)                                                                                 cAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTa                                                                    <  2:398450/150‑1 (MQ=255)                                                                                  aGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGGTGCGCTTTTTTAACATT‑TAGAATTGATCAATCCTTCTTGCAGATTAAAGCAAAATAACCGGTCGGCAATTCTTTAGtt                                                                     >  1:245362/1‑148 (MQ=255)                                                                                   gAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTAcgc                                                                 <  2:495863/151‑1 (MQ=255)                                                                                      aataaataTTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGaa                                                              <  1:480591/151‑1 (MQ=255)                                                                                       ataaataTTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGCTGCGCCTTTTTAACATT‑TAGAATTTATCAATCCCTTTGGCCGATTAAGCCAAAAAACCCCTTCGGCAATGCTTCAGTTAcgcg                                                                >  2:796324/1‑148 (MQ=255)                                                                                        taaataTTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TTGAATTTATCAATCCTTCTGGCAGATTAAGCCAAAAAAACCGTTCGGCAATGCTTTAGTTAcgcg                                                                >  2:246486/1‑147 (MQ=255)                                                                                                 tAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGGTGCGCTTATTTAACATT‑TAGAATTGATCAATCCTTCTTGCAGATTAAGCCAAAAAACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATgg                                                     >  2:607994/1‑149 (MQ=255)                                                                                                  aGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGGTGCGCTTATTTAACATT‑TAGGAATGATCAATCCTTCTGGCAGATTAAACCAAAATACCCGGTCCGCAATGCTTCAGTTTCCCGAAGAGGGAAGGt                                                    >  2:502381/1‑149 (MQ=255)                                                                                                  aGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGCTGCGCTTTTTTTACATT‑TTGAAATTATCAAACCTTCTGGCAGATTAAACCAAAAAAACCCGTCGGGAATGGGTCCGTTACGCGGAGAGGAATGGTCt                                                  >  1:702038/1‑151 (MQ=255)                                                                                                           ccAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCttt                                         <  1:278075/151‑1 (MQ=255)                                                                                                             aGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTTTTGAACATT‑TAGAATTGATCAAACCTTCTGGCAGATTAAGCCAAAATACCCCTTCGGCAATGCTTCAGGTACGCGGAGAGGAATGGTCTGAAtttttttt                                       >  2:646983/1‑151 (MQ=255)                                                                                                              gTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGGTGCGCTTATTTAAAATT‑TTGAATTTATCAATCCTTCTGGCAGATTAAACCAAAAAAACCGTTCGGCAATGGTTCAGTTAAGCGGAGAGGAATGGTCTGAATtttttttt                                      >  2:308429/1‑151 (MQ=255)                                                                                                                cAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGCTGCGCTTATTTAAAAAT‑TAGAAATTATTAAACCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGGAATGCTTCATTTACGCGAAGAGGAATGGTCTGAATTCttttttt                                     >  1:74847/1‑150 (MQ=255)                                                                                                                cAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATTCCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTg                                    <  1:705873/151‑1 (MQ=255)                                                                                                                 aGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTTGCTGCGCTTTTTTAACATT‑TTGAATTGATCAATCCTTCTGGCAGATTAAGCCCAAAAACCCGCTCGGCAAATCTTCAGTTTCGCGGAGAGGGAAGGGCTGAAtttttttttt                                     >  1:626276/1‑149 (MQ=255)                                                                                                                                                    ttGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTAtgtg  >  1:546062/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCGTTAAAAATGGGGCGGATAGTGAGGC‑GGTTTTGGCTGCGCTTATTGAACATT‑TAGAATTGATCAATCCTTCTGGCAGATTAAGCCAAAATACCCGTTCGGCAATGCTTCAGTTACGCGAAGAGGAATGGTCTGAATTCTTTTTTTGGTTGCTAAATTCTCTGGAGTGTCTGGATTATGTG  >  NC_000913/1245086‑1245381

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