Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,008,480	+A	coding (637/654 nt)	yedK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,008,476	1	.	A	92.1%	215.2 / 13.9	76	coding (633/654 nt)	yedK	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (5/1);  new base A (31/39);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 9.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TT‑A‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGT  >  NC_000913/2008401‑2008545                                                                              |                                                                         aGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTTATG‑tt                                                                           <  2:1562048/76‑1 (MQ=255) aGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑tt                                                                           <  1:4053425/76‑1 (MQ=255) aGAAATCGCGGCCAGTGGCTGGGTTCCAGCAAACCCAGTCTCCTCGCACCCCGTATCGCGCGCGGGGGGGGATG‑tt                                                                           >  1:2877188/1‑76 (MQ=255)  gAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TT‑a                                                                         <  1:2932651/76‑1 (MQ=255)   aaaTCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TT‑a‑a                                                                       <  2:4307867/76‑1 (MQ=255)    aaTCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TT‑a‑aa                                                                      <  1:4364116/76‑1 (MQ=255)       cgcgGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TT‑a‑aa‑aac                                                                  <  2:2526369/76‑3 (MQ=255)        gcgGCCAGTGGCTGTGTTCCAGCAAACCCGTTCTCCTGGCACCCGGTATCGCGCGCGGGGGGTAATG‑TT‑a‑aa‑aacc                                                                 >  2:3961133/1‑73 (MQ=255)        gcgGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TT‑a‑aa‑aacc                                                                 <  2:2909017/76‑4 (MQ=255)          ggCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTT‑A‑AA‑AC‑CAg                                                               >  1:3066329/1‑76 (MQ=255)           gCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAgg                                                              <  1:3873096/76‑1 (MQ=255)           gCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGGAATG‑TTAA‑AA‑AC‑CAgg                                                              >  1:2092846/1‑76 (MQ=255)            ccAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAggg                                                             <  1:2749310/76‑1 (MQ=255)            ccAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAggg                                                             <  1:1981740/76‑1 (MQ=255)             cAGTGGGTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGgc                                                            <  2:54608/76‑1 (MQ=255)             cAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGgc                                                            <  1:3807011/76‑1 (MQ=255)              ggTGGCGGGGGTCCAGCAAGCGAGTTCTCCGGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGgcg                                                           <  2:808680/75‑1 (MQ=255)              cgTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGgcg                                                           <  1:2124998/75‑1 (MQ=255)              aGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGgcg                                                           <  2:3239872/76‑1 (MQ=255)               gTGGCTGTGTTCCAGCAACACAGTTCTCCTGGCACCCCGTATCGCGCGCGGTGGGGAATG‑TTAA‑AA‑AC‑CAGGgcgc                                                          >  2:3544488/1‑76 (MQ=255)               gTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGgcgc                                                          <  1:171011/76‑1 (MQ=255)               gTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGgcgc                                                          <  2:1739908/76‑1 (MQ=255)                 ggCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCgg                                                        >  2:3916096/1‑76 (MQ=255)                 ggCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCgg                                                        >  2:4727491/1‑76 (MQ=255)                 ggCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCgg                                                        <  1:1110328/76‑1 (MQ=255)                  gCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGa                                                       >  2:2063715/1‑76 (MQ=255)                  gCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGa                                                       <  1:3722625/76‑1 (MQ=255)                  gCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGa                                                       <  2:1175710/76‑1 (MQ=255)                  gCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGa                                                       >  1:1370563/1‑76 (MQ=255)                   cTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAg                                                      >  2:361245/1‑76 (MQ=255)                    tgtgTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGt                                                     >  1:4108793/1‑76 (MQ=255)                     gtgtTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGtt                                                    <  2:2107166/76‑1 (MQ=255)                     gtgtTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGtt                                                    <  2:4830044/76‑1 (MQ=255)                     gtgtTCCAGCAAACCAGTGCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGtt                                                    <  1:2294057/76‑1 (MQ=255)                         tCCAGCAAACCAGTTCTCCTGGCACCCGGCATCGCGCGCGGTGGGGGATG‑TT‑A‑AA‑AAACAGGGCGCGGAGTTAAtt                                                >  2:1686531/1‑76 (MQ=255)                          ccAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTc                                               <  2:1395032/76‑1 (MQ=255)                          ccAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTc                                               <  2:2129222/76‑1 (MQ=255)                             gCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAAc                                            <  1:932117/76‑1 (MQ=255)                                 aCCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTT‑A‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGt                                        >  1:3649944/1‑76 (MQ=255)                                  ccAGTTCTCCGGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGtt                                       <  2:2146723/76‑1 (MQ=255)                                   cAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGttt                                      >  1:479758/1‑76 (MQ=255)                                     gTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGa                                    >  1:3492774/1‑76 (MQ=255)                                     gTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGa                                    <  2:981404/76‑1 (MQ=255)                                      ttCTCCTGGCACCCGGTATCGCGCGCGGTGGTTACTT‑TT‑ACAA‑AC‑CAGGGCGCGGAGTTAATTCAACCAGTTTgct                                   >  2:3344475/1‑74 (MQ=255)                                      ttCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGAt                                   >  1:2588333/1‑76 (MQ=255)                                      ttCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGAt                                   >  2:3312989/1‑76 (MQ=255)                                       tctcCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATa                                  <  1:3375457/76‑1 (MQ=255)                                       tctcCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATa                                  >  1:874914/1‑76 (MQ=255)                                         tcCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAg                                >  1:126758/1‑76 (MQ=255)                                          ccTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGc                               >  1:4469487/1‑76 (MQ=255)                                          ccTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGc                               <  1:4102263/76‑1 (MQ=255)                                            tGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGGTAATTCAACCTGTTTGATAAGCCt                             >  2:491618/1‑76 (MQ=255)                                             ggCACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCACCCTGTTTGATACGCCTg                            >  2:2994/1‑76 (MQ=255)                                               cACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGga                          <  1:3973267/76‑1 (MQ=255)                                               cACCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGga                          <  1:1873002/76‑1 (MQ=255)                                                aCCCGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGgag                         >  1:166128/1‑76 (MQ=255)                                                 cccGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGgaga                        <  1:3961133/76‑1 (MQ=255)                                                  ccGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGat                       <  1:2483324/76‑1 (MQ=255)                                                   cGGTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGatt                      <  1:4653372/76‑1 (MQ=255)                                                    ggTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTACTTCAACCTGTTTGATAAGACTGGAGatta                     >  1:3428390/1‑76 (MQ=255)                                                    ggTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGatta                     >  1:4431215/1‑76 (MQ=255)                                                    ggTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGatta                     >  1:1754349/1‑76 (MQ=255)                                                    ggTATCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGatta                     >  1:4299212/1‑76 (MQ=255)                                                    ggTATCGCGCGCGGTGGGTAATG‑TT‑A‑AACAC‑CCGTGCTCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGatta                     >  1:1780831/1‑76 (MQ=255)                                                     gTATCGCGCGCGGTGGGTCATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTACTTCAACCTGTTTTATAAGCCTGGAGattat                    <  2:2691760/76‑1 (MQ=255)                                                        tCGCGCGCGGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGa                 >  1:2690745/1‑76 (MQ=255)                                                        tCGCGCGCGGGGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGa                 <  2:3379304/76‑1 (MQ=255)                                                         cgcgcgcgGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGCTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGAt                >  1:1491872/1‑76 (MQ=255)                                                           ngcgcgGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATca              >  1:1934401/2‑76 (MQ=255)                                                           cgcgcgGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATca              <  2:3511939/76‑1 (MQ=255)                                                           cgcgcgGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATca              <  1:3660661/76‑1 (MQ=255)                                                             cgcgGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATcaac            <  2:4107356/76‑1 (MQ=255)                                                              gcgGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTCTTGATcacca           >  2:92454/1‑76 (MQ=255)                                                              gcgGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATcaaca           >  2:2884584/1‑76 (MQ=255)                                                              gcgGTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATcaaca           >  1:4388211/1‑76 (MQ=255)                                                                ggTGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAAt         >  2:3865449/1‑76 (MQ=255)                                                                  tGGGTACTG‑TT‑A‑AACGC‑CAGGGCGAGGAGTTAAGTCACTCTGTTTGATAAGCCTGGAGATTATTGATcaccaaggc       <  2:2202383/76‑4 (MQ=255)                                                                  tGGGTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATAc       <  2:1813245/76‑1 (MQ=255)                                                                  ggggTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATAc       <  2:699455/75‑1 (MQ=255)                                                                   gggTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACt      >  2:4247641/1‑76 (MQ=255)                                                                    ggTAATG‑TTAA‑AA‑AC‑CAGGTCGCGGAGTTAATTCAACCTGTTTGATCAGCCTGGATATTATTGATCAACAATACTg     <  1:1686532/76‑1 (MQ=255)                                                                     gTAATG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGc    <  1:2968194/76‑1 (MQ=255)                                                                       aaTG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAGTTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGt  <  1:4451462/76‑1 (MQ=255)                                                                       aaTG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGt  <  1:2480365/76‑1 (MQ=255)                                                                       aaTG‑TTAA‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGt  >  1:77065/1‑76 (MQ=255)                                                                              |                                                                         AGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATG‑TT‑A‑AA‑AC‑CAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGT  >  NC_000913/2008401‑2008545

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑