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

ATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCCCGAT  >  NC_000913/4007559‑4007822                                                                                                                                       |                                                                                                                                  aTGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                               <  1:989669/139‑1 (MQ=255) aTGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                               <  2:800905/139‑1 (MQ=255)         aCAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATaa                                                                                                                       >  2:1312973/1‑139 (MQ=255)           aGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATaaaa                                                                                                                     <  2:347385/139‑1 (MQ=255)              tGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAg                                                                                                                  >  1:1266234/1‑139 (MQ=255)               gCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGt                                                                                                                 >  1:1027868/1‑139 (MQ=255)                ccGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTg                                                                                                                <  2:694025/139‑1 (MQ=255)                ccGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTg                                                                                                                <  1:1312973/139‑1 (MQ=255)                        cGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATg                                                                                                        >  1:702375/1‑139 (MQ=255)                          ccAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGaa                                                                                                      <  2:486325/139‑1 (MQ=255)                          ccAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGaa                                                                                                      >  1:486325/1‑139 (MQ=255)                                   tGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTa                                                                                             <  2:522751/139‑1 (MQ=255)                                     cTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTaaa                                                                                           <  2:761669/139‑1 (MQ=255)                                               aCGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTccccc                                                                                 >  2:431795/1‑139 (MQ=255)                                                cGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                                <  2:1268798/139‑1 (MQ=255)                                                     ttGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGt                                                                           <  2:702375/139‑1 (MQ=255)                                                             cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                               <  2:616149/79‑1 (MQ=255)                                                             cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                               >  1:616149/1‑79 (MQ=255)                                                             cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgt                                                                    <  2:241631/138‑1 (MQ=255)                                                             cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgt                                                                    >  1:241631/1‑138 (MQ=255)                                                              gcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgtat                                                                  <  2:1266234/139‑1 (MQ=255)                                                               cAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                               >  2:338963/1‑77 (MQ=255)                                                               cAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                               >  1:1184060/1‑77 (MQ=255)                                                               cAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                               <  1:338963/77‑1 (MQ=255)                                                               cAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                               <  2:1184060/77‑1 (MQ=255)                                                               cAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGt                                                                           <  2:1243826/129‑1 (MQ=255)                                                               cAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGt                                                                           >  1:1243826/1‑129 (MQ=255)                                                                            tttGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAg                                                                                                                  <  2:847317/77‑1 (MQ=255)                                                                            tttGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAg                                                                                                                  >  1:847317/1‑77 (MQ=255)                                                                            tttGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTtc                                                    >  1:179949/1‑139 (MQ=255)                                                                              tGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTg                                                                                                                >  2:286155/1‑77 (MQ=255)                                                                              tGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTg                                                                                                                <  1:286155/77‑1 (MQ=255)                                                                                    aCCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTc                                            >  2:690993/1‑139 (MQ=255)                                                                                       gTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGcc                                         <  2:179949/139‑1 (MQ=255)                                                                                        tttATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGa                                                                             <  2:489065/102‑1 (MQ=255)                                                                                        tttATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGa                                                                             >  1:489065/1‑102 (MQ=255)                                                                                                gcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGcac                                >  2:1030760/1‑139 (MQ=255)                                                                                                 cTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGc                                                                                                               <  2:337638/59‑1 (MQ=255)                                                                                                 cTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGc                                                                                                               >  1:337638/1‑59 (MQ=255)                                                                                                 cTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGcaca                               >  1:90671/1‑139 (MQ=255)                                                                                                    aTTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGtt                                                         >  1:600245/1‑110 (MQ=255)                                                                                                    aTTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGtt                                                         <  2:600245/110‑1 (MQ=255)                                                                                                         tGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGaa                                                                                                      <  2:140843/60‑1 (MQ=255)                                                                                                         tGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGaa                                                                                                      >  1:140843/1‑60 (MQ=255)                                                                                                                gTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGt                                                                           >  2:401260/1‑80 (MQ=255)                                                                                                                gTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGt                                                                           <  1:401260/80‑1 (MQ=255)                                                                                                                    aTGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGccc            >  2:1017021/1‑139 (MQ=255)                                                                                                                          aTGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTgg                         <  2:543346/120‑1 (MQ=255)                                                                                                                          aTGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTgg                         >  1:543346/1‑120 (MQ=255)                                                                                                                          aTGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTcc      >  1:832946/1‑139 (MQ=255)                                                                                                                          aTGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCAGCGAGTCAAATGTATGTTGATGTTATTTCTCGCCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTcc      >  2:979930/1‑139 (MQ=255)                                                                                                                           tGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                                <  2:1127872/64‑1 (MQ=255)                                                                                                                           tGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                                >  1:1127872/1‑64 (MQ=255)                                                                                                                           tGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTccc     >  2:1346075/1‑139 (MQ=255)                                                                                                                             aCGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCCCGa   >  1:295329/1‑139 (MQ=255)                                                                                                                              cGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCCCGAt  <  1:431795/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                  ATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCCCGAT  >  NC_000913/4007559‑4007822

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