Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,704,171	G→A	D46D (GAC→GAT)	ydgJ ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,704,171	0	G	A	98.3%	195.5 / ‑5.6	59	D46D (GAC→GAT)	ydgJ	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (30/28);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCG  >  NC_000913/1704074‑1704265                                                                                                  |                                                                                               ggAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATc                                                                                              <  2:1413246/100‑1 (MQ=255)       gggTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCttt                                                                                        <  2:713420/100‑1 (MQ=255)          tATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTAc                                                                                     <  1:1411864/100‑1 (MQ=255)             cGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACttt                                                                                  <  2:144382/100‑1 (MQ=255)                 gggTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGtt                                                                              <  2:803388/100‑1 (MQ=255)                 gggTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGtt                                                                              <  1:794013/100‑1 (MQ=255)                  ggTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGttt                                                                             <  1:1175059/100‑1 (MQ=255)                   gtgtAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGttt                                                                             <  2:761444/99‑1 (MQ=255)                    tgtAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGACGCTTCAGCCCAGAGAAAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTtca                                                                           <  2:535489/100‑1 (MQ=255)                     gtAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTtcat                                                                          >  1:803815/1‑100 (MQ=255)                      tAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTtcatc                                                                         >  1:1748071/1‑100 (MQ=255)                      tAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTtcatc                                                                         <  1:900275/100‑1 (MQ=255)                      tAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTtcatc                                                                         >  2:363823/1‑100 (MQ=255)                        ggAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCAc                                                                       <  1:448675/100‑1 (MQ=255)                          aaTGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCActg                                                                     >  2:991322/1‑100 (MQ=255)                                 aTCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGgaga                                                              <  2:1777714/100‑1 (MQ=255)                                  tCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGAt                                                             >  1:185944/1‑100 (MQ=255)                                   cAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCAATGCTAGAGAtt                                                            >  2:1495778/1‑100 (MQ=255)                                      gTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGTTTTTACTTTTGTTTCATCACTGCTGGAGATTAcc                                                         <  1:6321/100‑1 (MQ=255)                                       tCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCg                                                        >  1:1482481/1‑100 (MQ=255)                                            gTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGt                                                   >  2:1532678/1‑100 (MQ=255)                                                ggATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCt                                               >  2:109768/1‑100 (MQ=255)                                                   tCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTTTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGcc                                            >  2:308834/1‑100 (MQ=255)                                                     gTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCg                                          >  2:222883/1‑100 (MQ=255)                                                      ttccgatcttGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCgg                                         <  2:555508/91‑1 (MQ=255)                                                       tAAACAGATGCTTGGGCTCAGAGACAACCGTAACCGTTAGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCggg                                        >  1:1121043/1‑100 (MQ=255)                                                         aaCAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCg                                      >  2:597848/1‑100 (MQ=255)                                                           cAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTg                                    <  2:607012/100‑1 (MQ=255)                                                             gATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCTAATCTGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGcc                                  <  2:180470/100‑1 (MQ=255)                                                               tGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCggagatcggaa                                >  1:555504/1‑91 (MQ=255)                                                                gCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGc                               <  1:801235/100‑1 (MQ=255)                                                                gCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGc                               <  1:403851/100‑1 (MQ=255)                                                                  ttCGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCaa                             >  1:1107608/1‑100 (MQ=255)                                                                    cGGCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCACGCCCgggcag                           <  2:102325/100‑6 (MQ=255)                                                                    cGGCTCAGAGACAACCGTAACCGTTGGCAAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATc                           <  1:1796904/100‑1 (MQ=255)                                                                     ggCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCa                          >  1:1258479/1‑100 (MQ=255)                                                                     ggCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCa                          >  1:1536250/1‑100 (MQ=255)                                                                     ggCTCAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCa                          >  1:1305901/1‑100 (MQ=255)                                                                         cAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACTGCCAGTTCCTGCACGAACGAACCCGCAATCAgg                        >  2:756680/1‑98 (MQ=255)                                                                         cAGAGACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAgggg                      <  1:268513/100‑1 (MQ=255)                                                                          agagACAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGgc                     >  2:1772930/1‑100 (MQ=255)                                                                           gagaCAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGgcg                    >  2:106117/1‑100 (MQ=255)                                                                            agaCAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGgcgc                   >  2:1526040/1‑100 (MQ=255)                                                                             gaCAACCGTAACCGTTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCa                  <  1:1683069/100‑1 (MQ=255)                                                                                aaccgtaaccgtTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATgg               <  1:991313/100‑1 (MQ=255)                                                                                 accgtaaccgtTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGTGCGCATgta              >  2:1346764/1‑98 (MQ=255)                                                                                 accgtaaccgtTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGa              >  1:1224905/1‑100 (MQ=255)                                                                                  ccgtaaccgtTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGaa             <  2:185944/100‑1 (MQ=255)                                                                                  ccgtaaccgtTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGaa             >  2:1438862/1‑100 (MQ=255)                                                                                    ttccgatctGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAgg           <  1:166992/92‑1 (MQ=255)                                                                                    gtaaccgtTGGCCAATAGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAgg           >  1:1615207/1‑100 (MQ=255)                                                                                      aactgtTGGCCAGTCGGCTTTTACTTTAGTTTGATAGCTGCTGGAAATTACCGCCAGTTCCTGCTCGGGCTTGCCAGCAATGAGGGGAGCATGGAAgg           >  1:639459/1‑98 (MQ=255)                                                                                      aaccgtTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGt          <  1:4752/99‑1 (MQ=255)                                                                                      aaccgtTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGAGCGTGCCCGCAATCAGGGGCGCATGAAAgg           >  2:626293/1‑98 (MQ=255)                                                                                        ccgtTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGttt        <  2:878918/99‑1 (MQ=255)                                                                                        ccgtTGGCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCACGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGtttt       >  1:1677949/1‑100 (MQ=255)                                                                                             ggCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAggagatcggaa  >  2:166992/1‑91 (MQ=255)                                                                                             ggCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCt    >  1:227747/1‑98 (MQ=255)                                                                                             ggCCAATCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCg  <  2:1355230/100‑1 (MQ=255)                                                                                                  |                                                                                               GGAAATGGGTATCGTTGGGTGTAGGAATGACAATCAGGTCTATGTTGGGATCGTTAAACAGATGCTTCGGCTCAGAGACAACCGTAACCGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCG  >  NC_000913/1704074‑1704265

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑