Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,961,917	G→C	G414A (GGC→GCC)	rep →	DNA helicase and single‑stranded DNA‑dependent ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,961,917	0	G	C	100.0%	230.3 / NA	65	G414A (GGC→GCC)	rep	DNA helicase and single‑stranded DNA‑dependent ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (38/27);  total (38/27)

CTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGACCCGCTTCACTC  >  NC_000913/3961791‑3962048                                                                                                                               |                                                                                                                                    cTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGATTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCt                                                                                                                         >  2:1134887/1‑139 (MQ=255)      ggtACGTCGTTTTTCTCTCGTCCTGAAATAAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGaaaa                                                                                                                    >  2:990851/1‑139 (MQ=255)        tACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAg                                                                                                                  >  2:1223685/1‑139 (MQ=255)                ttttCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGa                                                                                                          <  1:990851/139‑1 (MQ=255)                      ctcGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCga                                                                                                                                <  2:1327527/111‑2 (MQ=255)                      ctcGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCga                                                                                                                                >  1:1327527/1‑110 (MQ=255)                         gTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGAt                                                                                                 >  1:541276/1‑139 (MQ=255)                         gTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGAt                                                                                                 >  2:487407/1‑139 (MQ=255)                          tCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCTGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATg                                                                                                >  1:572741/1‑139 (MQ=255)                             tGAAATCGAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcg                                                                                             <  1:1223685/139‑1 (MQ=255)                              gAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCTGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcgc                                                                                            <  2:572741/139‑1 (MQ=255)                              gAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTAc                                                                                                                            <  1:1195507/107‑1 (MQ=255)                              gAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcgc                                                                                            <  2:541276/139‑1 (MQ=255)                              gAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTCCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTAc                                                                                                                            >  2:1195507/1‑107 (MQ=255)                                 aTCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCa                                                                                         >  2:356013/1‑139 (MQ=255)                                                     tCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGcc                                                                     >  2:1292944/1‑139 (MQ=255)                                                        gcgcGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcg                                                                                             <  2:749347/112‑1 (MQ=255)                                                        gcgcGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcg                                                                                             >  1:749347/1‑112 (MQ=255)                                                          gcgTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGa                                                                                                                       >  2:967278/1‑84 (MQ=255)                                                          gcgTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGa                                                                                                                       <  1:967278/84‑1 (MQ=255)                                                          gcgTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATa                                                                                      >  1:1120197/1‑117 (MQ=255)                                                          gcgTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATa                                                                                      <  2:1120197/117‑1 (MQ=255)                                                          gcgTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCtt                                                                <  1:887593/139‑1 (MQ=255)                                                           cgTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCttt                                                               >  1:90379/1‑139 (MQ=255)                                                              gCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGTTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGat                                                            >  2:113641/1‑139 (MQ=255)                                                              gCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGat                                                            >  2:855442/1‑139 (MQ=255)                                                              gCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATAGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAATGCATGTTTACCGCCAGCTTTGat                                                            >  2:993860/1‑139 (MQ=255)                                                               cTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGata                                                           <  1:356013/139‑1 (MQ=255)                                                                        ccGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGa                                                  >  1:133259/1‑139 (MQ=255)                                                                            aCGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCa                                              >  2:337031/1‑139 (MQ=255)                                                                            aCGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCa                                              >  2:489894/1‑139 (MQ=255)                                                                                 gACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCtt                                                                <  1:383509/116‑1 (MQ=255)                                                                                 gACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCtt                                                                >  2:383509/1‑116 (MQ=255)                                                                                      cgcATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCa                                              >  2:71652/1‑129 (MQ=255)                                                                                      cgcATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCa                                              <  1:71652/129‑1 (MQ=255)                                                                                      cgcATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGc                                    >  2:473516/1‑139 (MQ=255)                                                                                         aTTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGa                                 <  1:1134887/139‑1 (MQ=255)                                                                                           ttCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACg                               <  2:90379/139‑1 (MQ=255)                                                                                            tCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGat                                                            <  2:854519/109‑1 (MQ=255)                                                                                            tCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGat                                                            >  1:854519/1‑109 (MQ=255)                                                                                            tCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGt                              <  2:133259/139‑1 (MQ=255)                                                                                             cTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATa                                                                                      <  2:488685/82‑1 (MQ=255)                                                                                             cTGCGTATCGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATa                                                                                      >  1:488685/1‑82 (MQ=255)                                                                                                     cGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcgcg                                                                                           >  1:3985/1‑69 (MQ=255)                                                                                                     cGTTAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcgcg                                                                                           <  2:3985/69‑1 (MQ=255)                                                                                                        tAACACGCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCa                  <  1:1292944/139‑1 (MQ=255)                                                                                                              gCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAg                                     <  2:309205/114‑1 (MQ=255)                                                                                                              gCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAg                                     >  1:309205/1‑114 (MQ=255)                                                                                                              gCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGtt                          >  2:1235005/1‑125 (MQ=255)                                                                                                              gCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGtt                          <  1:1235005/125‑1 (MQ=255)                                                                                                              gCCGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGAcc            >  1:948076/1‑139 (MQ=255)                                                                                                               ccGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcg                                                                                             <  2:417353/57‑1 (MQ=255)                                                                                                               ccGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGAcg                                                                                             >  1:417353/1‑57 (MQ=255)                                                                                                               ccGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGc                                                                      <  2:1367353/80‑1 (MQ=255)                                                                                                               ccGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGc                                                                      >  1:1367353/1‑80 (MQ=255)                                                                                                               ccGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCtt                                       >  1:844400/1‑111 (MQ=255)                                                                                                               ccGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCtt                                       <  2:844400/111‑1 (MQ=255)                                                                                                               ccGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGAccc           >  2:1159578/1‑139 (MQ=255)                                                                                                               ccGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGAccc           >  1:1382317/1‑139 (MQ=255)                                                                                                               ccGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGAccc           >  2:1034232/1‑139 (MQ=255)                                                                                                                cGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACg                               >  1:1125341/1‑118 (MQ=255)                                                                                                                cGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGACCCg          >  2:1041136/1‑139 (MQ=255)                                                                                                                cGAAGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGACCCg          >  1:266494/1‑139 (MQ=255)                                                                                                                cGAAGCGAGAGATTGCCCCGACTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACg                               <  2:1125341/118‑1 (MQ=255)                                                                                                                   aGCGAGAGATTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGACCCGCtt       <  2:266494/139‑1 (MQ=255)                                                                                                                        gagaTTGCCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGACCCGCTtcactc  >  1:1475198/1‑139 (MQ=255)                                                                                                                               |                                                                                                                                    CTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGACAGCGCATTTCTGCGTATCGTTAACACGCCGAAGCGAGAGATTGGCCCGGCTACGCTGAAAAAGCTGGGTGAGTGGGCGATGACGCGCAATAAAAGCATGTTTACCGCCAGCTTTGATATGGGCCTGAGTCAGACGCTTAGCGGACGTGGTTATGAAGCATTGACCCGCTTCACTC  >  NC_000913/3961791‑3962048

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