Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,495,332	A→G	L293P (CTG→CCG)	ahr ←	broad specificity NADPH‑dependent aldehyde reductase, Zn‑containing

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,495,332	0	A	G	100.0%	178.8 / NA	50	L293P (CTG→CCG)	ahr	broad specificity NADPH‑dependent aldehyde reductase, Zn‑containing
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/25);  total (25/25)

GCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAA  >  NC_000913/4495214‑4495466                                                                                                                       |                                                                                                                                       gcgGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTa                                                                                                                    >  2:1106832/1‑139 (MQ=255)    gTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCa                                                                                                                 <  1:1030356/139‑1 (MQ=255)               tACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGAc                                                                                                      >  2:1468527/1‑139 (MQ=255)                aCCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACg                                                                                                     <  1:645586/139‑1 (MQ=255)                       cgcACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGAt                                                                                              <  1:1193265/139‑1 (MQ=255)                                    tGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTaa                                                                                 >  2:641592/1‑139 (MQ=255)                                      gCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAg                                                                                                                   >  2:1659493/1‑103 (MQ=255)                                      gCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAg                                                                                                                   <  1:1659493/103‑1 (MQ=255)                                      gCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACg                                                                               >  1:1579933/1‑139 (MQ=255)                                        gtcgtTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGcc                                                                                          >  2:1339990/1‑126 (MQ=255)                                        gtcgtTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGcc                                                                                          <  1:1339990/126‑1 (MQ=255)                                            ttAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAgg                                                                         <  1:875238/139‑1 (MQ=255)                                            ttAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAgg                                                                         <  1:395202/139‑1 (MQ=255)                                            ttAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAgg                                                                         <  1:1032384/139‑1 (MQ=255)                                                   tCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAac                                                                  >  2:378451/1‑139 (MQ=255)                                                                   ttCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGc                                                                                       <  2:1800784/102‑1 (MQ=255)                                                                   ttCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGc                                                                                       >  1:1800784/1‑102 (MQ=255)                                                                   ttCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTgaga                                                  >  2:280556/1‑139 (MQ=255)                                                                   ttCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTgaga                                                  >  1:935129/1‑139 (MQ=255)                                                                     cGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTa                                                                             <  2:1799759/110‑1 (MQ=255)                                                                     cGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTa                                                                             >  1:1799759/1‑110 (MQ=255)                                                                     cGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAAc                                                >  2:214454/1‑139 (MQ=255)                                                                     cGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAAc                                                <  1:641592/139‑1 (MQ=255)                                                                        tggtCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGc                                             >  1:766084/1‑139 (MQ=255)                                                                           tCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCacc                                          <  1:1468527/139‑1 (MQ=255)                                                                             ggCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGa                                                                    <  1:1525922/111‑1 (MQ=255)                                                                             ggCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGa                                                                    >  2:1525922/1‑111 (MQ=255)                                                                             ggCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCaccga                                        <  1:1628837/139‑1 (MQ=255)                                                                             ggCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCaccga                                        >  2:1724188/1‑139 (MQ=255)                                                                             ggCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCaccga                                        >  2:1683381/1‑139 (MQ=255)                                                                                       tgctgcGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGaa                              <  2:935129/139‑1 (MQ=255)                                                                                        gctgcGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAgg                                                                         >  1:270041/1‑95 (MQ=255)                                                                                        gctgcGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAgg                                                                         <  2:270041/95‑1 (MQ=255)                                                                                         ctgcGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACg                                                                                                     <  2:1178367/66‑1 (MQ=255)                                                                                         ctgcGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAAt                            >  1:796125/1‑139 (MQ=255)                                                                                         ctgaGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACg                                                                                                     >  1:1178367/1‑66 (MQ=255)                                                                                          tgcGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAAtt                           <  1:214454/139‑1 (MQ=255)                                                                                            cgggcggCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTAc                         <  1:280556/139‑1 (MQ=255)                                                                                              ggcggcAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTaaa                                                                           <  1:1034/87‑1 (MQ=255)                                                                                              ggcggcAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTaaa                                                                           >  2:1034/1‑87 (MQ=255)                                                                                                 ggcAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAgcggc                                                        <  2:1485061/103‑1 (MQ=255)                                                                                                 ggcAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAgcggc                                                        >  1:1485061/1‑103 (MQ=255)                                                                                                   cAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCa                                                                                                                 <  2:195986/44‑1 (MQ=255)                                                                                                   cAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCa                                                                                                                 >  1:195986/1‑44 (MQ=255)                                                                                                   cAAAACACATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATa                  >  1:1799005/1‑139 (MQ=255)                                                                                                         gCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAg            <  2:796125/139‑1 (MQ=255)                                                                                                         gCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAg            <  2:1579933/139‑1 (MQ=255)                                                                                                              aGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCt       <  1:1683381/139‑1 (MQ=255)                                                                                                                  tACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCaaa   >  1:758507/1‑139 (MQ=255)                                                                                                                   aCGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCACCTCaaaa  >  1:647243/1‑139 (MQ=255)                                                                                                                       |                                                                                                                                       GCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAA  >  NC_000913/4495214‑4495466

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