Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	231,694	0	C	T	42.6%	74.3 / 20.9	54	F191F (TTC→TTT)	yafD	endo/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base C (5/26);  new base T (23/0);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 9.29e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGC‑GTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGT  >  NC_000913/231555‑231841                                                                                                                                             |                                                                                                                                                     cTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGG‑GATTTTAATGCCTgggg                                                                                                                                           >  1:668834/1‑149 (MQ=255)  tGACACCCGCCTGTTGATGGGGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑AATTTTAATGCCttggg                                                                                                                                           >  1:351063/1‑146 (MQ=255)     cacCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GGTTTTAAAGGCCGGAGCCg                                                                                                                                        >  1:418443/1‑150 (MQ=255)       cccGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAg                                                                                                                                     <  1:875016/151‑1 (MQ=255)        ccGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGG‑AATTTTAATGCCtcg                                                                                                                                             >  2:967970/1‑140 (MQ=255)         cGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGaa                                                                                                                                   <  1:654506/151‑1 (MQ=255)           ccTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGGAGAAgg                                                                                                                                 >  2:493479/1‑151 (MQ=255)              gTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGa                                                                                                                                <  2:709877/149‑1 (MQ=255)                tGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTTAATGCCTGGGGCCGTAGAAgg                                                                                                                                 >  2:134794/1‑146 (MQ=255)                    ggtggtTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTTAATGCCTGGAGCCGGAGAAGGATGAAcgcg                                                                                                                        >  1:603606/1‑151 (MQ=255)                       ggtTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGG‑GATTTCAAAGCCTGGGGCCGGAGAAGGATGaaaccggt                                                                                                                      >  1:226732/1‑144 (MQ=255)                        gtTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTtat                                                                                                                    <  1:984312/151‑1 (MQ=255)                         ttAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGG‑GATTTTAAATCCTGGAGCCGGAGAAGGATGAACGCGTta                                                                                                                     >  1:175443/1‑149 (MQ=255)                            atatACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGG‑GATTTTAAAGCCCGGAGCCGTAGAAGGATGAAcgcggt                                                                                                                      >  1:842463/1‑143 (MQ=255)                              ataCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGc                                                                                                               <  1:769966/150‑1 (MQ=255)                              ataCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCt                                                                                                              <  2:738195/151‑1 (MQ=255)                                 cacGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCttt                                                                                                            <  2:787352/150‑1 (MQ=255)                                   cGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GGTTTTAATGCCTGGGGCCGGAGAAGGATGAACGCGGTATATTGCttt                                                                                                            >  2:438835/1‑148 (MQ=255)                                   cGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTgcg                                                                                                         <  1:528779/151‑1 (MQ=255)                                      cGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGG‑GATTTTAAAGCCTGGGGGCGGAGAAGGATGAAAGGGTTATATTGCTTTGTGCgg                                                                                                      >  2:18690/1‑151 (MQ=255)                                      cGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCgg                                                                                                      <  1:670556/151‑1 (MQ=255)                                      cGTCAAC‑TCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGG‑GATTTTAATGGCTGGGGCCGCAGAAGGGGGAAAGCGGTTTATCGCTTTGCGAggg                                                                                                     >  1:750509/1‑151 (MQ=255)                                          aaCTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTTAATGCCTGGGGCCGTAGAAGGGTGAAAGCGTTATTTCGCTTTTCGCggg                                                                                                     >  2:142722/1‑148 (MQ=255)                                           aCTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAAt                                                                                                 <  2:66665/151‑1 (MQ=255)                                           aCTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAAt                                                                                                 <  2:841197/151‑1 (MQ=255)                                              tCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTTAAAGCCCGGGGCCGGAGAAGGATGAAAGCGGTATATCGCTTTTCGCGGGgaa                                                                                                  >  2:113325/1‑147 (MQ=255)                                                   cTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGG‑GGTTTTAATGTCTGGGGCCGGAGAAGGATGAAAGCGGTTTATTGCTTTTCGCGGGGAATGTGGCTgg                                                                                         >  1:541483/1‑150 (MQ=255)                                                       gCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGG‑GATTTTAATGCCTGGAGCCGTAGAAGGGTGAACACGGTTTATTGCTTTTCGCGGGGAAAGTAGCTGCGcc                                                                                      >  1:446718/1‑150 (MQ=255)                                                       gCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTTAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATTGCTTTGCGCGGGGAAAGTCGCGGCGcc                                                                                      >  1:300890/1‑150 (MQ=255)                                                                   ataGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGG‑GATTTTAATTCCTGGGGCCGGAGAAGGATGAAAGCGGTATATTGCTTTTCGCGGGGAAAGTTGCTGCGCCAGGTGCGtttt                                                                           >  2:809558/1‑149 (MQ=255)                                                                   ataGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGttt                                                                            <  1:206523/148‑1 (MQ=255)                                                                    taGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTTAAATCCTGGGGGCGTAGAAGGGAGAAAGCGTTATATTCCTTTGCGCGGGGAAAGTCGCTGCGCCAGGGGCGGTTTTcc                                                                        >  2:27897/1‑151 (MQ=255)                                                                     aGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCg                                                                       <  2:810746/151‑1 (MQ=255)                                                                       tAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGG‑GATTTCAAATCCTGGGGCCGGAGAAAGATGAAAGCGGTTTATCGCTTTGCGCGGGGAATGTCGCCGCGCCAgggggggttttacggt                                                                     >  1:364995/1‑137 (MQ=255)                                                                        aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCgat                                                                     <  2:912535/150‑1 (MQ=255)                                                                         aGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTTAATGCCTGGGGCCGTAGAAGGATGAAAGCGTTTTATCGCTTTTCGCGGGGAAAGTCGCTGCGCCAGGGGGGGTTTTACgatg                                                                    >  2:228042/1‑150 (MQ=255)                                                                          gCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCgatgat                                                                  <  1:243564/151‑1 (MQ=255)                                                                            aGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTTAATGCCTGGAGCCGCAGAAGGATGAAAGCGGTATATCGCTTTTCGCGGGAAATGTCGCCGCGCCAGGGGCGttttttc                                                                        >  1:239899/1‑140 (MQ=255)                                                                                  ttCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTTAATTCCTGGGGCCGTAGAAGGATGAACGCGTTTTATCGCTTTGCGCGGGGAATGTCGCTGCGCCCGGTGGGGTTTTCCGATGATCAGCGcc                                                           >  2:233909/1‑150 (MQ=255)                                                                                    ccTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTTAATGCCTGGGGCCGTAGAAGGATGAAAGCGGTATATTGCTTTGCGCGGGGAAAGTCGCTGCGCCCGGGGCGGTTTTTCGATGATCAGCGcc                                                           >  2:329180/1‑148 (MQ=255)                                                                                    ccTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTc                                                        <  1:318429/151‑1 (MQ=255)                                                                                    ccTATTGGCGATCAGAGAGCTCTCCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGt                                                         <  2:195096/150‑1 (MQ=255)                                                                                             gATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGC‑GTTTg                                              <  2:351144/151‑1 (MQ=255)                                                                                                cAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGC‑GTTTGGTc                                           <  2:230069/151‑1 (MQ=255)                                                                                                   aTAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGGATTTC‑ATGTCTGGAGCCGTTGAAGGATGAAGGCGGTATATAGCTTTGGGCGGGGAATGTTGCCGCGCCAGGGGGGGTTTACCGATGATCAGTGCCGTCGGGCGGTTTggggtc                                         >  2:653919/1‑147 (MQ=255)                                                                                                      gCTCACCACAGCGGCCCGGTCATTATGGCGGGG‑GATTTTAAAGCCTGGGGCCGCAGAAGGGTGAACGCGGTTTTTCCCTTTTGGCGGGGAAATTGGCTGCGCCAGGGGGGGTTTTCCGCTGATTAGCGCCGTCCGGC‑GGTTGGGTGCCCGc                                     >  2:281809/1‑151 (MQ=255)                                                                                                                 cGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGC‑GTTTGGTCGCCCGCTCGATTTTGtt                          <  1:27893/151‑1 (MQ=255)                                                                                                                  ggCCCGGTCATTATGGCGGGA‑GATTTCAATTCCTGGGGCCGGAGAAGGGTGAACGCGGTATATCGCTTTGTGCGGGAAATGGCGCTGCGCCCGGGGCGGTTTTCCGATGATCAGCGCCCTCGGGC‑GGTTGGTCGCCCCCCCGGTTTTGttt                         >  1:250116/1‑151 (MQ=255)                                                                                                                     ccGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGC‑GTTTGGTCGCCCGCTCGATTTTGTTTTCt                      <  2:442479/151‑1 (MQ=255)                                                                                                                                     ga‑gaTTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGC‑GTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGt      <  2:609321/151‑1 (MQ=255)                                                                                                                                        gaTTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGC‑GTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCa    <  2:744695/151‑1 (MQ=255)                                                                                                                                        gaTTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGC‑GTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCa    <  1:146752/151‑1 (MQ=255)                                                                                                                                           ttCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGC‑GTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGt  <  1:218245/150‑1 (MQ=255)                                                                                                                                           ttCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGC‑GTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGt  <  1:583045/150‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                     CTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGA‑GATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGC‑GTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGT  >  NC_000913/231555‑231841

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