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%	194.8 / NA	55	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 (26/29);  total (26/29)

ATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGA  >  NC_000913/911007‑911243                                                                                                                              |                                                                                                                aTTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAAt                                                                                                    >  1:264271/1‑139 (MQ=255)     gAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGcc                                                                                                <  1:295796/139‑1 (MQ=255)            ttGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                         <  2:264271/139‑1 (MQ=255)             tGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaaa                                                                                        >  1:1023373/1‑139 (MQ=255)                cGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGa                                                                                     <  2:1023373/139‑1 (MQ=255)                    ataAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTaa                                                                                                             <  1:741986/111‑1 (MQ=255)                    ataAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTaa                                                                                                             >  2:741986/1‑111 (MQ=255)                    ataAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAatat                                                                                 >  1:1114075/1‑139 (MQ=255)                       aGCTGCAACCGTTTATTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGATCATCTGACGGAGGGGGAAGGGATCGGAGAGAAAGGAGGCGCTAACGGTTATATAGCCCGATGAAAGGGATatcat                                                                              <  1:922216/139‑1 (MQ=255)                           gCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCggg                                                                          <  2:1114075/139‑1 (MQ=255)                             aaCCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCa                                                                        >  1:433690/1‑139 (MQ=255)                             aaCCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCa                                                                        >  1:1238846/1‑139 (MQ=255)                               ccGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATa                                                                      <  2:433690/139‑1 (MQ=255)                                                    aTCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCaa                                                 >  1:226429/1‑139 (MQ=255)                                                     tctcCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCaaa                                                <  2:1238846/139‑1 (MQ=255)                                                        ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                               >  2:995811/1‑73 (MQ=255)                                                        ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                               <  1:995811/73‑1 (MQ=255)                                                                 gATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGa                                    >  1:22899/1‑139 (MQ=255)                                                                  aTAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGc                                                                  >  1:1211171/1‑108 (MQ=255)                                                                  aTAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGc                                                                  <  2:1211171/108‑1 (MQ=255)                                                                        gTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCTGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCg                                      <  1:900406/130‑1 (MQ=255)                                                                        gTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCTGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCg                                      >  2:900406/1‑130 (MQ=255)                                                                           aCAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                              >  2:658084/1‑87 (MQ=255)                                                                           aCAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                              <  1:658084/87‑1 (MQ=255)                                                                               gTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGc                      <  2:226429/139‑1 (MQ=255)                                                                                tttAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGcc                     <  2:1058320/139‑1 (MQ=255)                                                                                 ttAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCa                    <  1:311044/139‑1 (MQ=255)                                                                                  tAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                               <  2:1474496/47‑1 (MQ=255)                                                                                  tAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                               >  1:1474496/1‑47 (MQ=255)                                                                                  tAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAg                           <  1:407040/131‑1 (MQ=255)                                                                                  tAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAg                           >  2:407040/1‑131 (MQ=255)                                                                                  tAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAg                   >  1:320526/1‑139 (MQ=255)                                                                                   aGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                               >  2:160827/1‑46 (MQ=255)                                                                                   aGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                               <  1:799991/46‑1 (MQ=255)                                                                                   aGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                               <  1:547071/46‑1 (MQ=255)                                                                                   aGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                               <  1:160827/46‑1 (MQ=255)                                                                                   aGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTaaa                                                                                                     <  2:1262908/56‑1 (MQ=255)                                                                                   aGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTaaa                                                                                                     >  1:1262908/1‑56 (MQ=255)                                                                                     ttCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATcccc                                           <  2:1441770/112‑1 (MQ=255)                                                                                     ttCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATcccc                                           >  1:1441770/1‑112 (MQ=255)                                                                                     ttCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAg                           >  1:576747/1‑128 (MQ=255)                                                                                     ttCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAg                           <  2:576747/128‑1 (MQ=255)                                                                                       cATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCg                                                                                              <  1:331564/59‑1 (MQ=255)                                                                                       cATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCg                                                                                              >  2:331564/1‑59 (MQ=255)                                                                                         tCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                         <  1:617934/62‑1 (MQ=255)                                                                                         tCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                         >  2:617934/1‑62 (MQ=255)                                                                                         tCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCa                                                                        <  1:1099278/79‑1 (MQ=255)                                                                                         tCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCa                                                                        >  2:1099278/1‑79 (MQ=255)                                                                                          cTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                               >  2:449637/1‑39 (MQ=255)                                                                                          cTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                               <  1:449637/39‑1 (MQ=255)                                                                                                   ggggAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGa  >  1:1366633/1‑139 (MQ=255)                                                                                                                   agaAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                              <  1:231983/47‑1 (MQ=255)                                                                                                                   agaAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                              >  2:231983/1‑47 (MQ=255)                                                                                                                          aggCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCAt                                                                       >  2:642559/1‑47 (MQ=255)                                                                                                                          aggCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCAt                                                                       <  1:642559/47‑1 (MQ=255)                                                                                                                              |                                                                                                                ATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGA  >  NC_000913/911007‑911243

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: ‑