Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	410,479	A→G	S112S (TCA→TCG)	mak →	manno(fructo)kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	410,479	0	A	G	98.6%	237.2 / ‑4.7	70	S112S (TCA→TCG)	mak	manno(fructo)kinase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (38/31);  total (39/31)
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.22e-01

CAACCTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCAGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGG  >  NC_000913/410382‑410573                                                                                                  |                                                                                               cAACCTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGa                                                                                              <  1:47326/100‑1 (MQ=255)   aCCTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGaa                                                                                             <  1:842811/99‑1 (MQ=255)   aCCTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCACGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAg                                                                                            >  1:330036/1‑100 (MQ=255)     cTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCa                                                                                          <  2:1321393/100‑1 (MQ=255)      tGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGaagaag                                                                                         >  1:319629/1‑100 (MQ=255)      tGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGGAGCAg                                                                                         >  1:1206153/1‑100 (MQ=255)      tGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAg                                                                                         <  1:396585/100‑1 (MQ=255)        gCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTa                                                                                       >  2:506103/1‑100 (MQ=255)        gCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTa                                                                                       <  1:402664/100‑1 (MQ=255)          tCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGa                                                                                     >  2:1755815/1‑100 (MQ=255)           cAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGAt                                                                                    >  1:472213/1‑100 (MQ=255)           cAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGAt                                                                                    >  1:1534682/1‑100 (MQ=255)              cGGTCAGCCATTCGATAAAGTCTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGgc                                                                                 <  2:31502/100‑1 (MQ=255)              cGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGgc                                                                                 >  1:298799/1‑100 (MQ=255)                 tCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGcgcg                                                                              <  2:1568000/100‑1 (MQ=255)                   aGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAATCAGTAGATGGCGCGGc                                                                            <  1:813960/100‑1 (MQ=255)                   aGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGgcgc                                                                               >  1:212543/1‑97 (MQ=255)                    gCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCa                                                                           >  2:1002832/1‑100 (MQ=255)                       aTTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCg                                                                        >  1:1807246/1‑100 (MQ=255)                       aTTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCg                                                                        >  1:1441976/1‑100 (MQ=255)                        ttCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCAGAAGCAGTAGATGGCGCGGCAGCgg                                                                       >  1:79975/1‑100 (MQ=255)                         tCGATAAAGACTTAAGCGCGATGTTGCAGCGGGAAGTGCGGTTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCggg                                                                      <  2:548799/100‑1 (MQ=255)                           gATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAg                                                                    <  1:1181746/100‑1 (MQ=255)                           gATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCCGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCTGCGGGAg                                                                    >  1:972845/1‑100 (MQ=255)                            aTAAAGACTTAAGCGCCAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAgc                                                                   <  1:1596216/100‑1 (MQ=255)                             tAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAgcg                                                                  >  1:598344/1‑100 (MQ=255)                              aaaGACTTAAGCGCGAGGTTGCAGTGGGAAGTGCTGCTGGACAATGCCGCTAACTATCTGGCGGTTTCGGAAGCTGTAGATTGCGCGGCAGCGGGAgc                                                                   <  1:364023/98‑1 (MQ=255)                              aaaGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGTAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAgcg                                                                  <  1:539933/99‑1 (MQ=255)                              aaaGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAgcgc                                                                 >  2:1769548/1‑100 (MQ=255)                              aaaGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAgcgc                                                                 <  2:689336/100‑1 (MQ=255)                              aaaGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAgcgc                                                                 <  1:1095075/100‑1 (MQ=255)                                aGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAg                                                               >  1:386240/1‑100 (MQ=255)                                       aGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTAt                                                        >  2:1330549/1‑100 (MQ=255)                                         cgcgAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAACGCAGACGGTAttt                                                      >  2:309588/1‑100 (MQ=255)                                            gAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGcc                                                   <  1:16655/100‑1 (MQ=255)                                              ggTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGt                                                 >  1:862684/1‑100 (MQ=255)                                              ggTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGt                                                 <  2:1169284/100‑1 (MQ=255)                                               gTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTg                                                >  2:1477936/1‑100 (MQ=255)                                               gTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTATCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTg                                                >  1:734598/1‑100 (MQ=255)                                                 tGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGAt                                              >  1:832832/1‑100 (MQ=255)                                                    aGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTa                                            <  2:865794/99‑1 (MQ=255)                                                     gCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATc                                          >  2:1256735/1‑100 (MQ=255)                                                     gCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTAGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACAGTATTTGCCGTGATTATc                                          >  1:554066/1‑100 (MQ=255)                                                     gCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTATGGCAGTCTCGGATGTAGTAGATGGCTCAGCAGCGGGTGCGCTTATGGTATTTGCTGTCATTAtt                                          <  1:114925/100‑2 (MQ=255)                                                       gggAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCgg                                        <  2:985364/100‑1 (MQ=255)                                                        tgAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGc                                       <  2:1397120/99‑1 (MQ=255)                                                        ggAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGAGATTaacag                                        >  2:897679/1‑95 (MQ=255)                                                         gAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCa                                      >  1:1397109/1‑100 (MQ=255)                                                           aGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGc                                       >  1:1130320/1‑97 (MQ=255)                                                           aGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACg                                    >  2:1680617/1‑100 (MQ=255)                                                           aGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACg                                    <  2:862690/100‑1 (MQ=255)                                                            gTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACgg                                   >  2:310052/1‑100 (MQ=255)                                                             tGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACggg                                  <  1:1493011/100‑1 (MQ=255)                                                             tGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACggg                                  >  2:483975/1‑100 (MQ=255)                                                              gCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGa                                 >  2:1791510/1‑100 (MQ=255)                                                                 gCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATgc                              >  2:1002978/1‑100 (MQ=255)                                                                 gCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATgc                              <  2:563526/100‑1 (MQ=255)                                                                 gCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATgc                              <  2:710087/100‑1 (MQ=255)                                                                      cAAATGACGCTAACTGTCTGGCGGTATCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATgcggcg                          >  1:396199/1‑99 (MQ=255)                                                                             cGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTg                  <  2:1441987/100‑1 (MQ=255)                                                                              gCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTgg                 <  2:1080489/100‑1 (MQ=255)                                                                                     gTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCaa          >  2:1147761/1‑100 (MQ=255)                                                                                      tCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAAt         <  2:90372/100‑1 (MQ=255)                                                                                        tGGCGTTTCCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGTTTATCGGCACTGGATGCGGCGTGGGCGTGGCATTCAATgg       >  1:81294/1‑100 (MQ=255)                                                                                        tGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATgg       <  2:1141407/100‑1 (MQ=255)                                                                                           cGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATgggcg    <  1:1109887/100‑1 (MQ=255)                                                                                           cGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATgggcg    >  2:1427071/1‑100 (MQ=255)                                                                                           cGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATgggcg    <  1:1109904/100‑1 (MQ=255)                                                                                             gTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATgggcggg  >  2:844368/1‑100 (MQ=255)                                                                                             gTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATgggcggg  >  1:884687/1‑100 (MQ=255)                                                                                                  |                                                                                               CAACCTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCAGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGG  >  NC_000913/410382‑410573

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