Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 272,009 A→C pseudogene (246/427 nt) insN → IS911A regulator fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913272,0090AC88.2% 13.8 / ‑0.4 17pseudogene (246/427 nt)insNIS911A regulator fragment
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (12/3);  total (12/5)
Fisher's exact test for biased strand distribution p-value = 7.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

CGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCA  >  NC_000913/271764‑272215
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                              
cgcgCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTAcc                                                                                                                                                                                                            >  1:156408/1‑250 (MQ=18)
    ctctTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATg                                                                                                                                                                                                       >  2:79961/1‑251 (MQ=18)
         gATGTCAGACTCCCTGAACAGTGCCCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATgggacg                                                                                                                                                                                                  <  1:427201/251‑4 (MQ=17)
               aGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCt                                                                                                                                                                                            >  1:66770/1‑251 (MQ=17)
                          aCAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAggc                                                                                                                                                                                 >  1:495618/1‑249 (MQ=14)
                              ttCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCAACAATGGCACACCGGAGAGGCTACCAGATGAGACGCTTGCTTACTGGCAggcaca                                                                                                                                                                              >  2:131604/1‑245 (MQ=14)
                                ctcGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATg                                                                                                                                                                            >  1:425928/1‑250 (MQ=14)
                                           ggAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCGGGCTCATGAAAGAGCTgggg                                                                                                                                                                >  2:280533/1‑251 (MQ=11)
                                             aaaCTCAGAGCTCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGAAAGAGCTGGGGCt                                                                                                                                                              <  2:156408/251‑1 (MQ=11)
                                              aaCTCAGAGCGAATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAACACCGTACTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGAATCAGCCACGGTTCGGCCGGAGAAAGAAGCATCGCCACAATGGCAACCCGGAGCGGCTACCAGATGGGACGCTGGCTTGCTGGACGGCCCATGAAAGAGCTgggggtg                                                                                                                                                             >  2:459821/1‑248 (MQ=11)
                                               aCTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCGCATGAAAGAGCTGGGGCTgg                                                                                                                                                            >  1:143196/1‑251 (MQ=11)
                                                   agagcgcATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGAAAGAGCTGGGGCTGGTc                                                                                                                                                          >  1:481864/4‑249 (MQ=11)
                                                     agcgcATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGCGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGAAAGAGCTGGGGCTGGACAgct                                                                                                                                                      >  1:177892/2‑249 (MQ=11)
                                                      gcgcATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCATTCAGGTACTTGAGCTACATGGCAGCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCAGGCTCATGAAAGAGCTGGGGCTGGTCAgctg                                                                                                                                                     <  2:216174/251‑4 (MQ=11)
                                                           ttATCCTGTGGTCACACTCTGCCATGTGTTAGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGCAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCCCAATGGCAACCCGGAGAGGCTACCAGATGGGACGCTGGCTTGCTGGCCGGCTAATGCAAGAGCTGGGGCTGGtaagctggcagc                                                                                                                                                >  2:452652/1‑240 (MQ=11)
                                                                    ggCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAAGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCGACTc                                                                                                                                       <  2:251442/251‑1 (MQ=255)
                                                                                                                                                                                                            gCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCAAGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCa  <  1:48899/248‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                              
CGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCA  >  NC_000913/271764‑272215

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