Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,643,411	G→A	R83C (CGC→TGC)	relE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,643,411	0	G	A	97.7%	297.0 / ‑3.6	86	R83C (CGC→TGC)	relE	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (34/50);  total (35/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

AAGCAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGTTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAA  >  NC_000913/1643322‑1643500                                                                                          |                                                                                          aaGCAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATtc                                                                                        <  2:392415/93‑1 (MQ=255) aaGCAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATtc                                                                                        <  1:486386/93‑1 (MQ=255) aaGCAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATtc                                                                                        <  1:2312925/93‑1 (MQ=255)  aGCAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATtct                                                                                       >  1:1569368/1‑93 (MQ=255)   gCAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCAATATACTTCCGAGCATtctc                                                                                      <  2:667261/93‑1 (MQ=255)     aGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCtt                                                                                    <  2:980664/93‑1 (MQ=255)      gAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCttt                                                                                   >  2:2178932/1‑93 (MQ=255)           gtCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCa                                                                              >  1:3769726/1‑93 (MQ=255)             cACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAAc                                                                            <  2:1101423/93‑1 (MQ=255)             cACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAAc                                                                            >  1:3509746/1‑93 (MQ=255)             cACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCATTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAAc                                                                            <  2:2520987/93‑1 (MQ=255)              aCCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACa                                                                           <  2:3086595/93‑1 (MQ=255)                cTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGa                                                                         >  2:3668971/1‑93 (MQ=255)                 ttCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGaa                                                                        >  2:2182816/1‑93 (MQ=255)                   cGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAAt                                                                      <  1:633687/93‑1 (MQ=255)                    ggTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGAACGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATc                                                                     >  2:735109/1‑93 (MQ=255)                    cgTGCCAAACTGAGATTTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATc                                                                     <  2:1433001/92‑1 (MQ=255)                         gAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTATCTTTTCCCAACAGAAATCACGaa                                                                >  1:3030481/1‑93 (MQ=255)                           aaCAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGaaaa                                                              <  1:4419968/93‑1 (MQ=255)                            aCAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAAc                                                             <  2:1365667/93‑1 (MQ=255)                            aCAGAGATGTCATGCTGTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAAc                                                             <  2:2576872/93‑1 (MQ=255)                             cAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACg                                                            >  1:1781287/1‑93 (MQ=255)                               gagaTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGAc                                                          >  2:642609/1‑93 (MQ=255)                                 gaTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACaa                                                        <  2:852119/93‑1 (MQ=255)                                 gaTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACaa                                                        <  2:1983113/93‑1 (MQ=255)                                  aTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGTTCTCTTTTCCCAACAGAAATCACGAAAACGACAAc                                                       <  2:825985/93‑1 (MQ=255)                                  aTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAAc                                                       >  2:273975/1‑93 (MQ=255)                                  aTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAAc                                                       <  2:832240/93‑1 (MQ=255)                                  aTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAAc                                                       >  2:1417454/1‑93 (MQ=255)                                  aTGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAAc                                                       <  1:1386644/93‑1 (MQ=255)                                    gTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATCCTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACtt                                                     <  2:4324175/93‑1 (MQ=255)                                     tCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACttt                                                    <  2:1637706/93‑1 (MQ=255)                                     tCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACttt                                                    <  1:1384709/93‑1 (MQ=255)                                     tCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACttt                                                    >  2:328667/1‑93 (MQ=255)                                     tCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACttt                                                    >  1:4290038/1‑93 (MQ=255)                                        tGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTtctc                                                 <  1:961597/93‑1 (MQ=255)                                         gCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCg                                                <  1:1227518/93‑1 (MQ=255)                                             tGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCta                                            >  1:584879/1‑93 (MQ=255)                                             tGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCta                                            <  1:1716094/93‑1 (MQ=255)                                              ggTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCtat                                           <  2:2991989/93‑1 (MQ=255)                                                   agagAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCt                                      <  2:4242230/93‑1 (MQ=255)                                                   agagAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCt                                      <  2:2066192/93‑1 (MQ=255)                                                    gagaATGCGTTTGACCGCCTCTCTATATACTGCCGAGCATGCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTg                                     <  1:2933362/93‑1 (MQ=255)                                                    gagaATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCACTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTg                                     >  1:4660833/1‑93 (MQ=255)                                                     agaATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTgg                                    <  2:3655687/93‑1 (MQ=255)                                                      gaATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGt                                   >  2:3880529/1‑93 (MQ=255)                                                      gaATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGt                                   >  1:4116517/1‑93 (MQ=255)                                                      gaATGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGt                                   <  1:1967223/93‑1 (MQ=255)                                                       aaTGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGta                                  >  2:3896687/1‑93 (MQ=255)                                                       aaTGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGta                                  <  1:3443345/93‑1 (MQ=255)                                                       aaTGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGta                                  <  1:3929355/93‑1 (MQ=255)                                                        aTGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGtat                                 <  1:1349125/93‑1 (MQ=255)                                                         tGCGTTTGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGtata                                <  1:2106848/93‑1 (MQ=255)                                                          gCGTTTGACCGCCTCGCTATATACTTCCGAGCGTTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATAc                               >  1:1685752/1‑93 (MQ=255)                                                             tttGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAg                            >  2:1294816/1‑93 (MQ=255)                                                             tttGACCGCCTCGCTATATACTACCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACCACAACTTTCTCGCCTATAACCTGGTATACAAg                            <  2:3508689/93‑1 (MQ=255)                                                               tGACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGc                          >  1:3902187/1‑93 (MQ=255)                                                                gACCGCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCg                         <  1:3232116/93‑1 (MQ=255)                                                                    gCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAg                     >  2:4227718/1‑93 (MQ=255)                                                                    gCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAg                     <  1:1519175/93‑1 (MQ=255)                                                                    gCCTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAg                     >  1:3351899/1‑93 (MQ=255)                                                                     ccTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGc                    >  2:1776617/1‑93 (MQ=255)                                                                      cTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGcc                   <  2:3496557/93‑1 (MQ=255)                                                                      cTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGcc                   <  1:2405835/93‑1 (MQ=255)                                                                      cTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGcc                   >  1:3310784/1‑93 (MQ=255)                                                                      cTCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGcc                   <  1:3384606/93‑1 (MQ=255)                                                                       tCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCt                  <  1:2859119/93‑1 (MQ=255)                                                                       tCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCt                  >  2:3581188/1‑93 (MQ=255)                                                                       tCGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCt                  <  1:102594/93‑1 (MQ=255)                                                                        cGCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCTg                 <  2:1953885/93‑1 (MQ=255)                                                                         gCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCTGa                <  1:3737437/93‑1 (MQ=255)                                                                         gCTATATACTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCTGa                >  2:86400/1‑93 (MQ=255)                                                                           tatataCTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCTGAAg              >  2:18386/1‑93 (MQ=255)                                                                            atataCTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGa             <  1:3848363/93‑1 (MQ=255)                                                                            atataCTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCTGAAGa             <  1:4212201/93‑1 (MQ=255)                                                                             tataCTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCTGAAGAc            >  2:1156931/1‑93 (MQ=255)                                                                             tataCTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCTGAAGAc            <  1:136367/93‑1 (MQ=255)                                                                               taCTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCg          <  2:4308115/93‑1 (MQ=255)                                                                               taCTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCg          <  2:1444853/93‑1 (MQ=255)                                                                                aCTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCgg         >  1:3095254/1‑93 (MQ=255)                                                                                aCTTCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCgg         <  1:4004909/93‑1 (MQ=255)                                                                                  ttCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAg       >  2:3013187/1‑93 (MQ=255)                                                                                  ttCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAg       >  1:431093/1‑93 (MQ=255)                                                                                   tCCGAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGc      <  1:3952764/93‑1 (MQ=255)                                                                                      gAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTa   >  1:2800822/1‑93 (MQ=255)                                                                                      gAGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTa   >  1:4166680/1‑93 (MQ=255)                                                                                       aGCATTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGCCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTaa  <  1:897187/93‑1 (MQ=255)                                                                                          |                                                                                          AAGCAGAAGTGTCACCTTCGGTGCGAAACAGAGATGTCATGCTTTGGTTCAGAGAATGCGTTTGACCGCCTCGCTATATACTTCCGAGCGTTCTCTTTTCCCAACAGAAATCACGAAAACGACAACTTTCTCGTCTATAACCTGGTATACAAGGCGATAGCCTGAAGACCGGAGCTTAA  >  NC_000913/1643322‑1643500

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