Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,496,351 C→T 23.8% intergenic (‑142/‑54) ahr ← / → leuX aldehyde reductase, NADPH‑dependent/tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,496,3510CT23.8% 24.9 / 6.1 17intergenic (‑142/‑54)ahr/leuXaldehyde reductase, NADPH‑dependent/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base C (6/7);  new base T (2/2);  total (8/9)
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.65e-01

CATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCA  >  NC_000913/4496119‑4496607
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                
cctCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTg                                                                                                                                                                                                                                                       <  1:75761/242‑1 (MQ=255)
cATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTg                                                                                                                                                                                                                                                       >  2:75761/1‑244 (MQ=255)
        ggCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGTATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTcc                                                                                                                                                                                                                            >  2:88492/1‑263 (MQ=255)
        ggCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGTATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTcc                                                                                                                                                                                                                            <  1:88492/263‑1 (MQ=255)
                 agacagggATCGTACTCATAAACTTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGcc                                                                                                                                                                                                          <  1:224256/267‑1 (MQ=255)
                                     aaCTTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACATCTGCACGGCCGAAGTGGCTGAATccctgcct                                                                                                                                                                                     >  2:224256/1‑266 (MQ=255)
                                             gTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATc                                                                                                                                                                                            >  1:49150/1‑258 (MQ=255)
                                                                                       aCATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAAc                                                                                                                                                                                                                                   <  1:73661/177‑1 (MQ=255)
                                                                                        cATTTTTTGGTCTTCTTGTGGTTTTGTTTGGGAGGTGTAGAGCATGGCAGGGGGGGGTGCTTGGAGTGTGACAAGGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCTACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGcc                                                                                                                                  <  1:89517/273‑1 (MQ=255)
                                                                                                                                                  gCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGTATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAACACCCAAATGAGGTTATTTTTTACTCCTGaaaa                                                                        >  2:212217/1‑272 (MQ=255)
                                                                                                                                                                            tCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCta                                              <  1:281643/273‑1 (MQ=255)
                                                                                                                                                                                tGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATAtt                                          >  1:162781/1‑273 (MQ=255)
                                                                                                                                                                                 gTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATAtt                                          >  1:26752/1‑272 (MQ=255)
                                                                                                                                                                                                   ttCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAg                                                                                                                    <  2:202207/180‑1 (MQ=255)
                                                                                                                                                                                                   ttCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAg                                                                                                                    >  1:202207/1‑180 (MQ=255)
                                                                                                                                                                                                                      cGAGCAAGCTGGCGCTTGTATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATc       <  1:212217/270‑1 (MQ=255)
                                                                                                                                                                                                                        aGCATGCTTGCGCTGTCATGGTTGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTTCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTTCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCa  <  1:260545/273‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                
CATCTTGTGGCCTCAGCTCACCGGGATCGTACTCATAAACTTCCAGTTCGCCGCCCGCTTCTTTTGCGGCATAGCTTTTTATCATCGACATTTTTTGGTCCTTCTCTGGTGTTGTTTGGGAAGTGTAGAGCATGGCAGGGCGGGGTGCCTGGAGTGTGACAAAGGTTACACATCGCTGTATGCAATGCTGAAAATTTCAGCACTTAGCGAGGTGCGAGCAAGCTGGCGCTTGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAAGTATGTAAATAGACCTCAACTGAGGTCTTTTTTTATGCCTGAAATCCAGTGTTTATCTATCTTTCCCGCTATATTAACTCTCTCAAGGTCAACCGATATCAACGTACATCTACCA  >  NC_000913/4496119‑4496607

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: