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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913271,9010GA83.3% 23.1 / ‑4.4 13pseudogene (138/427 nt)insNIS911A regulator fragment
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (6/4);  total (8/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTC  >  NC_000913/271673‑272042
                                                                                                                                                                                                                                    |                                                                                                                                             
tctccgataacaccagaacaaatcgaaatacgtaagctgaggcaaaagctacagcgcattgaagtagagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCACTGAACAGTTCTCGATAATCGGGCAACTCAGGGCGCATTATCCTGTGGTAACACTCTGCCGTGTGATCGGAGTGCATCGCAGCAGCTACAGATACTGGAAAAGCCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTAt                                                                                                                         <  1:64543‑M1/162‑1 (MQ=255)
     gataacaccagaacaaatcgaaatacgtaagctgaggaaaacgctaccacgcattgaaatggagaatgacatattaacaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGc                                                                                                                    <  1:239261‑M1/167‑1 (MQ=255)
              agaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTa                                                                                                           >  1:87414‑M1/76‑251 (MQ=255)
              agaacaaatcgaaatacggaagcagaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGTAAAAACCGCCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTa                                                                                                           >  2:122177‑M1/76‑251 (MQ=255)
                          aatacgtaagcggaggaaaaagcgacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCTCTCTTGATGTCAGACTCCCTGAACAGTTCTCTATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTTGAAAAACCGTCCTGAAAAGCCATACGGCAGACGGGCCGTAGTACGCAGTCTGttacttgagctacat                                                                                               >  2:171125‑M1/64‑244 (MQ=255)
                                            aaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTTGGCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATGCTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGccacggttcg                                                                             >  1:89523‑M1/46‑247 (MQ=255)
                                            aaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGccacggttcg                                                                             <  1:2174‑M1/206‑5 (MQ=255)
                                            aaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGccacggttcg                                                                             >  2:2174‑M1/46‑247 (MQ=255)
                                                     acgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGc                                                                     <  1:59693‑M1/214‑1 (MQ=255)
                                                          ttgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCaaga                                                                 <  1:41345‑M1/218‑1 (MQ=255)
                                                                                aaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCCACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAacccg                                          >  1:109730‑M1/10‑249 (MQ=255)
                                                                                          ccGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCCCACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTACTGAAAACCCAGCCGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAGCTACATGGCATCAGCAACGGTTCGGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTAcc                               >  1:49575/1‑251 (MQ=18)
                                                                                                                        gTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGAAAAACAAGACGGCAGACGGGCTGTATTACGCAGTCAGGGACTTGAGCTACATGACATCAGCAACGGATCGGCCGAAGCAAGAAGCATCGCCACAATGGCAACCCGGAGAGGCTACCAGCTGGGACGCTGGCTTGCTGGCAgggtc  >  1:225738/1‑246 (MQ=17)
                                                                                                                                                                                                                                    |                                                                                                                                             
GCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTC  >  NC_000913/271673‑272042

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