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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913271,8170GA95.5% 61.8 / ‑6.0 22pseudogene (54/427 nt)insNIS911A regulator fragment
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (12/9);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAACAAATGAGAATACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCAATATACTCAACATGAACTAGATCTGGTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTC  >  NC_000913/271572‑272042
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 
cggtggcagatgccgcgaaagctatggatgttggcctttccacaatgacaagatgggtcaaacaacttcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGagcgca                                                                                                                                                                                                                              >  2:18001‑M1/191‑250 (MQ=255)
                                                 aagatgggtcaaacaactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCagcagc                                                                                                                                                                             >  1:80774‑M1/142‑251 (MQ=255)
                                                       ggtcaaacaactgcgtgatgagagtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgcaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCACTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGATCATCGCAGCAGCTAc                                                                                                                                                                          <  1:45779‑M1/113‑1 (MQ=255)
                                                               aactgcgtgatgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaaacgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGaa                                                                                                                                                               >  2:208541‑M1/128‑251 (MQ=255)
                                                                        atgagcgtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTcc                                                                                                                                                      >  1:260631‑M1/119‑251 (MQ=255)
                                                                              gtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgacatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGaaa                                                                                                                                                 <  1:113516‑M1/138‑1 (MQ=255)
                                                                              gtcagggcaaaacaccaaaagcctctccgataacaccagaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAGATACTGGAAAAACCGTCCTGaaa                                                                                                                                                 <  1:41181‑M1/138‑1 (MQ=255)
                                                                                  gggcaaaacaccaaaagcctctccgataaaaccagaacaaatcgaaatacgtaagctgaggaaaaagatacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATACTCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGAAGCTACAGATACTGGAAAAACCGtcctgaagaacc                                                                                                                                             <  2:109730‑M1/142‑6 (MQ=255)
                                                                                                     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)
                                                                                                                    gaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCt                                                                                                                                                                            <  1:46444‑M1/111‑1 (MQ=255)
                                                                                                                    gaacaaatcgaaatacgtaagctgaggaaaaagctacaacgcattgaaatggagaatgaaatattaaaaaggctACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGAGCGCATTATCCTGTGGTCACACTCTGCCATGTGTTCGGGGTTCAGCGCAGCAGCt                                                                                                                                                                            >  2:46444‑M1/75‑185 (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)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 
GGAACAAATGAGAATACAAATGGGCTAATTCGGCAGTACTTTCCTAAAAAGACATGTCTTGCCCAATATACTCAACATGAACTAGATCTGGTTGCTGCTCAGCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAACCGCGCTCTTGATGTCAGACTCCCTGAACAGTTCTCGATAATCGGGAAACTCAGGGCGCGTTATCCTGTGGCCACTCTCTGCCATGTGTTCGGGGTTCATCGCAGCAGCTACAAATACTGGAAAAACCGTCCTGAAAAGCCAGACGGCAGACGGGCTGTATTACGCAGTCAGGTACTTGAACTGCATGGCATCAGCCACGGCTCTGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTC  >  NC_000913/271572‑272042

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