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%	153.6 / NA	48	*610* (TAG→TAA)	recQ	ATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/24);  total (24/24)

GATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCCC  >  NC_000913/4007558‑4007819                                                                                                                                        |                                                                                                                               gaTGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATca                                                                                                                             >  1:1014290/1‑139 (MQ=255) gaTGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATca                                                                                                                             >  1:765505/1‑139 (MQ=255)  aTGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                            <  1:494714/139‑1 (MQ=255)    ggCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCACGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcagca                                                                                                                          <  1:874442/139‑1 (MQ=255)      cTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCagcagc                                                                                                                        >  1:837043/1‑139 (MQ=255)               tGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAg                                                                                                               >  1:348497/1‑139 (MQ=255)                gCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGt                                                                                                              >  2:621576/1‑139 (MQ=255)                gCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGt                                                                                                              >  1:920118/1‑139 (MQ=255)                  cGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCAT‑AAAAAGTGcc                                                                                                           <  2:1014290/139‑1 (MQ=255)                   gATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                           <  2:348497/139‑1 (MQ=255)                        ccGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTAt                                                                                                      >  2:894352/1‑139 (MQ=255)                           ccAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGaa                                                                                                   >  1:1053240/1‑139 (MQ=255)                             aGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGa                                                                                                 <  1:775093/139‑1 (MQ=255)                             aGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGa                                                                                                 <  1:247594/139‑1 (MQ=255)                               cGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACt                                                                                               <  2:38402/139‑1 (MQ=255)                               cGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACt                                                                                               >  2:729291/1‑139 (MQ=255)                                         aGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGt                                                                                     >  2:674992/1‑139 (MQ=255)                                               aaCGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTcccc                                                                               <  2:837043/139‑1 (MQ=255)                                                 cGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                             <  1:361220/139‑1 (MQ=255)                                                  ggCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTCATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcgc                                                                            <  1:729291/139‑1 (MQ=255)                                                      ttGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGt                                                                        <  2:1203681/139‑1 (MQ=255)                                                        gggATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCa                                                                      <  2:920118/139‑1 (MQ=255)                                                              cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGa                                                                                                 <  2:145197/106‑1 (MQ=255)                                                              cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGa                                                                                                 >  1:145197/1‑106 (MQ=255)                                                              cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgta                                                                >  1:650141/1‑139 (MQ=255)                                                                 aaGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgtatgt                                                             <  2:765505/139‑1 (MQ=255)                                                                  aGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGtt                                                            <  1:546212/139‑1 (MQ=255)                                                                        aaCGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcgc                                                                            <  1:493124/117‑1 (MQ=255)                                                                        aaCGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcgc                                                                            >  2:493124/1‑117 (MQ=255)                                                                             tttGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTtc                                                 >  2:255066/1‑139 (MQ=255)                                                                                     aCCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTc                                         >  2:948728/1‑139 (MQ=255)                                                                                                cgcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGca                              >  2:1034136/1‑139 (MQ=255)                                                                                                cgcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGca                              <  2:308103/139‑1 (MQ=255)                                                                                                  cTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                           <  1:61134/60‑1 (MQ=255)                                                                                                  cTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                           >  2:61134/1‑60 (MQ=255)                                                                                                  cTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTtctc                                               <  2:112579/120‑1 (MQ=255)                                                                                                  cTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTtctc                                               >  1:112579/1‑120 (MQ=255)                                                                                                     aTTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATaaaa                                                                                                                  <  1:1028152/50‑1 (MQ=255)                                                                                                     aTTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATaaaa                                                                                                                  >  2:1028152/1‑50 (MQ=255)                                                                                                       tCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATtgtg                       <  1:674992/139‑1 (MQ=255)                                                                                                          tGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGtt                                                            <  2:732129/99‑1 (MQ=255)                                                                                                          tGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGtt                                                            >  1:732129/1‑99 (MQ=255)                                                                                                                 gTTGATGGCGATGACGAAGAATAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGa             >  1:716307/1‑139 (MQ=255)                                                                                                                         cGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCg                                        >  2:532915/1‑104 (MQ=255)                                                                                                                         cGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCg                                        <  1:532915/104‑1 (MQ=255)                                                                                                                         cGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGt     >  1:756825/1‑139 (MQ=255)                                                                                                                           aTGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTcc   >  2:670246/1‑139 (MQ=255)                                                                                                                            tGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTccc  <  2:952608/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                               GATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCCC  >  NC_000913/4007558‑4007819

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