Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,926,782	G→A	*221* (TAG→TAA)  D686D (GAC→GAT)	exoX → ptrB ←	exodeoxyribonuclease 10, DNA exonuclease X protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,926,782	0	G	A	100.0%	205.9 / NA	64	*221* (TAG→TAA)  D686D (GAC→GAT)	exoX ptrB	exodeoxyribonuclease 10, DNA exonuclease X protease II
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/35);  total (29/35)

TGGCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGA  >  NC_000913/1926647‑1926909                                                                                                                                        |                                                                                                                                tGGCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATcc                                                                                                                              >  1:117616/1‑139 (MQ=255)          ccGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTa                                                                                                                     <  2:988036/139‑1 (MQ=255)             tGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGt                                                                                                                      <  2:1200626/135‑1 (MQ=255)             tGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGt                                                                                                                      >  1:1200626/1‑135 (MQ=255)             tGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCg                                                                                                                  >  2:619762/1‑139 (MQ=255)               gCAAAGCGGTTTCCGACGTTGCCGAATGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCggg                                                                                                                <  2:1526167/139‑1 (MQ=255)               gCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCggg                                                                                                                >  1:1526167/1‑139 (MQ=255)                              aCGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTgcgc                                                                                                 >  1:739705/1‑139 (MQ=255)                                       aaCGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                        <  2:1011853/139‑1 (MQ=255)                                       aaCGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAAGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                        >  2:613353/1‑139 (MQ=255)                                           cgATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGa                                                                                    <  2:26205/139‑1 (MQ=255)                                            gATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGaa                                                                                   <  2:117616/139‑1 (MQ=255)                                                   gCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAAGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCata                                                                            <  1:613353/139‑1 (MQ=255)                                                       tCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTcc                                                                        >  2:61974/1‑139 (MQ=255)                                                       tCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTcc                                                                        <  2:848267/139‑1 (MQ=255)                                                                gTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACg                                                               <  1:619762/139‑1 (MQ=255)                                                                          aaCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCgg                                                                                                                 >  1:1301466/1‑79 (MQ=255)                                                                          aaCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCgg                                                                                                                 <  2:1301466/79‑1 (MQ=255)                                                                              tGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCa                                                                                                                           <  2:883983/65‑1 (MQ=255)                                                                              tGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCa                                                                                                                           >  1:883983/1‑65 (MQ=255)                                                                               ggACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTa                                                >  2:398646/1‑139 (MQ=255)                                                                                gACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTaa                                               >  1:1206120/1‑139 (MQ=255)                                                                                 aCAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTaaa                                              <  2:739705/139‑1 (MQ=255)                                                                                     cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAACTACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGa                                          >  2:1085356/1‑139 (MQ=255)                                                                                     cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATcc                                                                                                                              >  2:713036/1‑55 (MQ=255)                                                                                     cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATcc                                                                                                                              <  1:713036/55‑1 (MQ=255)                                                                                     cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTa                                                                                                                     >  2:615398/1‑64 (MQ=255)                                                                                     cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTa                                                                                                                     <  1:615398/64‑1 (MQ=255)                                                                                     cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                        >  1:343776/1‑93 (MQ=255)                                                                                     cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                        <  2:343776/93‑1 (MQ=255)                                                                                     cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGAt                                                   <  2:661194/130‑1 (MQ=255)                                                                                     cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGAt                                                   >  1:661194/1‑130 (MQ=255)                                                                                     cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGa                                          >  2:486788/1‑139 (MQ=255)                                                                                     cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGa                                          >  1:19889/1‑139 (MQ=255)                                                                                     cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGa                                          >  2:1372659/1‑139 (MQ=255)                                                                                     cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGa                                          >  1:1517170/1‑139 (MQ=255)                                                                                     cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGa                                          >  1:1277583/1‑139 (MQ=255)                                                                                       tGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGAcc                                        <  2:1277583/139‑1 (MQ=255)                                                                                          gCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTA‑‑GGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGAtt                                   >  2:1117738/1‑139 (MQ=255)                                                                                              ggAGCTGCGTTTAACACTGAAACATTATCGGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTa                                 <  1:263139/139‑1 (MQ=255)                                                                                               gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTa                                                                                                                     >  2:1052117/1‑54 (MQ=255)                                                                                               gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTa                                                                                                                     <  1:1052117/54‑1 (MQ=255)                                                                                               gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGc                                                                                                                   <  2:743160/56‑1 (MQ=255)                                                                                               gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGc                                                                                                                   >  1:743160/1‑56 (MQ=255)                                                                                               gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTAc                                >  1:1511178/1‑139 (MQ=255)                                                                                                  cTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTA‑‑GGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCa                           <  1:1117738/139‑1 (MQ=255)                                                                                                   tGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGcc                            <  1:61974/139‑1 (MQ=255)                                                                                                     cGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                        >  2:1248366/1‑77 (MQ=255)                                                                                                     cGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                        <  1:1248366/77‑1 (MQ=255)                                                                                                        ttAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGt                                                        >  1:955544/1‑106 (MQ=255)                                                                                                        ttAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGt                                                        <  2:955544/106‑1 (MQ=255)                                                                                                        ttAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGc                       <  1:1085356/139‑1 (MQ=255)                                                                                                          aaCACTGAAACATTATCTGGAAAATACTTAATCCGCGGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGcc                      <  2:1383627/138‑1 (MQ=255)                                                                                                          aaCACTGAAACATTATCTGGAAAATACTTAATCCGCGGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGcc                      >  1:1383627/1‑138 (MQ=255)                                                                                                          aaCACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCt                     <  1:1372659/139‑1 (MQ=255)                                                                                                               tGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                <  2:1206120/139‑1 (MQ=255)                                                                                                               tGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                <  1:302528/139‑1 (MQ=255)                                                                                                               tGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                <  2:19889/139‑1 (MQ=255)                                                                                                               tGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                <  2:1517170/139‑1 (MQ=255)                                                                                                                    cATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGa                                                                                    >  1:1183322/1‑66 (MQ=255)                                                                                                                    cATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGa                                                                                    <  2:1183322/66‑1 (MQ=255)                                                                                                                             gAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACaga  <  1:398646/139‑1 (MQ=255)                                                                                                                                   aCTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCa                                                                              >  1:1078715/1‑57 (MQ=255)                                                                                                                                   aCTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCa                                                                              <  2:1078715/57‑1 (MQ=255)                                                                                                                                     ttAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                        <  2:615029/45‑1 (MQ=255)                                                                                                                                        |                                                                                                                                TGGCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGA  >  NC_000913/1926647‑1926909

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