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%	199.5 / NA	62	*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 (32/30);  total (32/30)

TGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATA  >  NC_000913/1926659‑1926915                                                                                                                            |                                                                                                                                      tGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCg                                                                                                                        >  2:1137423/1‑139 (MQ=255)  ggCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCgg                                                                                                                       <  2:1448308/139‑1 (MQ=255)      aaGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTaa                                                                                                                   >  1:108333/1‑139 (MQ=255)                      tGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCg                                                                                                                                   >  1:71304/1‑107 (MQ=255)                      tGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCg                                                                                                                                   <  2:71304/107‑1 (MQ=255)                                            cTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCAT‑TTCCAt                                                                            >  1:1647886/1‑139 (MQ=255)                                                 cTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCt                                                                        <  2:108333/139‑1 (MQ=255)                                                 cTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATACGCAGGCGTAGCGGGTAATGTTCACTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCt                                                                        >  2:1080755/1‑139 (MQ=255)                                                    gTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACg                                                                     >  2:417828/1‑139 (MQ=255)                                                     ttATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGc                                                                    >  1:662516/1‑139 (MQ=255)                                                     ttATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGc                                                                    >  2:710904/1‑139 (MQ=255)                                                       aTTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCt                                                                  <  1:1080755/139‑1 (MQ=255)                                                       aTTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCt                                                                  <  1:417828/139‑1 (MQ=255)                                                          taataaCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCAT‑TTCCATCGCTACGCCTTCGt                                                              <  2:1647886/139‑1 (MQ=255)                                                              aaCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGc                                                                                                                              <  2:1028028/72‑1 (MQ=255)                                                              aaCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGc                                                                                                                              >  1:1028028/1‑72 (MQ=255)                                                              aaCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACg                                                           >  2:957895/1‑139 (MQ=255)                                                               aCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTcc                                                                                                             <  1:1234769/88‑1 (MQ=255)                                                               aCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTcc                                                                                                             >  2:1234769/1‑88 (MQ=255)                                                                  tGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATg                                                                                                                 <  1:1222871/81‑1 (MQ=255)                                                                  tGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATg                                                                                                                 >  2:1222871/1‑81 (MQ=255)                                                                  tGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGAttt                                                       >  1:586400/1‑139 (MQ=255)                                                                  tGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGAttt                                                       >  2:357700/1‑139 (MQ=255)                                                                   ggACGGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTcc                                                                                                             <  2:1026339/84‑1 (MQ=255)                                                                   ggACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTa                                                      <  1:957895/139‑1 (MQ=255)                                                                   ggACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCAGGTAATGTTcc                                                                                                             >  1:1026339/1‑84 (MQ=255)                                                                         cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTgcg                                                                                                        <  2:1523148/83‑1 (MQ=255)                                                                         cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTgcg                                                                                                        >  1:1523148/1‑83 (MQ=255)                                                                         cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCt                                                                  >  1:695062/1‑121 (MQ=255)                                                                         cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCt                                                                  <  2:695062/121‑1 (MQ=255)                                                                           tGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGAcc                                              <  2:1303819/139‑1 (MQ=255)                                                                                 cGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTc                                                                               <  1:406024/100‑1 (MQ=255)                                                                                 cGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTc                                                                               >  2:406024/1‑100 (MQ=255)                                                                                   gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                              <  2:1108700/83‑1 (MQ=255)                                                                                   gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                              >  1:1108700/1‑83 (MQ=255)                                                                                   gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGa                                                                                          >  1:260461/1‑87 (MQ=255)                                                                                   gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGa                                                                                          <  2:260461/87‑1 (MQ=255)                                                                                   gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTAc                                      >  1:926822/1‑139 (MQ=255)                                                                                        gCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCa                                 >  2:668192/1‑139 (MQ=255)                                                                                         cGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTcc                                                                              >  2:1269545/1‑93 (MQ=255)                                                                                         cGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTcc                                                                              <  1:1269545/93‑1 (MQ=255)                                                                                         cGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGAcc                                              <  2:977857/125‑1 (MQ=255)                                                                                         cGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGAcc                                              >  1:977857/1‑125 (MQ=255)                                                                                         cGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCAt                                >  1:506775/1‑139 (MQ=255)                                                                                         cGTTTAACACTGAAACATGATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCAt                                <  1:1137423/139‑1 (MQ=255)                                                                                              aaCACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCt                           <  1:710904/139‑1 (MQ=255)                                                                                                   tGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                      >  1:682207/1‑139 (MQ=255)                                                                                                   tGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                      <  1:357700/139‑1 (MQ=255)                                                                                                             tCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa            <  2:926822/139‑1 (MQ=255)                                                                                                                ggAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACag         <  1:668192/139‑1 (MQ=255)                                                                                                                ggAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACag         <  1:1656571/139‑1 (MQ=255)                                                                                                                 gAAAATACTTAATCCGCAGGCGTAGCGGGTAATGtt                                                                                                               <  2:1524776/36‑1 (MQ=255)                                                                                                                 gAAAATACTTAATCCGCAGGCGTAGCGGGTAATGtt                                                                                                               >  1:1524776/1‑36 (MQ=38)                                                                                                                 gAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCAt                                <  2:260398/115‑1 (MQ=255)                                                                                                                 gAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCAt                                >  1:260398/1‑115 (MQ=255)                                                                                                                   aaaTACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTcc                                                                              <  2:994173/67‑1 (MQ=255)                                                                                                                   aaaTACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTcc                                                                              >  1:994173/1‑67 (MQ=255)                                                                                                                   aaaTACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                      >  2:660163/1‑123 (MQ=255)                                                                                                                   aaaTACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                      <  1:660163/123‑1 (MQ=255)                                                                                                                       aCTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGAcc                                              <  1:1064513/95‑1 (MQ=255)                                                                                                                       aCTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGAcc                                              >  2:1064513/1‑95 (MQ=255)                                                                                                                       aCTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATa  >  1:1270698/1‑139 (MQ=255)                                                                                                                            |                                                                                                                                      TGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATA  >  NC_000913/1926659‑1926915

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