Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,413,735	0	A	T	38.0%	68.3 / 30.8	51	Y2* (TAT→TAA)	ydaC	DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (1/30);  new base T (19/0);  total (20/31)
Fisher's exact test for biased strand distribution p-value = 6.58e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

CAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG  >  NC_000913/1413588‑1413885                                                                                                                                                    |                                                                                                                                                         cAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTtata                                                                                                                                                        <  1:730942/150‑1 (MQ=255)   gTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAt                                                                                                                                                     <  2:652082/151‑1 (MQ=255)      ttACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGa                                                                                                                                                   <  2:49655/150‑1 (MQ=255)      ttACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGa                                                                                                                                                   <  2:407122/150‑1 (MQ=255)        aCTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTaa                                                                                                                                                  >  2:924872/1‑149 (MQ=255)               gtcctctGGAAACGCCAGCTTCTTTTTTATCTGAGAGCATTTTTTCGCATTCTGGTTTCGTTAGTATAGATTTTGAATATCTTGTCCAGTTAGGCGGTGTGCCTCACTCCGTGCCAATAGTGGCGGTAATATCATAAATGAACAtctc                                                                                                                                            <  1:107092/145‑1 (MQ=255)               gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTc                                                                                                                                            <  1:482024/148‑1 (MQ=255)                ctctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCat                                                                                                                                         <  2:510919/150‑1 (MQ=255)                    ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAACACCCCCcttatt                                                                                                                                     >  2:495941/1‑150 (MQ=255)                      ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattatt                                                                                                                                     <  2:1036562/148‑1 (MQ=255)                        aaaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAATACCCCCACttttttttt                                                                                                                                >  2:18608/1‑150 (MQ=255)                          aCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCa                                                                                                                              <  2:799506/151‑1 (MQ=255)                           cGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCTCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAAGTTTTATATTCATACCTCCATCATTATTTcc                                                                                                                               >  1:1079366/1‑149 (MQ=255)                           cGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGAAAATTTTTTAATGGACAACCCCCTTAttttttt                                                                                                                                >  2:428139/1‑147 (MQ=255)                              cTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCATGAACACCTCCCCTATTAttttttgtg                                                                                                                           >  2:72016/1‑144 (MQ=255)                              cTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATGAAAACCCCCCTTATTTTTTTCCGTgg                                                                                                                          >  2:274836/1‑151 (MQ=255)                               tGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAAAACCCCCCTTATTTTTTTCAGTgg                                                                                                                          >  2:954617/1‑150 (MQ=255)                                gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTgg                                                                                                                          <  1:79788/149‑1 (MQ=255)                                gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGt                                                                                                                         <  2:327064/150‑1 (MQ=255)                                        tttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTAt                                                                                                                 <  2:882910/150‑1 (MQ=255)                                               gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCt                                                                                                          <  2:551512/150‑1 (MQ=255)                                                      ttttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTATTAACAACCCCCTTTTTTTTTTCAGTGGGTTGtttttttttttttt                                                                                                        >  1:1065248/1‑137 (MQ=255)                                                       tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTg                                                                                                  <  2:847482/150‑1 (MQ=255)                                                              cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGcttctg                                                                                            <  1:805444/149‑2 (MQ=255)                                                                  cTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTGATTTTTTACAAGAAAACCCCCCTTTTTTTTTTCAGCGGGTGGTTTATTCCTTCTTTTTTGTGGtttttttt                                                                                         >  2:170076/1‑148 (MQ=255)                                                                       ttCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTTATTAACACCTCCCTTATTTTTTTTACTGGGTTGTTTTTTTCATCTTTTTTGTGGttttttttt                                                                                        >  1:808361/1‑143 (MQ=255)                                                                        tCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCATTAAAACCTCCCCTTTTTTTTTCTGTGGTTCGTTTTTTTCCTCTTTTTCGTGCTTCTTTTTTCC‑TTctcc                                                                               >  2:1124610/1‑148 (MQ=255)                                                                             aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGATTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑a                                                                          <  1:706102/150‑1 (MQ=255)                                                                                ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AAcc                                                                       <  1:926381/150‑1 (MQ=255)                                                                                   gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCGGtt                                                                   <  2:1141549/151‑1 (MQ=255)                                                                                   gATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCGGtt                                                                   <  2:300756/151‑1 (MQ=255)                                                                                      tttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCGGTTCg                                                                 <  1:130301/150‑1 (MQ=255)                                                                                        tGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAACAACCCCCCTTTTTTTTTTAGTGGGTTGTTTTTTTCCTTTTTTCAGTGCTTTTTTTTTTC‑TTTCACCAT‑AAACCGtt                                                                   >  1:831093/1‑146 (MQ=255)                                                                                          aaTATCTTGTCCAGTTAGTAGGTGTGCCCCCTTCCTTTTCAATAGTGGCGGTAATTTTCTCCATCAAGACCTCCATTATTAATTCCAGTGGTTCGCTTATTCCACCTTTCGAGTGCCTCTCTTTCTG‑TTCCTCCAT‑AACCGGTTCGGGAtt                                                            <  1:895762/151‑1 (MQ=255)                                                                                          aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCGGTTCGGGAtt                                                            <  1:1054018/151‑1 (MQ=255)                                                                                          aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTCATAAACACCCCCCTTTTTTTTTTCCGTGGGTTGTTTTTTTCATTTTTTTTGTGGTTCTTTTTTTC‑TTTCCCCCT‑AACCCGGTCgggg                                                              >  2:204877/1‑148 (MQ=255)                                                                                           atatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAAAACCTCCCTTATTTTTTTCCGTGGGTTGTTTTTTTCATCTTTTTAGTGCTTTTTTTTTTCTTTCCCCCTA‑AACCGGTTggg                                                                >  2:381800/1‑147 (MQ=255)                                                                                              tCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACACTAACACCCCCCTTATTATTTTCCGTGGTTTGTTTTTTTCCTCTTTTTAGTGGTTTTTTTTTTC‑TTCCCCCCT‑AACCCGTTCgggggttgcc                                                        >  2:604004/1‑144 (MQ=255)                                                                                               cTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATGAACAACCCCCTTATTTTTTTCAGTGGGTTGTTTTTTTCCTTTTTTTGGTGGTTTTTTTTTTC‑TTTCCCCCT‑AACCCGGTCGGGGtt                                                            >  2:407048/1‑146 (MQ=255)                                                                                                 tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCGGTTCGGGATTCGCa                                                        <  2:349048/148‑1 (MQ=255)                                                                                                 tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCGGTTCGGGATTCGCATc                                                      <  2:626316/150‑1 (MQ=255)                                                                                                  gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCTCCAT‑AACCGGTTCGGGATTCGCATCCg                                                    <  2:1065646/151‑1 (MQ=255)                                                                                                    ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCGGTTCGGGATTCGCATCCGtt                                                  <  2:95575/151‑1 (MQ=255)                                                                                                            tagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTc                                           <  1:104591/150‑1 (MQ=255)                                                                                                               gAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTc                                       <  2:219883/151‑1 (MQ=255)                                                                                                                   gCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTACATGAAAACCTCCCTTATTATTTTCAGTGGGTTGGTTATTTCATCTTTTTAGGGGTTTTTTTTTAC‑TTTCACCAT‑AACCGGTTCGGGGTTTGCATCCCTTACACTTCGCCCgg                                     >  2:410649/1‑149 (MQ=255)                                                                                                                    ccACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAATAACACCCCCCTTTTTTTTTTCAGTGGGTTGTTTTTTTCCTCTTTTTAGTGCTTTTTTTTTTC‑TTTCACCATAAACCGGTTGGGGGTT‑GCATCCCTTACACTTTGccc                                       >  2:42290/1‑144 (MQ=255)                                                                                                                           ccttTTCAATAGTGGCGGTAATTTTATAAATTAACACCCCCATTATTATTTTTACTGGGTCTTTTTTTTCATTTTTCGTCTTCTTCTTTTTTTC‑TTTCCCCCC‑CACCCGTTCGGGGTTCGCATTCGTTAATTTTCTCTCCGTTatat                               >  1:219437/1‑147 (MQ=255)                                                                                                                                 cAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATgg                     <  2:425479/151‑1 (MQ=255)                                                                                                                                     aGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCtt                  <  1:215225/150‑1 (MQ=255)                                                                                                                                          cGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCg            <  2:203926/151‑1 (MQ=255)                                                                                                                                            gTAATTTTATACATTAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTGCGAGTGCTTCTTTTTCAC‑TTCCACTAT‑AACCGGTGCGGGATTCGCATCCGTTACGCTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGtt          <  1:671085/151‑1 (MQ=255)                                                                                                                                                   tataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATg   >  2:802950/1‑151 (MQ=255)                                                                                                                                                    ataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATgg  <  2:856011/151‑1 (MQ=255)                                                                                                                                                    ataCATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATgg  <  1:285574/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                         CAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG  >  NC_000913/1413588‑1413885

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: ‑