Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,496,168 G→T G14G (GGC→GGA ahr ← NADPH‑dependent aldehyde reductase Ahr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,496,1680GT90.0% 50.0 / NA 21G14G (GGC→GGAahrNADPH‑dependent aldehyde reductase Ahr
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (8/10);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCCAATCCCGACACGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGT  >  NC_000913/4495938‑4496404
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
catccAATCCCGTCACGCTGACCGACCTGCAATCCTTTATCCTTCGCGGCGCTCCCGATTGCCACCACTCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTACCGATCATCGACATATCGTAATGTCAGATCCCGCAGTAATCCACCTGCACTTCAACATCGTGTGTCCTCAGCTCACCTGGATCTTACTCATAATCTTCCAGTTCTCCGCCCGCTTCTCTTGCGGc                                                                                                                                                                                                                          <  1:224166/248‑1 (MQ=255)
         ccGACACGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCTCCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCGTTTGCGGCATAGCtttt                                                                                                                                                                                                                 <  2:36514/251‑1 (MQ=255)
                gCTGACCGACCTGCACACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTatcatc                                                                                                                                                                                                          <  2:28724/251‑1 (MQ=255)
                     ccgaccTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCCCCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACAt                                                                                                                                                                                                     >  2:285/1‑251 (MQ=255)
                              aaaCCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTCTCGCTCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTc                                                                                                                                                                                            >  2:114323/1‑251 (MQ=255)
                                                  gcTCCCGAGTGCCACCACGCGCCCAAGCACCTCAGGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGGTCTCCGCCCGCTTCTTTTGCGTCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTggg                                                                                                                                                                        <  2:262116/251‑1 (MQ=255)
                                                                 cacGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCAGTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGGGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCGGCTCTGGTGTTGTTTGGGAAGTGTAGAGCATgg                                                                                                                                                         <  1:103418/251‑1 (MQ=255)
                                                                        ccAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGGGCATGGCAGGGCg                                                                                                                                                  >  2:3396/1‑251 (MQ=255)
                                                                               ccTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGTCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCGCCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCGCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGTGTGcc                                                                                                                                           <  2:190662/251‑1 (MQ=255)
                                                                                                                         cGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACCTCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTg                                                                                                  >  2:221853/1‑250 (MQ=255)
                                                                                                                            tATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGaaa                                                                                               <  1:169940/250‑1 (MQ=255)
                                                                                                                                      cGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAgg                                                                                                                                                     >  2:188711/1‑186 (MQ=255)
                                                                                                                                      cGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAgg                                                                                                                                                     <  1:188711/186‑1 (MQ=255)
                                                                                                                                         cAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTCCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATGTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGTTTACACATCGCTGTATGCAATGGTGAAAATTTCAGCACTTAg                                                                                 <  1:100990/251‑1 (MQ=255)
                                                                                                                                            aTCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCg                                                                               >  1:127329/1‑250 (MQ=255)
                                                                                                                                               ggAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGt                                                                           <  1:285/251‑1 (MQ=255)
                                                                                                                                                        gATCCCGCAGTAATCCACCTGCACTCCAACATCGTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACGTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTGTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTTGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAg                                                                  <  2:269005/251‑1 (MQ=255)
                                                                                                                                                                                ttCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATACCCTTCCCGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTCTCATCGACATTTTTTGGTCCTTCTCTGGTGTTGATTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGCAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCCTGGTGGCGTGCg                                          >  1:246118/1‑251 (MQ=255)
                                                                                                                                                                                  caccATCTTGGGGCCTCGGCGCGCCGGGATCGTACTCATAAGATTCCAGTTCCCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGAc                                        <  2:229671/248‑1 (MQ=255)
                                                                                                                                                                                                             gATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCAt             >  1:253430/1‑251 (MQ=255)
                                                                                                                                                                                                                         aaaCTTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGt  <  2:253430/250‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
CAGCCAATCCCGACACGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACAGATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGT  >  NC_000913/4495938‑4496404

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