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,1680GT93.9% 70.1 / ‑3.1 33G14G (GGC→GGA) ahrNADPH‑dependent aldehyde reductase Ahr
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (13/18);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

ACGCGCCGTCCAGCCAATCCCGACACGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACA‑GATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCC  >  NC_000913/4495928‑4496407
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
aCGCGCCGTCCAGCCAATCCCGACACGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCGCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACA‑GATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGctt                                                                                                                                                                                                                                       <  1:358480/251‑1 (MQ=255)
                      acacGCTGACCGACCTGCAAACCTTTATCCTACGCGGCGCTCCCGAGCGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACA‑GATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTat                                                                                                                                                                                                                 >  1:101019/1‑251 (MQ=255)
                      aaacGCAGACCGACCTGCAAACCTTTATCCTGCACGGCGCTCCCGAGCGCCACCACGCGCACAATCACCTCATGCCCGGCAACCAGCGGATATTGAGAAAATCCCCATTCGTTATCGATCATCGACA‑GATCGGAATGACAGATCCCGCAGTAATACACTTTCACTTCAACCTCTTGTGGCCGCAGCTCACCAGGGTCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATTGCTtattat                                                                                                                                                                                                                 >  1:101028/3‑251 (MQ=255)
                         cGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCACCGACA‑GATCGGAATGGCAGATCCCGCAGGAATCCACCTGCACTTCAACATCTTGTTGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTatcat                                                                                                                                                                                                              <  1:240989/251‑1 (MQ=255)
                                   cccacAAACCTTTACTCGGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCCCCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCTTTATCGATCACCGACA‑GATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCTTACTCATAAACTTCCCGTTCTCCGCCCGCTTCTTATGCGGCATAGCTTTTTATCATCGACAttttt                                                                                                                                                                                                    <  2:294011/247‑1 (MQ=255)
                                        aaaCCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACA‑GATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTc                                                                                                                                                                                               <  1:258769/251‑1 (MQ=255)
                                              ttATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTAGCGATCATCGACA‑GATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACAGCTTGTGGCCTCAGCTCACCGGGATCGGACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTtc                                                                                                                                                                                           <  1:60960/249‑1 (MQ=255)
                                               tATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACA‑GATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTctct                                                                                                                                                                                        >  1:337021/1‑251 (MQ=255)
                                                  ggttaaaggcgCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCTGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACC‑GATCGGAATGGCAGATCCCGCAGTAATCCACATGCACTTCAACATCTTGTGGCCTCAGCTAACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATCACTTTTTATCATCGACATTTTTTGGTCCTTCTCTggg                                                                                                                                                                                     >  2:302620/8‑250 (MQ=255)
                                                        cgtcgcTCCCGAGGGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTAGCGATCATCGACA‑GATCGGAATGGCAGATCCCGCAGTAATCCAACTGCACTTCAACATCTTATGGCCTCAGCTCACCGGGATCGTACTAATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCAACATTTTTTGGTCCTTCTCTGTGgttgtt                                                                                                                                                                               >  2:335676/4‑251 (MQ=255)
                                                                                     aataCCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACA‑GATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCgggg                                                                                                                                                  >  1:125026/4‑251 (MQ=255)
                                                                                       cacctaatgccaaacAACCATCGGATATTGTTAAAATCCCCTTTCGTTATCGATCATCGACA‑GTTCGGAATGGCATATCACGCAGTAGTCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTTCAGTTCTCAGCCCGCTTCATTTGCGGCATAGCTTTTTATCCTCGACCTTATTTGGTCCTTCTCTGGAGGTGTTTGGGTAGTGTAGAGCCAGGTAGGGCGGGGtt                                                                                                                                                >  2:24398/15‑250 (MQ=255)
                                                                                         ccTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTAACGATCATCGACA‑GATCGGAATGGCAGATCCCGCCGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCACGCTTCTTTTGCGGCATCGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGATGcc                                                                                                                                              >  1:221772/1‑251 (MQ=255)
                                                                                           tCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACA‑GATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTg                                                                                                                                            >  1:261572/1‑251 (MQ=255)
                                                                                                         cagtggATATTGTTCAAATCCCCATTCGTTATCGATCATCGACA‑GATCGGAAGGGCAGATCTCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCTTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCGCCGCCCGCTGCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTTGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGTATTGTGACAAAgg                                                                                                                              <  2:191797/247‑1 (MQ=255)
                                                                                                         cagcggctgttgataaaaTCCCCATTCGTATTCGACAATCGACAAGA‑CGGAATGCCAGTCCCCGCAGTCATCCACCTGCACTCCAACATCTTGTGTCCTCAGCTCACCTGTAGCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCGCTGGTGGTGTTTGGGAAGCGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAgg                                                                                                                              <  2:128865/237‑1 (MQ=255)
                                                                                                           gCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACA‑GATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGtt                                                                                                                            <  1:27700/251‑1 (MQ=255)
                                                                                                                       aaaaTCCCCATGCGTTATCGACCTTCGACA‑GATCGGAAGGGCAGATCCCGCAGTAATCCACCTGCACTGCACCATCTTGTGGCCTCAGCTCACCGGGAGCGTACTCATAACCTTCCAGTTCTCGGCCCGCTTCTGTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCAGGGGTTGTTTGGGAAGTGTAGAGCATGGCAGTGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGt                                                                                                                <  1:145592/251‑1 (MQ=255)
                                                                                                                                              atcGACA‑GATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCa                                                                                         <  1:302620/251‑1 (MQ=255)
                                                                                                                                                                     tCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCt                                                                   >  1:240920/1‑251 (MQ=255)
                                                                                                                                                                         gCAGTACTCCCCCTGCACTTCAACATCCTGCTGCCTCTGCTCACCGTGATCCTACTCACAAACTTCCATTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTTGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGTTGCGAGCAAGCTGgcg                                                               <  2:261572/251‑1 (MQ=255)
                                                                                                                                                                              aaTCCACCTGCACTTCAACAACTTGTGGCCTCAGCTCTCCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGc                                                          <  1:82426/251‑1 (MQ=255)
                                                                                                                                                                               aTCCACCTGCACTTCCACATCTTGTGGCCTCAGCTCACCGGGATCGAAATAATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCCTAGATTTTTATAATCGACATTTTTTGGTCCTTCTCTGGAGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACACAGGTTAAACATCGCTGTATGCAATGCTGAAATTTTCAACACTTAGCGCGGTGTTAGTAAGCTGGCGCTTGc                                                          >  2:236537/1‑250 (MQ=255)
                                                                                                                                                                                tCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCAACGGGATCGTACTCATACACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTAAAAATCGCTGTATGCAGAGCTGGAACTTTCGGCACTTAGAGAGTTGCGAGCAAGCTGTCGCTTGcct                                                        >  2:41499/1‑249 (MQ=255)
                                                                                                                                                                                            tCAACATCTTGTGGCCTCAGCTCCCCGGGTTCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCg                                             <  2:63061/250‑1 (MQ=255)
                                                                                                                                                                                                          ccTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCcac                              >  1:104439/1‑251 (MQ=255)
                                                                                                                                                                                                                  cACCGGGATCGTACCAATGAACTACCTGTTATCCGGCCGCTTCTTTTGCGGCTAAGCTTTTTATCATCTACATTTTTTGGGCCGTCGCTGGTGTTGTGTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGTTTACACAGCGCTGTATGCAATGCTGCAAATGTCAGCACTTAGCGAGGTACGAGCAAGCTGTCGCTTTCATGGTGGTGGGCGACAGGTATAATCCACTACGtt                        <  2:245590/249‑1 (MQ=255)
                                                                                                                                                                                                                         aTCGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATa               <  1:378631/251‑1 (MQ=255)
                                                                                                                                                                                                                           cGTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATAcc             >  1:373277/1‑251 (MQ=255)
                                                                                                                                                                                                                            gTACTCATAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACct            <  1:373709/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aTAAACTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAg      <  1:208138/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aaCTTCCATTTCTCCTCCCTCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTTGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGTTGCCTGGAGTTTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGc   <  1:77094/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aCTTCCAGTTCTCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATCCCTCTTCAGTGcc  <  1:212236/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
ACGCGCCGTCCAGCCAATCCCGACACGCTGACCGACCTGCAAACCTTTATCCTGCGCGGCGCTCCCGAGTGCCACCACGCGCCCAATCACCTCATGCCCGGCAACCAGCGGATATTGTGAAAATCCCCATTCGTTATCGATCATCGACA‑GATCGGAATGGCAGATCCCGCAGTAATCCACCTGCACTTCAACATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCC  >  NC_000913/4495928‑4496407

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