Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,507,424 G→C pseudogene (228/255 nt) insO → IS911B regulator fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,507,4240GC93.5% 83.0 / ‑4.8 31pseudogene (228/255 nt)insOIS911B regulator fragment
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (14/15);  minor base A (0/2);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CCAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTGTTCTCATTTGTTCCCCGCTGCCAAGGACTCTGAGGA  >  NC_000913/4507194‑4507670
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
ccAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGaa                                                                                                                                                                                                                                      <  1:84181/249‑1 (MQ=255)
  aaTGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAatat                                                                                                                                                                                                                                  >  2:87414/1‑251 (MQ=255)
              aaaTTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTAc                                                                                                                                                                                                                      >  2:80774/1‑251 (MQ=255)
                          agagTTTAAACGCGACTACGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCCGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATATGTAAGCTGAGGAAAAAGCTAAAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgat                                                                                                                                                                                                          <  2:89523‑M2/251‑15 (MQ=255)
                                   aCGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGAATTGACATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactc                                                                                                                                                                                                 <  1:171125‑M2/251‑24 (MQ=255)
                                    cgcgAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactcc                                                                                                                                                                                                <  1:208541‑M2/251‑25 (MQ=255)
                                                  aCTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcg                                                                                                                                                                                  >  2:129864‑M2/1‑213 (MQ=255)
                                                           tgACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacaggtagcgataatcgg                                                                                                                                                                          >  2:1896‑M2/1‑204 (MQ=255)
                                                           tgACAAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACACTGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACAACAGAACAAATCGAAATACGGAAGCTGAGGAAAAAGCTACAACGAATTGAAATGGAGAATGCAATATTAACAAGGCTaccgcgctattgatgacagactccctgaacagttatcgataatcggg                                                                                                                                                                         <  2:105581‑M2/251‑48 (MQ=255)
                                                             aCCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGACATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcggga                                                                                                                                                                        <  2:260631‑M2/250‑49 (MQ=255)
                                                             aCCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGACTGAAATATTAAAAAGGCTaccgcgatcttgatgtcagactccctgaacagttctcgataatcggga                                                                                                                                                                        <  2:7244‑M2/250‑49 (MQ=255)
                                                                         cGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTTCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgca                                                                                                                                                           >  2:18001‑M2/1‑190 (MQ=255)
                                                                               cAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGAAAGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgaattatcc                                                                                                                                                     <  1:135522‑M2/251‑68 (MQ=255)
                                                                                                                          aaGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagc                                                                                                          >  1:80774‑M2/1‑141 (MQ=255)
                                                                                                                                ggTCAAACAACTGCGTGATGAGAGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGCAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactcactgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggatcatcgcagcagctac                                                                                                       <  1:45779‑M2/248‑114 (MQ=255)
                                                                                                                                        aacTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAAACGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaa                                                                                            >  2:208541‑M2/1‑127 (MQ=255)
                                                                                                                                                 atgaGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcc                                                                                   >  1:260631‑M2/1‑118 (MQ=255)
                                                                                                                                                       gTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGACATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaa                                                                              <  1:113516‑M2/250‑139 (MQ=255)
                                                                                                                                                       gTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaa                                                                              <  1:41181‑M2/250‑139 (MQ=255)
                                                                                                                                                           gggCAAAACACCAAAAGCCTCTCCGATAAAACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGATACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgatactcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagaagctacagatactggaaaaaccgtcctgaagaacc                                                                          <  2:109730‑M2/250‑143 (MQ=255)
                                                                                                                                                                              tctcCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGCAAAAGCTACAGCGCATTGAAGTAGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactcactgaacagttctcgataatcgggcaactcagggcgcattatcctgtggtaacactctgccgtgtgatcggagtgcatcgcagcagctacagatactggaaaagccgtcctgaaaagccagacggcagacgggctgtat                                                      <  1:64543‑M2/251‑163 (MQ=255)
                                                                                                                                                                                   gATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAACGCTACCACGCATTGAAATGGAGAATGACATATTAACAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgc                                                 <  1:239261‑M2/251‑168 (MQ=255)
                                                                                                                                                                                            aGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgcagtcaggta                                        >  1:87414‑M2/1‑75 (MQ=255)
                                                                                                                                                                                            aGAACAAATCGAAATACGGAAGCAGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactgtaaaaaccgccctgaaaaaccagacggcagacgggctgtattacgcagtcaggta                                        >  2:122177‑M2/1‑75 (MQ=255)
                                                                                                                                                                                             gAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagct                                                                                                         <  1:46444‑M2/185‑112 (MQ=255)
                                                                                                                                                                                             gAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcagcgcagcagct                                                                                                         >  2:46444‑M2/1‑74 (MQ=255)
                                                                                                                                                                                                        aaTACGTAAGCGGAGGAAAAAGCGACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgctctcttgatgtcagactccctgaacagttctctataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatacttgaaaaaccgtcctgaaaagccatacggcagacgggccgtagtacgcagtctgttacttgagctacat                            >  2:171125‑M2/1‑63 (MQ=255)
                                                                                                                                                                                                                          aaaGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgttggcacactctgccatgtgttcggggttcatcgcagcagctacagatgctggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgcagtcaggtacttgagctacatggcatcagccacggttcg          >  1:89523‑M2/1‑45 (MQ=255)
                                                                                                                                                                                                                          aaaGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgcagtcaggtacttgagctacatggcatcagccacggttcg          >  2:2174‑M2/1‑45 (MQ=255)
                                                                                                                                                                                                                          aaaGCTACAACGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgcagtcaggtacttgagctacatggcatcagccacggttcg          <  1:2174‑M2/251‑207 (MQ=255)
                                                                                                                                                                                                                                   aCGCATTGAAATGGAGAATGAAATATTAAAAAGGCTaccgcgctcttgatgtcagactccctgaacagttctcgataatcgggaaactcagagcgcattatcctgtggtcacactctgccatgtgttcggggttcatcgcagcagctacagatactggaaaaaccgtcctgaaaaaccagacggcagacgggctgtattacgcagtcaggtacttgagctacatggcatcagccacggttcggccggagc  <  1:59693‑M2/250‑215 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
CCAATGAAAAAAAGAAATTTCAGCGCAGAGTTTAAACGCGAATCCGCTCAACTGGTTGTTGACCAGAAATACACGGTGGCAGATGCCGCGAAAGCTATGGATGTTGGCCTTTCCACAATGACAAGATGGGTCAAACAACTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACATGTCTTTTTAGGAAAGTACTGCCGAATTAGCCCATTTGTGTTCTCATTTGTTCCCCGCTGCCAAGGACTCTGAGGA  >  NC_000913/4507194‑4507670

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: