Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,294,583 G→A 8.1% pseudogene (410/567 nt) insZ ← pseudogene, transposase homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,294,5830GA 8.1% 146.4 / 3.3 49pseudogene (410/567 nt)insZpseudogene, transposase homolog
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base A (2/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

AGGAGCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGTTCAGCGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCGGGAT  >  NC_000913/1294435‑1294728
                                                                                                                                                    |                                                                                                                                                 
aggagCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCg                                                                                                                                                   >  4:335572/1‑149 (MQ=255)
     cGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGttt                                                                                                                                              >  7:281369/1‑149 (MQ=255)
             tGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGctc                                                                                                                                      <  2:91117/149‑1 (MQ=255)
               tttCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACGCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCGTGGTGACGCCGTGTTTCCCAGCTCTa                                                                                                                                    <  7:93120/149‑1 (MQ=255)
               tttCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTa                                                                                                                                    >  5:27681/1‑149 (MQ=255)
                 tCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTAt                                                                                                                                   <  1:343407/148‑1 (MQ=255)
                  cacaGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCATGTTTCCCAGCTCTATCa                                                                                                                                 >  6:291805/1‑149 (MQ=255)
                    caGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCAtt                                                                                                                               >  4:106262/1‑149 (MQ=255)
                     aGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCATGTTTCCCAGCTCTATCAttt                                                                                                                              <  5:291805/149‑1 (MQ=255)
                         tAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGAGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGa                                                                                                                          >  3:79853/1‑149 (MQ=255)
                          aaCGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGAt                                                                                                                         >  7:307631/1‑149 (MQ=255)
                          aaCGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGAt                                                                                                                         >  3:98595/1‑149 (MQ=255)
                                tttttGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGata                                                                                                                   >  5:304447/1‑149 (MQ=255)
                                      atatCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATAttctt                                                                                                              <  6:353612/148‑1 (MQ=255)
                                       tatCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTcc                                                                                                            <  8:281369/149‑1 (MQ=255)
                                          cTTAACCGACCTCGGACTTGATTGTCGGGTGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTTGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCAc                                                                                                          <  4:79853/148‑1 (MQ=255)
                                           ttAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACg                                                                                                         >  1:143095/1‑148 (MQ=255)
                                                cGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACgcaggca                                                                                                   <  6:27681/149‑1 (MQ=255)
                                                   ccTCGGCCTTGATTGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTa                                                                                                <  5:16503/149‑1 (MQ=255)
                                                    cTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAg                                                                                               >  1:64417/1‑149 (MQ=255)
                                                        gCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAg                                                                                           <  4:98595/149‑1 (MQ=255)
                                                          cTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGAAGGTCCCTGGTGACGCCATGTTTCACATCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGACAGTg                                                                                         >  7:326028/1‑149 (MQ=255)
                                                           ttGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTgg                                                                                        >  2:167517/1‑149 (MQ=255)
                                                             gATGGTCGGGGGCGTTTTGTTATGAACACCAGTTACAGAGCAGCGCGAAGATGCTCCAGTCGTTTCGGGATGGTGCCTGGTGACGCCATGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCAAGCAGGCAGTAGCCAGTg                                                                                         <  8:326028/146‑1 (MQ=255)
                                                             gATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTgg                                                                                        <  3:335572/147‑1 (MQ=255)
                                                                 gTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTa                                                                                  >  3:19595/1‑149 (MQ=255)
                                                                       ggCGTTTTGTTATGAACAACACTTCGAGAGCACCCCGAAGATGCTCCAGTGGTTTCGGGATGGTCCCTGGTCACGACGTGTTTCGCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACATccc                                                                            <  8:354360/149‑1 (MQ=255)
                                                                         cGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCtt                                                                          <  1:286547/149‑1 (MQ=255)
                                                                           ttttGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCAtt                                                                                                                               <  2:322890/94‑1 (MQ=255)
                                                                           ttttGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCAtt                                                                                                                               >  1:322890/1‑94 (MQ=255)
                                                                           ttttGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCt                                                                        <  3:106262/149‑1 (MQ=255)
                                                                                ttATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGtt                                                                    >  1:246405/1‑148 (MQ=255)
                                                                                       caccacTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGTTCAGCGt                                                              >  8:197367/1‑147 (MQ=255)
                                                                                       caccacTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGTTCAGCGTca                                                            <  8:307631/149‑1 (MQ=255)
                                                                                          cacTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGCTGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTGCGGGTTCAGCGTcagcc                                                         >  8:358788/1‑148 (MQ=255)
                                                                                             ttCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCCCGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGTTCAGCGTcagcag                                                        >  5:55500/1‑147 (MQ=255)
                                                                                                cAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGTTCAGCGTCAGCAGCAGGt                                                   >  2:298063/1‑149 (MQ=255)
                                                                                                        cccGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGTTCAGCGTCAGCAGTAGGTCTTCGCTg                                           >  5:140572/1‑149 (MQ=255)
                                                                                                        cccGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGTTCAGCGTCAGCAGCAGGTCTTCGCTg                                           <  4:19595/149‑1 (MQ=255)
                                                                                                         ccGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGTTCAGCGTCAGCAGCAGGTCTTCGCTGt                                          <  2:362127/149‑1 (MQ=255)
                                                                                                                   tCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGTTCAGCGTCAGCAGCAGGTCTTCGCTGTAAAAAAGtt                                 >  7:33383/1‑148 (MQ=255)
                                                                                                                     cAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGa                                                                                                                          <  2:187936/57‑1 (MQ=255)
                                                                                                                     cAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGa                                                                                                                          >  1:187936/1‑57 (MQ=255)
                                                                                                                      aGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGTTCAGCGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTAt                              >  4:314232/1‑148 (MQ=255)
                                                                                                                                     gTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGTTCAGCGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATaga              <  1:167517/149‑1 (MQ=255)
                                                                                                                                     gTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGTTCAGCGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATaga              <  3:314232/149‑1 (MQ=255)
                                                                                                                                     gTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGTTCAGCGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATaga              <  2:143095/149‑1 (MQ=255)
                                                                                                                                       cccTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGTTCAGCGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAGAGt            >  8:67570/1‑149 (MQ=255)
                                                                                                                                            gTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGTTCAGCGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATc        <  2:64417/148‑1 (MQ=255)
                                                                                                                                                 gCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGTTCAGCGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCGGGAt  <  7:197367/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AGGAGCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGCACCCCGAAGATGCTCCAGTCGTTTCGGGATGGTCCCTGGTGACGCCGTGTTTCCCAGCTCTATCATTTCTGATGCGATATTCTTCCACGCAGGCAGTAGCCAGTGGCGGTTACAACCCTTCTGGTTCAGCGTCAGCAGCAGGTCTTCGCTGTAAAAAAGTTTATCAAACAACGTAATAGAGTTATCCGGGAT  >  NC_000913/1294435‑1294728

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑