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

GATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTC  >  NC_000913/4007558‑4007817                                                                                                                                        |                                                                                                                             gaTGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATca                                                                                                                           >  1:1362945/1‑139 (MQ=255)  gtGGCTGAACAGATGCCGAGCAGCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCGTTGGCAAACCGATTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGAGGAAGAGTAATcag                                                                                                                          <  2:947114/138‑1 (MQ=255)  aTGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                          <  2:1036938/139‑1 (MQ=255)           cAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCAGGTTGATGGCGATGACGAAGAGTAATCAGCAGCATaaa                                                                                                                 >  2:977207/1‑139 (MQ=255)               tGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAg                                                                                                             >  2:1478209/1‑139 (MQ=255)                         cGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATg                                                                                                   >  2:910890/1‑139 (MQ=255)                                      cTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                         >  1:434898/1‑120 (MQ=255)                                      cTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                         <  2:434898/120‑1 (MQ=255)                                           cGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                         <  1:626263/115‑1 (MQ=255)                                           cGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                         >  2:626263/1‑115 (MQ=255)                                                      ttGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGt                                                                      <  2:1362945/139‑1 (MQ=255)                                                      ttGGGATGCGCAAGCTGGAACGATTTGGCAAAGCGTTTATGTCACTGATTCGTGCGCATGTTGATGTCGATGACGAGGAGCAATCAGCAGTATAAAAAAGAGCCAGTATGAAGACTCCGTAAACGCTTCCCCCGCGAGt                                                                      <  2:1229227/139‑1 (MQ=255)                                                              cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                          >  2:1013592/1‑79 (MQ=255)                                                              cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                          <  1:1013592/79‑1 (MQ=255)                                                              cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgta                                                              >  1:87566/1‑139 (MQ=255)                                                               gcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGATTCCGTAAACGTTTCCCCCGCGAGTCAAatgtat                                                             <  2:126805/139‑1 (MQ=255)                                                                cAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcagca                                                                                                                        >  1:986759/1‑79 (MQ=255)                                                                cAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAatgtatg                                                            <  2:87566/139‑1 (MQ=255)                                                                cAAACTGGAACGCTTTCGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATctgca                                                                                                                        <  2:986759/79‑1 (MQ=255)                                                                      ggAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGTGATGACGAAGAGTAATCAGCAGCCTAAAAAAGTGCCAGTATGAAGACTCCCTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATg                                                      <  1:977207/139‑1 (MQ=255)                                                                        aaCGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                          >  1:30666/1‑69 (MQ=255)                                                                        aaCGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                          <  2:30666/69‑1 (MQ=255)                                                                             tttGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTtc                                               >  1:202741/1‑139 (MQ=255)                                                                             tttGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTtc                                               >  1:1484264/1‑139 (MQ=255)                                                                                     aCCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                           <  2:1168271/103‑1 (MQ=255)                                                                                     aCCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                           >  1:1168271/1‑103 (MQ=255)                                                                                        gTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGcc                                    <  1:1137029/139‑1 (MQ=255)                                                                                        gTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGcc                                    <  1:910890/139‑1 (MQ=255)                                                                                                cgcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACt                                                                                             >  1:762152/1‑74 (MQ=255)                                                                                                cgcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACt                                                                                             <  2:762152/74‑1 (MQ=255)                                                                                                cgcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTa                                                                                        <  1:1341838/79‑1 (MQ=255)                                                                                                cgcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTa                                                                                        >  2:1341838/1‑79 (MQ=255)                                                                                                 gcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                         >  1:1415162/1‑61 (MQ=255)                                                                                                 gcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                         <  2:1415162/61‑1 (MQ=255)                                                                                                 gcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                           >  2:154343/1‑91 (MQ=255)                                                                                                 gcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                           <  1:154343/91‑1 (MQ=255)                                                                                                 gcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGcac                           >  1:216122/1‑139 (MQ=255)                                                                                                  cTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcagcag                                                                                                                       <  2:1296092/46‑1 (MQ=255)                                                                                                  cTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcagcag                                                                                                                       >  1:1296092/1‑46 (MQ=255)                                                                                                  cTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGcc                                    <  1:151966/129‑1 (MQ=255)                                                                                                  cTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGcc                                    >  2:151966/1‑129 (MQ=255)                                                                                                       tCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTg                                                                                                           <  2:607791/53‑1 (MQ=255)                                                                                                       tCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTg                                                                                                           >  1:607791/1‑53 (MQ=255)                                                                                                          tGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                         >  1:427081/1‑52 (MQ=255)                                                                                                          tGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                         <  2:427081/52‑1 (MQ=255)                                                                                                                tGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATg            <  2:1064519/139‑1 (MQ=255)                                                                                                                         cGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGa                                                                                               <  1:1295403/47‑1 (MQ=255)                                                                                                                         cGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGa                                                                                               >  2:1295403/1‑47 (MQ=255)                                                                                                                         cGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATtg                       <  2:609515/119‑1 (MQ=255)                                                                                                                         cGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATtg                       >  1:609515/1‑119 (MQ=255)                                                                                                                         cGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGt   >  1:375256/1‑139 (MQ=255)                                                                                                                          gATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTc  <  2:202741/139‑1 (MQ=255)                                                                                                                           aTGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCa                                                                    >  2:689708/1‑72 (MQ=255)                                                                                                                           aTGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCa                                                                    <  1:689708/72‑1 (MQ=255)                                                                                                                           aTGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGtt                                                          >  2:559277/1‑82 (MQ=255)                                                                                                                           aTGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGtt                                                          <  1:559277/82‑1 (MQ=255)                                                                                                                           aTGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATtg                       >  1:1402837/1‑117 (MQ=255)                                                                                                                           aTGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATtg                       <  2:1402837/117‑1 (MQ=255)                                                                                                                            tGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGtt                                                          >  2:989426/1‑81 (MQ=255)                                                                                                                            tGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGtt                                                          <  1:989426/81‑1 (MQ=255)                                                                                                                                        |                                                                                                                             GATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTC  >  NC_000913/4007558‑4007817

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