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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913272,0020AG88.9% 16.2 / ‑1.2 18pseudogene (239/427 nt)insNIS911A regulator fragment
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (12/4);  total (12/6)
Fisher's exact test for biased strand distribution p-value = 9.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

TCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCA  >  NC_000913/271756‑272215
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                     
aaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCggaga                                                                                                                                                                                                                   >  1:291973‑M1/7‑251 (MQ=255)
     tACCGCGCTCTTGATGTCAGACTCCATGAACAGTTCTCGATAATCGGGAAAATCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGGTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAAACAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTAAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGGAACCCGGAgaggct                                                                                                                                                                                                               <  1:351208‑M1/249‑3 (MQ=255)
       ccGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGACACTCAGAGCGAATTATCCTGTAGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCTTGAAAAAACAGACGGCAGACGGGCTGTATTCCGCAGTAAGGTACTTGCACTACATGGCATAAGCCACGGTTTGGCCGGAGCAAGAAGCATCGCCACACGGGCAACCCagagcggcttcc                                                                                                                                                                                                            >  2:60439/1‑243 (MQ=255)
        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)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                     
TCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGTAACGGAACCAAATCAGGTGTGGTGCGGTGATGTGACCTATAGTGTGCCCGGAGTTCAGGGCGGGCA  >  NC_000913/271756‑272215

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