Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 280,171 G→A pseudogene (192/249 nt) insX ← pseudogene, IS3 family transposase family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913280,1710GA97.1% 125.0 / ‑4.2 34pseudogene (192/249 nt)insXpseudogene, IS3 family transposase family
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (11/22);  total (12/22)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

GCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGGCGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTT  >  NC_000913/280023‑280317
                                                                                                                                                    |                                                                                                                                                  
gCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGacgacg                                                                                                                                                  <  1:520768/151‑1 (MQ=255)
    aacaacAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGTAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGAtc                                                                                                                                              >  1:219430/1‑151 (MQ=255)
     acaacaGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGAtct                                                                                                                                             <  2:482967/151‑1 (MQ=255)
      caacaGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGAtctc                                                                                                                                            >  2:252565/1‑151 (MQ=255)
       aacaGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCg                                                                                                                                           <  1:541568/151‑1 (MQ=255)
       aacaGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCg                                                                                                                                           <  1:313129/151‑1 (MQ=255)
        acaGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCgg                                                                                                                                          >  2:273260/1‑151 (MQ=255)
         caGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGc                                                                                                                                         <  1:240170/151‑1 (MQ=255)
              cGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTg                                                                                                                                      >  2:362370/1‑149 (MQ=255)
                agaAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGc                                                                                                                                   >  1:121209/1‑150 (MQ=255)
                 gaAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCAt                                                                                                                                 <  1:33791/151‑1 (MQ=255)
                  aaaGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATc                                                                                                                                >  1:247139/1‑151 (MQ=255)
                       cacTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACt                                                                                                                           <  1:92329/151‑1 (MQ=255)
                             aGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGCTCTCGGCCTGAGCATCTGACTGTTCtt                                                                                                                     >  2:72779/1‑151 (MQ=255)
                                   aaaCACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTttatt                                                                                                                 >  2:330032/1‑149 (MQ=255)
                                         cGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGgaag                                                                                                         <  1:590130/151‑1 (MQ=255)
                                                gaTAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCgg                                                                                                  <  2:360999/151‑1 (MQ=255)
                                                   aaTTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGAcc                                                                                               <  2:156033/151‑1 (MQ=255)
                                                           ggggTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTActtct                                                                                       <  1:153746/151‑1 (MQ=255)
                                                                gtTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTAt                                                                                   <  2:79768/150‑1 (MQ=255)
                                                                 ttGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATcc                                                                                 >  2:231465/1‑151 (MQ=255)
                                                                           aGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAg                                                                       <  1:273267/151‑1 (MQ=255)
                                                                                gAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGCCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCtg                                                                   >  1:110036/1‑150 (MQ=255)
                                                                                     tACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGGCGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTgtgggtg                                                             >  1:489171/1‑151 (MQ=255)
                                                                                      aCAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTgtgg                                                                <  2:399875/147‑1 (MQ=255)
                                                                                        agcagcTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGgtggtg                                                          >  2:599286/1‑151 (MQ=255)
                                                                                         gcagcTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGa                                                         <  2:246712/151‑1 (MQ=255)
                                                                                         gcagcTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGa                                                         <  2:271440/151‑1 (MQ=255)
                                                                                                                  cgtcgGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTgg                                 <  2:554452/150‑1 (MQ=255)
                                                                                                                     cgGTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAAc                             <  1:518279/151‑1 (MQ=255)
                                                                                                                      ggTAACCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACa                            <  1:330039/151‑1 (MQ=255)
                                                                                                                         aaCCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGa                          <  1:599301/150‑1 (MQ=255)
                                                                                                                         aaCCCGTTTCAGCTCTTTCTGGAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGa                          <  2:4761/150‑1 (MQ=255)
                                                                                                                                             ggAGACGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCt      <  2:126354/150‑1 (MQ=255)
                                                                                                                                              gagaCGACGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGt    <  2:67373/151‑1 (MQ=255)
                                                                                                                                                 acgacgGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGCGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGttt  >  1:531484/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCTAAACAACAGACCGAGAAAGACACTGAAGTTCAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACAGCAGCTTTTTTTAATATGTCCCGTTCGTCGGTAACCCGTTTCAGCTCTTTCTGGAGACGGCGGATCTCGGCCTGAGCATCTGACTGTTCTTTATTAGTGGAAGAATCCGGACCGTACTTCTTTATCCAGGCATAAAGGCTGTGGGTGGTGATATCGAGACGTGTTGCAACGCTGGCAACAGAATAACCGAGATCAACAACCTGTTT  >  NC_000913/280023‑280317

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