Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,273,440	0	T	G	49.0%	6.5 / 74.5	51	K54Q (AAA→CAA)	ychN	putative sulfur relay protein
Reads supporting (aligned to +/- strand):  ref base T (17/9);  new base G (14/11);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

TTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTCGCTCCCGTAAGGTGCGCCATTGGCAACGATCAC  >  NC_000913/1273291‑1273587                                                                                                                                                      |                                                                                                                                                    ttCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGtttt                                                                                                                                                    >  2:477291/1‑151 (MQ=255)       aaGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCcg                                                                                                                                              >  2:371749/1‑151 (MQ=255)        agagTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTGTTGCCcgc                                                                                                                                             <  2:190620/151‑1 (MQ=255)          agTACCGATTTCCGCCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCgcgc                                                                                                                                           <  1:248756/151‑1 (MQ=255)             aCCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTGTTGCCCGCGCaa                                                                                                                                         <  1:727735/150‑1 (MQ=255)                  tttCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGc                                                                                                                                   >  1:601793/1‑151 (MQ=255)                       aCCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCtgtg                                                                                                                               >  2:873025/1‑150 (MQ=255)                        ccccATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGAc                                                                                                                             >  1:359943/1‑151 (MQ=255)                            atcaatcaGAGTAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGc                                                                                                                          <  2:314540/150‑1 (MQ=255)                               aatcaGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATc                                                                                                                       >  1:245939/1‑150 (MQ=255)                                atcaGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATc                                                                                                                       >  2:503909/1‑149 (MQ=255)                                        aaGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAgg                                                                                                             >  1:549770/1‑151 (MQ=255)                                             tACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAAc                                                                                                          >  1:553493/1‑149 (MQ=255)                                             tACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACa                                                                                                         <  1:377401/150‑1 (MQ=255)                                              aCTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAAc                                                                                                          >  2:755281/1‑148 (MQ=255)                                               cTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGAc                                                                                                      >  1:513540/1‑151 (MQ=255)                                                  aTCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGc                                                                                                    >  2:389463/1‑150 (MQ=255)                                                     ccGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGc                                                                                                    >  1:559983/1‑147 (MQ=255)                                                      cgcgCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATc                                                                                               <  1:700834/151‑1 (MQ=255)                                                            cGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGAtt                                                                                         <  2:255957/151‑1 (MQ=255)                                                                    cAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCt                                                                                 >  2:284153/1‑151 (MQ=255)                                                                    cAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCt                                                                                 <  1:314494/151‑1 (MQ=255)                                                                            gCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGt                                                                         >  1:806724/1‑151 (MQ=255)                                                                                       ccGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGc                                                              <  1:56921/151‑1 (MQ=255)                                                                                           tACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGc                                                          >  2:675663/1‑151 (MQ=255)                                                                                              aTTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGc                                                       >  1:126585/1‑151 (MQ=255)                                                                                                  tGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGc                                                   >  2:668240/1‑151 (MQ=255)                                                                                                   gAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAg                                                    <  2:920985/149‑1 (MQ=255)                                                                                                      cGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGtt                                               <  1:304420/151‑1 (MQ=255)                                                                                                      cGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGt                                                >  1:10572/1‑150 (MQ=255)                                                                                                         tAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTaa                                             >  2:471167/1‑150 (MQ=255)                                                                                                          aaGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTaaa                                            <  2:485387/150‑1 (MQ=255)                                                                                                           aGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACa                                          <  2:562177/151‑1 (MQ=255)                                                                                                             gATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACa                                          >  1:719293/1‑149 (MQ=255)                                                                                                               tctcCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAgg                                       >  1:777357/1‑150 (MQ=255)                                                                                                                  ccAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGa                                      <  1:504792/148‑1 (MQ=255)                                                                                                                   cAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTc                                   >  2:111128/1‑150 (MQ=255)                                                                                                                    aGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTCGc                                 <  1:495242/151‑1 (MQ=255)                                                                                                                     gCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTCGCt                                <  2:775971/151‑1 (MQ=255)                                                                                                                     gCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTGTTGCCCGCGAAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTCCTCTCCTGCTCTCGTAACGAAATGGCCAGCCGCAAGCTGTTAAACAAGGAtt                                    >  2:107975/1‑147 (MQ=255)                                                                                                                      cATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTCGCTc                               >  2:82315/1‑151 (MQ=255)                                                                                                                          tGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTCGCTCCCGt                           >  2:688375/1‑151 (MQ=255)                                                                                                                                  gttgtAGCCTTCCCCTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTCGCTCCCGTAAGGTgc                    <  1:32512/150‑1 (MQ=255)                                                                                                                                    tgtAGCCTTCCCCTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTCGCTCCCGTAAGGTGCGcc                 >  2:125345/1‑151 (MQ=255)                                                                                                                                     gtAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTCGCTCCCGTAAGGTGCGcc                 >  1:34335/1‑150 (MQ=255)                                                                                                                                         ccTTCCCCTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTCGCTCCCGTAAGGTGCGCCATTgg            >  2:106620/1‑151 (MQ=255)                                                                                                                                             ccccTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTCGCTCCCGTAAGGTGCGCCATTGGCAAc        >  2:344066/1‑151 (MQ=255)                                                                                                                                             ccccTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTCGCTCCCGTAAGGTGCGCCATTGGCAAc        >  2:345904/1‑151 (MQ=255)                                                                                                                                              cccTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTCGCTCCCGTAAGGTGCGCCATTGGCAACg       <  1:665042/151‑1 (MQ=255)                                                                                                                                               ccTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTCGCTCCCGTAAGGTGCGCCATTGGCAACGa      <  1:779972/151‑1 (MQ=255)                                                                                                                                                cTGGTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTCGCTCCCGTAAGGTGCGCCATTGGCAACGAt     >  1:569521/1‑151 (MQ=255)                                                                                                                                                  ggTTGTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTCGCTCCCGTAAGGTGCGCCATTGGCAACGATCa   <  1:279649/151‑1 (MQ=255)                                                                                                                                                   gttgttGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTCGCTCCCGTAAGGTGCGCCATTGGCAACGATCAc  >  1:775348/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                    TTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAATGGCCAGCCGCAAGCTGTTAAACAAGGATTCGCTCCCGTAAGGTGCGCCATTGGCAACGATCAC  >  NC_000913/1273291‑1273587

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