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%	134.7 / NA	42	*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 (19/23);  total (19/23)

GCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGC  >  NC_000913/1926649‑1926911                                                                                                                                      |                                                                                                                                  gCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGc                                                                                                                              <  1:175294/139‑1 (MQ=255)      tACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCg                                                                                                                         >  1:1332390/1‑139 (MQ=255)       aCCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGt                                                                                                                        <  2:422941/139‑1 (MQ=255)          gTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGc                                                                                                                     <  2:1332390/139‑1 (MQ=255)          gTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGc                                                                                                                     <  1:717072/139‑1 (MQ=255)                 aGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAAt                                                                                                              >  2:1130230/1‑139 (MQ=255)                 aGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAAt                                                                                                              >  1:1373892/1‑139 (MQ=255)                      tttCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTcc                                                                                                         <  2:330165/139‑1 (MQ=255)                                                  cTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCatat                                                                             >  2:765445/1‑139 (MQ=255)                                                    aTCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAgcgc                                                                                             >  2:746180/1‑121 (MQ=255)                                                    aTCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAgcgc                                                                                             <  1:746180/121‑1 (MQ=255)                                                      cTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCa                                                                         >  1:1497234/1‑139 (MQ=255)                                                           cTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCt                                                                    <  1:765445/139‑1 (MQ=255)                                                           cTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCt                                                                    <  1:257037/139‑1 (MQ=255)                                                            tGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTa                                                                   <  1:157240/139‑1 (MQ=255)                                                            tGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTa                                                                   <  2:1497234/139‑1 (MQ=255)                                                                   ttAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTc                                                            >  2:1468043/1‑139 (MQ=255)                                                                        aaCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGt                                                                                                                 >  2:1311209/1‑81 (MQ=255)                                                                        aaCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGt                                                                                                                 <  1:1311209/81‑1 (MQ=255)                                                                            tGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATGTGGAAAATACTTAATCCg                                                                                                                               >  2:1068545/1‑63 (MQ=255)                                                                            tGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCg                                                                                                                               <  1:1068545/63‑1 (MQ=255)                                                                               aCAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTaaa                                                >  1:184673/1‑139 (MQ=255)                                                                                   cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGcgcg                                                                                            >  2:263572/1‑91 (MQ=255)                                                                                   cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGcgcg                                                                                            <  1:263572/91‑1 (MQ=255)                                                                                   cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGa                                            >  1:831120/1‑139 (MQ=255)                                                                                   cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGa                                            >  1:101867/1‑139 (MQ=255)                                                                                   cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGa                                            >  2:1144688/1‑139 (MQ=255)                                                                                        gCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGc                                                                                                                          >  2:1379408/1‑56 (MQ=255)                                                                                        gCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGc                                                                                                                          <  1:1379408/56‑1 (MQ=255)                                                                                        gCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGaaaa                                                                                   <  2:615202/95‑1 (MQ=255)                                                                                        gCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGaaaa                                                                                   >  1:615202/1‑95 (MQ=255)                                                                                                 tGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGcc                              <  2:320027/139‑1 (MQ=255)                                                                                                  gCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCa                             <  2:184673/139‑1 (MQ=255)                                                                                                        aaCACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCt                       <  2:1373892/139‑1 (MQ=255)                                                                                                        aaCACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCt                       <  1:1130230/139‑1 (MQ=255)                                                                                                                  cATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGt             <  2:831120/139‑1 (MQ=255)                                                                                                                             aaaTACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGc  <  1:1144688/139‑1 (MQ=255)                                                                                                                                 aCTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTa                                                                   <  1:411698/70‑1 (MQ=255)                                                                                                                                 aCTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTa                                                                   >  2:411698/1‑70 (MQ=255)                                                                                                                                 aCTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACag     >  1:702109/1‑132 (MQ=255)                                                                                                                                 aCTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACag     <  2:702109/132‑1 (MQ=255)                                                                                                                                   ttAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACg                                                                 >  1:287757/1‑70 (MQ=255)                                                                                                                                      |                                                                                                                                  GCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGC  >  NC_000913/1926649‑1926911

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