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_0009134,495,3320AG97.4% 135.0 / NA 39L293P (CTG→CCG) ahrbroad specificity NADPH‑dependent aldehyde reductase, Zn‑containing
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (17/21);  minor base T (1/0);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCC  >  NC_000913/4495214‑4495460
                                                                                                                      |                                                                                                                                
gcgGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTa                                                                                                              <  2:1278182/139‑1 (MQ=255)
 cgGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAg                                                                                                             <  2:638528/139‑1 (MQ=255)
  ggTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGc                                                                                                                       >  1:1006566/1‑128 (MQ=255)
  ggTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGc                                                                                                                       <  2:1006566/128‑1 (MQ=255)
   gTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCa                                                                                                           <  2:1070311/139‑1 (MQ=255)
          gcCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCag                                                                                                    >  1:1101537/1‑139 (MQ=255)
          gcCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCag                                                                                                    >  2:331022/1‑139 (MQ=255)
           ccTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCaga                                                                                                   <  2:1101537/139‑1 (MQ=255)
           ccTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCaga                                                                                                   <  2:1224020/139‑1 (MQ=255)
                      cgcACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGAt                                                                                        >  2:615448/1‑139 (MQ=255)
                        cacaTGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCg                                                                                      >  2:206904/1‑139 (MQ=255)
                           aTGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGTCGGTAGCAGAACCAGAGACGCTGCGATCGccc                                                                                   >  1:522107/1‑139 (MQ=255)
                               tGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCaa                                                                               <  1:331022/139‑1 (MQ=255)
                                     gCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACg                                                                         >  1:1191906/1‑139 (MQ=255)
                                       gtcgtTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTa                                                                       >  2:184969/1‑139 (MQ=255)
                                         cgtTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGGCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTaaa                                                                     >  2:958266/1‑139 (MQ=255)
                                            tAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGc                                                                  <  1:615448/139‑1 (MQ=255)
                                                                  ttCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAg                                                                                                             <  2:873085/74‑1 (MQ=255)
                                                                  ttCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAg                                                                                                             >  1:873085/1‑74 (MQ=255)
                                                                  ttCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAgg                                                                   >  2:1146247/1‑116 (MQ=255)
                                                                  ttCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAgg                                                                   <  1:1146247/116‑1 (MQ=255)
                                                                          tCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCacc                                    <  2:634198/139‑1 (MQ=255)
                                                                          tCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCacc                                    >  1:1309443/1‑139 (MQ=255)
                                                                            ggCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGTCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCaccga                                  <  2:522107/139‑1 (MQ=255)
                                                                            ggCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGAt                                                                                        <  1:657265/85‑1 (MQ=255)
                                                                            ggCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGAt                                                                                        >  2:657265/1‑85 (MQ=255)
                                                                            ggCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGCCTCATAAGGCGTGCCGGTAGCAGAACCAGAGTCGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCaccga                                  <  1:889881/139‑1 (MQ=255)
                                                                                 aaCCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTa                             <  1:206904/139‑1 (MQ=255)
                                                                                         tgcGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGGCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAAtt                     <  1:958266/139‑1 (MQ=255)
                                                                                         tgcGGGCGGCAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCaccga                                  <  2:794256/126‑1 (MQ=255)
                                                                                         tgagagcggcAAAACGCATCAGCTTACGCTGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCaccga                                  >  1:794256/6‑126 (MQ=255)
                                                                                             ggcggcAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGc                                                                                                            <  2:1132388/48‑1 (MQ=255)
                                                                                             ggcggcAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGc                                                                                                            >  1:1132388/1‑48 (MQ=255)
                                                                                             ggcggcAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGcc                                                                 >  1:310189/1‑91 (MQ=255)
                                                                                             ggcggcAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGcc                                                                 <  2:310189/91‑1 (MQ=255)
                                                                                             ggcggcAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTAc                   <  2:654727/137‑1 (MQ=255)
                                                                                             ggcggcAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTAc                   >  1:654727/1‑137 (MQ=255)
                                                                                             ggcggcAAAACGCATCAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCg                 >  2:1331629/1‑139 (MQ=255)
                                                                                                            cAGCTTACGCGGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGcc  <  2:1309443/139‑1 (MQ=255)
                                                                                                                      |                                                                                                                                
GCGGTAACGCGCCTTACCGTCGCGCACATGCTGGATGGCGTCGTTAATTTTCGACATCGGGAACAGTTCGGTGGTCGGCGCAACCTTGCTGCGGGCGGCAAAACGCATCAGCTTACGCAGCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCC  >  NC_000913/4495214‑4495460

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