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%	201.4 / NA	57	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 (28/29);  total (28/29)

CACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAG  >  NC_000913/910996‑911259                                                                                                                                         |                                                                                                                                cACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                               <  1:42583/139‑1 (MQ=255)     ccTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACg                                                                                                                           <  1:497939/139‑1 (MQ=255)               cGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAA‑GAATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGcc                                                                                                                <  1:1048432/139‑1 (MQ=255)                   tgtTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGa                                                                                                             >  1:605447/1‑139 (MQ=255)                         ggCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAg                                                                                                       <  2:605447/139‑1 (MQ=255)                             atatAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAat                                                                                                   <  1:103197/139‑1 (MQ=255)                                taAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGc                                                                                                                 <  1:659552/122‑1 (MQ=255)                                taAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGAGCTAACGGTTAAATAGc                                                                                                                 >  2:659552/1‑122 (MQ=255)                                             tttgtttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAAc                                                                                                                            <  1:133153/98‑1 (MQ=255)                                             tttgtttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAAc                                                                                                                            >  2:133153/1‑98 (MQ=255)                                             tttgtttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGcc                                                                                                                >  2:1140964/1‑110 (MQ=255)                                             tttgtttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGcc                                                                                                                <  1:1140964/110‑1 (MQ=255)                                                                   ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                               <  2:1319503/73‑1 (MQ=255)                                                                   ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTa                                                                             <  1:850331/123‑1 (MQ=255)                                                                   ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTa                                                                             >  2:850331/1‑123 (MQ=255)                                                                       aaTGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGAt                                                                                                            >  1:822582/1‑88 (MQ=255)                                                                       aaTGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGAt                                                                                                            <  2:822582/88‑1 (MQ=255)                                                                       aaTGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaaa                                                                                                        <  1:1258348/92‑1 (MQ=255)                                                                       aaTGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaaa                                                                                                        >  2:1258348/1‑92 (MQ=255)                                                                         tgtgATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGc                                                       >  1:1123826/1‑139 (MQ=255)                                                                            gATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                         <  1:1444572/86‑1 (MQ=255)                                                                            gATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                         >  2:1444572/1‑86 (MQ=255)                                                                            gATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGa                                                    >  1:1228968/1‑139 (MQ=255)                                                                              tAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                               <  1:252977/62‑1 (MQ=255)                                                                              tAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATg                                                  >  1:1449179/1‑139 (MQ=255)                                                                              tAACGGTAACAAGTTTAGTTCATCTGACGAAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                               >  2:252977/1‑62 (MQ=255)                                                                                 cGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCa                                               >  1:1410231/1‑139 (MQ=255)                                                                                   gTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCa                                                                                        <  2:1352292/96‑1 (MQ=255)                                                                                   gTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCa                                                                                        >  1:1352292/1‑96 (MQ=255)                                                                                     aaCAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAg                                           <  2:1449179/139‑1 (MQ=255)                                                                                       cAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCa                                         <  1:1203400/139‑1 (MQ=255)                                                                                            ttAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                               <  1:1084491/48‑1 (MQ=255)                                                                                            ttAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                               >  2:1084491/1‑48 (MQ=255)                                                                                              aGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGt                                  >  1:661147/1‑139 (MQ=255)                                                                                              aGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGt                                  >  1:648930/1‑139 (MQ=255)                                                                                                        aCGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCggg                                                                                          >  2:1300907/1‑73 (MQ=255)                                                                                                        aCGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCggg                                                                                          <  1:1300907/73‑1 (MQ=255)                                                                                                          ggAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGcc                                                                                                                >  2:1165760/1‑49 (MQ=255)                                                                                                          ggAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGcc                                                                                                                <  1:1165760/49‑1 (MQ=255)                                                                                                               gggAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAg                                                                                    >  1:51496/1‑72 (MQ=255)                                                                                                               gggAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAg                                                                                    <  2:51496/72‑1 (MQ=255)                                                                                                                 gAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGAttt               <  2:1410231/139‑1 (MQ=255)                                                                                                                        tGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTc        >  1:1136777/1‑139 (MQ=255)                                                                                                                          ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAgg                                                                                   <  2:1055065/62‑1 (MQ=255)                                                                                                                          ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAgg                                                                                   >  1:1055065/1‑62 (MQ=255)                                                                                                                          ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTa                                                                             >  1:171327/1‑68 (MQ=255)                                                                                                                          ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTa                                                                             <  2:171327/68‑1 (MQ=255)                                                                                                                          ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCgaga                                                                      <  2:214529/75‑1 (MQ=255)                                                                                                                          ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCgaga                                                                      >  1:214529/1‑75 (MQ=255)                                                                                                                          ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGaa                                                                     >  2:1475118/1‑76 (MQ=255)                                                                                                                          ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGaa                                                                     <  1:1475118/76‑1 (MQ=255)                                                                                                                              agaAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                                              <  2:449440/47‑1 (MQ=255)                                                                                                                              agaAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                                              >  1:449440/1‑47 (MQ=255)                                                                                                                              agaAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAg  <  2:1136777/139‑1 (MQ=255)                                                                                                                              agaAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAg  >  2:1362403/1‑139 (MQ=255)                                                                                                                                aaaGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGc                                          >  1:881446/1‑97 (MQ=255)                                                                                                                                aaaGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGc                                          <  2:881446/97‑1 (MQ=255)                                                                                                                                         |                                                                                                                                CACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAG  >  NC_000913/910996‑911259

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