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%	229.3 / NA	64	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 (30/34);  total (30/34)

CGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGT  >  NC_000913/911005‑911256                                                                                                                                |                                                                                                                              cGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTaa                                                                                                                    <  2:446186/139‑1 (MQ=255)     tCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTaa                                                                                                                           <  1:1300811/128‑1 (MQ=255)     tCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTaa                                                                                                                           >  2:1300811/1‑128 (MQ=255)       gAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTGCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGcc                                                                                                              <  2:1158050/139‑1 (MQ=255)                         aGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATatcat                                                                                           <  1:1274988/139‑1 (MQ=255)                         aGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATatcat                                                                                           <  2:734090/139‑1 (MQ=255)                             gCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCggg                                                                                       <  2:1421284/139‑1 (MQ=255)                             gCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCggg                                                                                       <  2:595964/139‑1 (MQ=255)                               aaCCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATatcat                                                                                           >  2:1556078/1‑133 (MQ=255)                               aaCCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATatcat                                                                                           <  1:1556078/133‑1 (MQ=255)                                     ttgtttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGc                                                                               >  2:1328428/1‑139 (MQ=255)                                        tttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                             <  1:468644/91‑1 (MQ=255)                                        tttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                             >  2:468644/1‑91 (MQ=255)                                         ttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGAtt                                                                           <  2:1299548/139‑1 (MQ=255)                                            aTGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                       >  1:691514/1‑109 (MQ=255)                                            aTGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                       <  2:691514/109‑1 (MQ=255)                                                   tcCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAAc                                                                 >  2:106806/1‑139 (MQ=255)                                                       tctcCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCaaa                                                             >  1:1692937/1‑139 (MQ=255)                                                       tctcCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCaaa                                                             <  1:106806/139‑1 (MQ=255)                                                          ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                             <  2:814125/73‑1 (MQ=255)                                                          ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                             >  1:814125/1‑73 (MQ=255)                                                              aaTGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCt                                                      >  1:1671739/1‑139 (MQ=255)                                                                   gATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATatcat                                                                                           <  1:1597508/97‑1 (MQ=255)                                                                   gATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATatcat                                                                                           >  2:1597508/1‑97 (MQ=255)                                                                     tAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATatcat                                                                                           <  1:1725030/95‑1 (MQ=255)                                                                     tAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATatcat                                                                                           >  2:1725030/1‑95 (MQ=255)                                                                        cGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATatcat                                                                                           <  2:1224234/92‑1 (MQ=255)                                                                        cGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATatcat                                                                                           >  1:1224234/1‑92 (MQ=255)                                                                          gTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCa                                                                                     <  1:360356/96‑1 (MQ=255)                                                                          gTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCa                                                                                     >  2:360356/1‑96 (MQ=255)                                                                          gTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAgg                                          >  1:1192367/1‑139 (MQ=255)                                                                           tAACAAGTTTAGTTCATCTGACGGGGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGa                                         >  2:371753/1‑139 (MQ=255)                                                                           tAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGa                                         >  1:1023430/1‑139 (MQ=255)                                                                              cAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCa                                      <  1:410448/139‑1 (MQ=255)                                                                                 gTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGc                                   <  2:1192367/139‑1 (MQ=255)                                                                                  tttAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAA‑GAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGTAGCAGGCCa                                 <  1:1162759/139‑1 (MQ=255)                                                                                       ttCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAg                                                                                                                >  1:1530052/1‑57 (MQ=255)                                                                                       ttCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAg                                                                                                                <  2:1530052/57‑1 (MQ=255)                                                                                          aTCTGACGGGGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTa                          <  1:371753/139‑1 (MQ=255)                                                                                           tCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATatcat                                                                                           >  1:184519/1‑73 (MQ=255)                                                                                           tCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATatcat                                                                                           <  2:184519/73‑1 (MQ=255)                                                                                           tCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTaa                         >  1:691463/1‑139 (MQ=255)                                                                                            cTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGc                                       <  2:651659/124‑1 (MQ=255)                                                                                            cTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGc                                       >  1:651659/1‑124 (MQ=255)                                                                                                     ggggAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGcc                                                                                                              <  1:1573906/45‑1 (MQ=255)                                                                                                     ggggAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGcc                                                                                                              >  2:1573906/1‑45 (MQ=255)                                                                                                     ggggAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGc                                                                                                               <  1:205768/44‑1 (MQ=255)                                                                                                     ggggAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGc                                                                                                               >  2:205768/1‑44 (MQ=255)                                                                                                     ggggAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCgag                                                                    >  2:1152172/1‑86 (MQ=255)                                                                                                     ggggAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCgag                                                                    <  1:1152172/86‑1 (MQ=255)                                                                                                      gggAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGAt                                                                                                          <  1:1046564/48‑1 (MQ=255)                                                                                                      gggAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGAt                                                                                                          >  2:1046564/1‑48 (MQ=255)                                                                                                             gATGGGAGAGAAAGGAGGTGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCg       >  1:170850/1‑139 (MQ=255)                                                                                                             gATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCg                                                                                                            <  2:984253/39‑1 (MQ=255)                                                                                                             gATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCg                                                                                                            >  1:984253/1‑39 (MQ=255)                                                                                                                 ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGc                                             >  2:1017161/1‑97 (MQ=255)                                                                                                                 ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGc                                             <  1:1017161/97‑1 (MQ=255)                                                                                                                  gagagaAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGt  <  2:1023430/139‑1 (MQ=255)                                                                                                                      gaAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAg                                        <  2:1397910/98‑1 (MQ=255)                                                                                                                      gaAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAg                                        >  1:1397910/1‑98 (MQ=255)                                                                                                                       aaaGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGAttt            >  2:461155/1‑124 (MQ=255)                                                                                                                       aaaGGAGGCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGAttt            <  1:461155/124‑1 (MQ=255)                                                                                                                            aggCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCgg                                                  <  2:517011/81‑1 (MQ=255)                                                                                                                            aggCGCTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCgg                                                  >  1:517011/1‑81 (MQ=255)                                                                                                                                |                                                                                                                              CGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAA‑GGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGT  >  NC_000913/911005‑911256

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