Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,007,693	G→A	*610* (TAG→TAA)	recQ →	ATP‑dependent DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,007,693	0	G	A	100.0%	209.1 / NA	65	*610* (TAG→TAA)	recQ	ATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/35);  total (30/35)

ATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCC  >  NC_000913/4007559‑4007818                                                                                                                                       |                                                                                                                              aTGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                           <  1:754783/139‑1 (MQ=255) aTGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCCGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                           <  2:1001349/139‑1 (MQ=255) aTGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCACAAGCTGGAACGCTTTGGCAAACCGTTTATGGCTCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                           <  1:1223705/139‑1 (MQ=255)          cAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATaaa                                                                                                                  >  2:1229233/1‑139 (MQ=255)             aTGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATaaaaaa                                                                                                               >  2:1688366/1‑139 (MQ=255)             aTGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATaaaaaa                                                                                                               >  1:259229/1‑139 (MQ=255)                  gATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          <  1:268117/139‑1 (MQ=255)                  gATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          <  1:1688366/139‑1 (MQ=255)                       ccGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          <  2:1370502/134‑1 (MQ=255)                       ccGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          >  1:1370502/1‑134 (MQ=255)                           cAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAg                                                                                                 <  1:1229233/139‑1 (MQ=255)                                  aTGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGt                                                                                          >  1:1101704/1‑139 (MQ=255)                                          cGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTAt                                                                                                     >  2:126839/1‑120 (MQ=255)                                          cGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTAt                                                                                                     <  1:126839/120‑1 (MQ=255)                                                cGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                            <  2:259229/139‑1 (MQ=255)                                                    gTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          >  1:1592657/1‑105 (MQ=255)                                                    gTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          <  2:1592657/105‑1 (MQ=255)                                                     ttGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGATTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGt                                                                       <  2:266443/139‑1 (MQ=255)                                                       gggATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCa                                                                     >  2:1278205/1‑139 (MQ=255)                                                        ggATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCaa                                                                    <  2:1101704/139‑1 (MQ=255)                                                             cttgcgcaggAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgta                                                               >  1:769014/9‑139 (MQ=255)                                                               cAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          >  1:1692680/1‑94 (MQ=255)                                                               cAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          <  2:1692680/94‑1 (MQ=255)                                                                 aGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGtt                                                           <  2:689094/139‑1 (MQ=255)                                                                     ggAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTgttg                                                       <  2:1316529/139‑4 (MQ=255)                                                                     ggAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATg                                                       <  1:444574/139‑1 (MQ=255)                                                                     ggAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATg                                                       <  1:1278205/139‑1 (MQ=255)                                                                              tGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAGGAGTAATCAGCAGCATAAAAAAGTg                                                                                                            <  1:249708/77‑1 (MQ=255)                                                                              tGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAGGAGTAATCAGCAGCATAAAAAAGTg                                                                                                            >  2:249708/1‑77 (MQ=255)                                                                                gCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                           <  2:964054/60‑1 (MQ=255)                                                                                gCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                           >  1:964054/1‑60 (MQ=255)                                                                                 cAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCAcc                                           >  1:483437/1‑139 (MQ=255)                                                                                    aCCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTc                                        >  2:1112886/1‑139 (MQ=255)                                                                                       gTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGcc                                     <  1:186040/139‑1 (MQ=255)                                                                                        tttATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          >  1:638542/1‑69 (MQ=255)                                                                                        tttATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          <  2:638542/69‑1 (MQ=255)                                                                                          tATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGa                                                         >  2:1124850/1‑116 (MQ=255)                                                                                          tATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGa                                                         <  1:1124850/116‑1 (MQ=255)                                                                                                  tGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTa                                                                                         >  1:875928/1‑76 (MQ=255)                                                                                                  tGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTa                                                                                         <  2:875928/76‑1 (MQ=255)                                                                                                    aTTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGca                             <  2:894351/134‑1 (MQ=255)                                                                                                    aTTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGca                             >  1:894351/1‑134 (MQ=255)                                                                                                         tGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          <  2:135373/52‑1 (MQ=255)                                                                                                         tGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          >  1:135373/1‑52 (MQ=255)                                                                                                            gcATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCTAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                            <  1:406413/79‑1 (MQ=255)                                                                                                            gcATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCTAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                            >  2:406413/1‑79 (MQ=255)                                                                                                              aTGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGgcg                   >  2:88121/1‑134 (MQ=255)                                                                                                              aTGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGgcg                   <  1:88121/134‑1 (MQ=255)                                                                                                                gTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgt                                                                >  2:316270/1‑87 (MQ=255)                                                                                                                gTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgt                                                                <  1:316270/87‑1 (MQ=255)                                                                                                                   gATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGcc         >  2:433057/1‑139 (MQ=255)                                                                                                                      ggCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATg             <  2:1490978/132‑1 (MQ=255)                                                                                                                      ggCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATg             >  1:1490978/1‑132 (MQ=255)                                                                                                                        cGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGgcg                   <  1:816629/124‑1 (MQ=255)                                                                                                                        cGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGgcg                   >  2:816629/1‑124 (MQ=255)                                                                                                                        cGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGt    <  2:483437/139‑1 (MQ=255)                                                                                                                          aTGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTcc  >  2:561666/1‑139 (MQ=255)                                                                                                                           tGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATg                                  <  2:854293/106‑1 (MQ=255)                                                                                                                           tGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATg                                  >  1:854293/1‑106 (MQ=255)                                                                                                                            gACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGa                                                                                                <  2:396012/43‑1 (MQ=255)                                                                                                                            gACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGa                                                                                                >  1:396012/1‑43 (MQ=255)                                                                                                                            gACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTa                                                                                         <  1:330500/50‑1 (MQ=255)                                                                                                                            gACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTa                                                                                         >  2:330500/1‑50 (MQ=255)                                                                                                                             aCGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgtatg                                                             >  1:354933/1‑77 (MQ=255)                                                                                                                             aCGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgtatg                                                             <  2:354933/77‑1 (MQ=255)                                                                                                                                       |                                                                                                                              ATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCC  >  NC_000913/4007559‑4007818

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