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%	131.0 / NA	37	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 (19/18);  total (19/18)

CACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAG  >  NC_000913/910996‑911259                                                                                                                                         |                                                                                                                                cACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                               <  1:440113/139‑1 (MQ=255)               cGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATCGGAGAGAAAGGAGGCGCTAACGGTTAAATAGc                                                                                                                 <  1:599116/139‑1 (MQ=255)                        tGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaaa                                                                                                        <  1:592769/139‑1 (MQ=255)                               ataAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                         <  2:71631/131‑1 (MQ=255)                               ataAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                         >  1:71631/1‑131 (MQ=255)                                   gCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcatc                                                                                             >  2:933360/1‑139 (MQ=255)                                           cGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATaa                                                                                     <  1:933360/139‑1 (MQ=255)                                            gtttgtttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCa                                                                                        <  2:118426/135‑1 (MQ=255)                                            gtttgtttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCa                                                                                        >  1:118426/1‑135 (MQ=255)                                                           ctcCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                         <  1:163759/103‑1 (MQ=255)                                                           ctcCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                         >  2:163759/1‑103 (MQ=255)                                                            tcCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGAGGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAAc                                                                    >  1:632269/1‑139 (MQ=255)                                                                   ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAata                                                                                                  <  2:451166/102‑1 (MQ=255)                                                                   ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAata                                                                                                  >  1:451166/1‑102 (MQ=255)                                                                   ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCga                                                                        <  2:347033/128‑1 (MQ=255)                                                                   ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCga                                                                        >  1:347033/1‑128 (MQ=255)                                                                   ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATcc                                                             >  1:611782/1‑139 (MQ=255)                                                                   ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATcc                                                             >  2:337237/1‑139 (MQ=255)                                                                                 cGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCa                                               >  2:268882/1‑139 (MQ=255)                                                                                  ggTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAg                                              <  2:611782/139‑1 (MQ=255)                                                                                     aaCAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCGCCTGCGGATGGCAGGAg                                           <  2:632269/139‑1 (MQ=255)                                                                                      aCAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAg                                           >  1:363840/1‑138 (MQ=255)                                                                                      aCAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAg                                           <  2:363840/138‑1 (MQ=255)                                                                                            ttAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCa                                    <  1:268882/139‑1 (MQ=255)                                                                                              aGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGc                                                <  1:797740/125‑1 (MQ=255)                                                                                              aGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGc                                                >  2:797740/1‑125 (MQ=255)                                                                                                                        tGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTc        >  2:581039/1‑139 (MQ=255)                                                                                                                        tGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTc        >  2:853776/1‑139 (MQ=255)                                                                                                                         gggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCg       >  2:760360/1‑139 (MQ=255)                                                                                                                          ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCgg      >  2:3490/1‑139 (MQ=255)                                                                                                                          ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCgg      >  2:236589/1‑139 (MQ=255)                                                                                                                           gagagaAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGt     <  1:853776/139‑1 (MQ=255)                                                                                                                              agaAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAg  <  1:581039/139‑1 (MQ=255)                                                                                                                                aaaGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGAtt                                                                              <  1:734418/61‑1 (MQ=255)                                                                                                                                aaaGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGAtt                                                                              >  2:734418/1‑61 (MQ=255)                                                                                                                                    gaggCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACg                               <  2:55083/104‑1 (MQ=255)                                                                                                                                    gaggCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACg                               >  1:55083/1‑104 (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 < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑