Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	746,063	A→G	Q43Q (CAA→CAG)	nei →	endonuclease VIII and 5‑formyluracil/5‑hydroxymethyluracil DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	746,063	0	A	G	100.0%	142.6 / NA	38	Q43Q (CAA→CAG)	nei	endonuclease VIII and 5‑formyluracil/5‑hydroxymethyluracil DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/21);  total (17/21)

TCACCATGCCTGAAGGCCCGGAGATCCGCCGTGCAGCGGATAACCTGGAGGCGGCGATCAAAGGCAAACCACTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAACACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGACCACGCGA  >  NC_000913/745930‑746198                                                                                                                                      |                                                                                                                                        tCACCATGCCTGAAGGCCCGGAGATCCGCCGTGCAGCGGATAACCTGGAGGCGGCGATCAAAGGCAAACCACTAACTGATGTCTGGTTTGCCTGCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGt                                                                                                                                    >  1:714411/1‑139 (MQ=255)   aCCATGCCTGAAGGCCCGGAGATCCGCCGTGCAGCGGATAACCTGGAGGCGGCGATCAAAGGCAAACCACTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTa                                                                                                                                  <  2:714411/139‑1 (MQ=255)    ccATGCCTGAAGGCCCGGAGATCCGCCGTGCAGCGGATAACCTGGAGGCGGCGATCAAAGGCAAAGCACTAACTGATGTCTGGTTTGCCTTCCCACAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTAc                                                                                                                                 <  1:804910/139‑1 (MQ=255)    ccATGCCTGAAGGCCCGGAGATCCGCCGTGCAGCGGATAACCTGGAGGCGGCGATCAAAGGCAAACCACTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTAc                                                                                                                                 <  2:615315/139‑1 (MQ=255)    ccATGCCTGAAGGCCCGGAGATCCGCCGTGCAGCGGATAACCTGGAGGCGGCGATCAAAGGCAAACCACTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTAc                                                                                                                                 <  2:786054/139‑1 (MQ=255)                  ccGGAGATCCGCCGTGCAGCGGATAACCTGGAGGCGGCGATCAAAGGCAAACCACTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAAcgc                                                                                                                   >  2:77930/1‑139 (MQ=255)                   cGGAGATCCGCCGTGCAGCGGATAACCTGGAGGCGGCGATCAAAGGCAAACCACTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAAcgcg                                                                                                                  >  2:463426/1‑139 (MQ=255)                   cGGAGATCCGCCGTGCAGCGGATAACCTGGAGGCGGCGATCAAAGGCAAACCACTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAAcgcg                                                                                                                  <  2:732098/139‑1 (MQ=255)                   cGGAGATCCGCCGTGCAGCGGATAACCTGGAGGCGGCGATCAAAGGCAAACCACTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAAcgcg                                                                                                                  <  1:463426/139‑1 (MQ=255)                           cgccgTGCAGCGGATAACCTGGAGGCGGCGATCAAATGCAAACCACTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGt                                                                                                              >  1:462070/1‑135 (MQ=255)                           cgccgTGCAGCGGATAACCTGGAGGCGGCGATCAAAGGCAAACCACTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGt                                                                                                              <  2:462070/135‑1 (MQ=255)                                  cAGCGGATAACCTGGAGGCGGCGATCAAAGGCAAACCACTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCgttgtt                                                                                                   >  2:954399/1‑139 (MQ=255)                                   aGCGGATAACCTGGAGGCGGCGATCAAAGGCAAACCACTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTa                                                                                                  >  1:786627/1‑139 (MQ=255)                                                      gcgATCAAAGGCAAACCACTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACt                                                                               <  1:954399/139‑1 (MQ=255)                                                      gcgATCAAAGGCAAACCACTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACt                                                                               >  2:183293/1‑139 (MQ=255)                                                           cAAAGGCAAACCACTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACg                                                                          <  2:786627/139‑1 (MQ=255)                                                                    aCCACTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAgca                                                                                                                                       >  1:916918/1‑69 (MQ=255)                                                                       aCTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCAt                                                              <  2:375268/139‑1 (MQ=255)                                                                        cTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGt                                                                                                              >  1:813728/1‑90 (MQ=255)                                                                        cTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGt                                                                                                              <  2:813728/90‑1 (MQ=255)                                                                         tAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTAc                                                                                                                                 <  1:885333/70‑1 (MQ=255)                                                                         tAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTAc                                                                                                                                 >  2:885333/1‑70 (MQ=255)                                                                         tAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATaa                                                            <  2:294771/139‑1 (MQ=255)                                                                          aaCTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACg                                                                          >  1:850099/1‑124 (MQ=255)                                                                          aaCTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACg                                                                          <  2:850099/124‑1 (MQ=255)                                                                          aaCTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAAt                                                           <  1:183293/139‑1 (MQ=255)                                                                                      gTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGc                                               >  2:1066411/1‑139 (MQ=255)                                                                                                       tAAAACCTTATCAATCACAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGa                              <  2:1046156/139‑1 (MQ=255)                                                                                                                atcaatcaCAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGctc                                                                       >  1:820813/1‑89 (MQ=255)                                                                                                                atcaatcaCAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGctc                                                                       <  2:820813/89‑1 (MQ=255)                                                                                                                atcaatcaCAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTAc                                                  <  2:581150/110‑1 (MQ=255)                                                                                                                atcaatcaCAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTAc                                                  >  1:581150/1‑110 (MQ=255)                                                                                                                 tcaatcaCAACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGaa                    <  2:586771/139‑1 (MQ=255)                                                                                                                       acaACTTATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCg              >  1:439703/1‑139 (MQ=255)                                                                                                                             tATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGaa                    >  1:937449/1‑127 (MQ=255)                                                                                                                             tATCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGaa                    <  2:937449/127‑1 (MQ=255)                                                                                                                              aTCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGACCa       <  2:51605/139‑1 (MQ=255)                                                                                                                               tCGGTCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGACCAc      >  2:655889/1‑139 (MQ=255)                                                                                                                                   tCAGCACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGACCACGCGa  >  1:729103/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                        TCACCATGCCTGAAGGCCCGGAGATCCGCCGTGCAGCGGATAACCTGGAGGCGGCGATCAAAGGCAAACCACTAACTGATGTCTGGTTTGCCTTCCCGCAGTTAAAACCTTATCAATCACAACTTATCGGTCAACACGTTACCCATGTGGAAACGCGTGGTAAGGCGTTGTTAACTCATTTTTCCAACGACTTAACGCTCTACAGCCATAATCAGCTTTACGGCGTCTGGCGCGTGGTTGATACCGGCGAAGAGCCGCAGACCACGCGA  >  NC_000913/745930‑746198

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