Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	911,132	A→G	intergenic (‑83/+50)	poxB ← / ← hcr	pyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding/HCP oxidoreductase, NADH‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	911,132	0	A	G	100.0%	236.3 / NA	66	intergenic (‑83/+50)	poxB/hcr	pyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding/HCP oxidoreductase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/32);  total (34/32)

GCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTT‑GTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGT  >  NC_000913/911003‑911262                                                                                                                                   |                                                                                                                                   gCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTT‑GTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGtt                                                                                                                           >  1:211232/1‑139 (MQ=255)            tgtTTTGGCGATATAAGCTGCAACCGTTTGGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCg                                                                                                                 >  1:207426/1‑139 (MQ=255)                ttGGCGATATAAGCTGCAACCGTTT‑GTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGc                                                                                                                    <  2:79866/131‑1 (MQ=255)                ttGGCGATATAAGCTGCAACCGTTT‑GTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGc                                                                                                                    >  1:79866/1‑131 (MQ=255)                ttGGCGATATAAGCTGCAACCGTTT‑GTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                            <  2:1265052/139‑1 (MQ=255)                               gCAACCGTTT‑GTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCggg                                                                                             <  1:1089061/139‑1 (MQ=255)                                 aaCCGTTT‑GTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                                  <  2:359221/100‑1 (MQ=255)                                 aaCCGTTT‑GTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                                  >  1:359221/1‑100 (MQ=255)                                 aaCCGTTT‑GTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCa                                                                                           <  2:211232/139‑1 (MQ=255)                                            ttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGAtt                                                                                 <  2:276530/139‑1 (MQ=255)                                                    tctcCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCgaga                                                                         >  2:281012/1‑139 (MQ=255)                                                           ctcCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAAt                                                                  <  1:1181643/139‑1 (MQ=255)                                                             ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGGAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATcc                                                                >  1:628092/1‑139 (MQ=255)                                                             ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                                                 <  1:166613/106‑1 (MQ=255)                                                             ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                                                 >  2:166613/1‑106 (MQ=255)                                                             ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCg                                                                                               <  2:1341958/108‑1 (MQ=255)                                                             ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCg                                                                                               >  1:1341958/1‑108 (MQ=255)                                                                  aTGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTg                                                           <  2:167838/139‑1 (MQ=255)                                                                   tgtgATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                                                 <  2:1569633/100‑1 (MQ=255)                                                                   tgtgATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                                                 >  1:1569633/1‑100 (MQ=255)                                                                     tgATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGGAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCgg                                                        <  2:628092/139‑1 (MQ=255)                                                                        tAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                            >  1:673128/1‑84 (MQ=255)                                                                        tAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                            <  2:673128/84‑1 (MQ=255)                                                                        tAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                            >  2:962267/1‑84 (MQ=255)                                                                        tAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAAGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                            <  1:962267/84‑1 (MQ=255)                                                                           cGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGTGGTGCTAACGGTTAACTAGCCCGATGAAAGGAATATCATCGGGCa                                                                                           <  1:930354/98‑1 (MQ=255)                                                                           cGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGTGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCa                                                                                           >  2:930354/1‑98 (MQ=255)                                                                             gTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAgg                                                >  1:761847/1‑139 (MQ=255)                                                                               aaCAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCTGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAg                                              <  1:365937/139‑1 (MQ=255)                                                                               aaCAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                            <  2:773424/77‑1 (MQ=255)                                                                               aaCAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                            >  1:773424/1‑77 (MQ=255)                                                                                   aGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAgg                                          <  2:693412/139‑1 (MQ=255)                                                                                        aGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATg                                                                                                              <  1:1436258/66‑1 (MQ=255)                                                                                        aGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATg                                                                                                              >  2:1436258/1‑66 (MQ=255)                                                                                        aGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCa                                                                     <  2:1016086/107‑1 (MQ=255)                                                                                        aGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGt                                     >  1:1049489/1‑139 (MQ=255)                                                                                        aGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGAGAGAACCa                                                                     >  1:1016086/1‑107 (MQ=255)                                                                                          ttCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACg                                                                                                                              <  1:854390/48‑1 (MQ=255)                                                                                          ttCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACg                                                                                                                              >  2:854390/1‑48 (MQ=255)                                                                                            cATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCgag                                                                          <  2:1530297/98‑1 (MQ=255)                                                                                            cATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCgag                                                                          >  1:1530297/1‑98 (MQ=255)                                                                                            cATCTGAAGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGACAGGAGCAGGCCAGTACGt                                 >  1:1053388/1‑139 (MQ=255)                                                                                              tCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTaa                               <  1:210001/139‑1 (MQ=255)                                                                                               cTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                                                 <  1:201357/72‑1 (MQ=255)                                                                                               cTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                                                 >  2:201357/1‑72 (MQ=255)                                                                                               cTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAAc                              >  2:244584/1‑139 (MQ=255)                                                                                                      aGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATACCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGc                       >  2:1433884/1‑139 (MQ=255)                                                                                                        ggggAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGa                     >  1:605993/1‑139 (MQ=255)                                                                                                                gATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCg             >  1:846502/1‑139 (MQ=255)                                                                                                                gATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGCACGTAACCTTCAGCGATTTCGGCg             <  1:573381/139‑1 (MQ=255)                                                                                                                  tGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                            <  1:190243/42‑1 (MQ=255)                                                                                                                  tGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                            >  2:190243/1‑42 (MQ=255)                                                                                                                   gggAGAGAAAGGAGTCGCTAACGGTTAAATAGCCCGATGaa                                                                                                            >  2:979121/1‑41 (MQ=255)                                                                                                                   gggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                            <  1:979121/41‑1 (MQ=255)                                                                                                                   gggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATa                                                                                         >  2:613959/1‑60 (MQ=255)                                                                                                                   gggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATa                                                                                         <  1:613959/60‑1 (MQ=255)                                                                                                                    ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTg                                                           >  2:408994/1‑89 (MQ=255)                                                                                                                    ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTg                                                           <  1:408994/89‑1 (MQ=255)                                                                                                                         gaAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGt        >  2:540466/1‑135 (MQ=255)                                                                                                                         gaAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGt        <  1:540466/135‑1 (MQ=255)                                                                                                                          aaaGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCgag                                                                          <  2:584300/68‑1 (MQ=255)                                                                                                                          aaaGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCgag                                                                          >  1:584300/1‑68 (MQ=255)                                                                                                                          aaaGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGc                                                   <  2:1101902/91‑1 (MQ=255)                                                                                                                          aaaGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGc                                                   >  1:1101902/1‑91 (MQ=255)                                                                                                                          aaaGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCg   >  1:1018804/1‑139 (MQ=255)                                                                                                                           aaGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGt  >  2:545235/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                   GCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTT‑GTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGT  >  NC_000913/911003‑911262

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