Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,519,014 G→A pseudogene (294/504 nt) insB1 → IS1 transposase B,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,519,0140GA98.4% 185.4 / ‑6.8 64pseudogene (294/504 nt)insB1IS1 transposase B,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (40/23);  minor base T (1/0);  total (41/23)
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

CGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGA  >  NC_000913/4518889‑4519144
                                                                                                                               |                                                                                                                                  
cgcTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                                       >  1:182490/1‑139 (MQ=255)
     gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                  >  2:544534/1‑139 (MQ=255)
     gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                  <  1:502047/139‑1 (MQ=255)
     gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                  <  1:4656/139‑1 (MQ=255)
                    tGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGc                                                                                                   <  1:496272/139‑1 (MQ=255)
                          gATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAg                                                                                             >  2:302677/1‑139 (MQ=255)
                               ggaggaCGGTTGTGGCGCACGTCTTCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGaaa                                                                                        >  1:139531/1‑139 (MQ=255)
                                    aCGGTTGTGGCGCACGTCTTCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGc                                                                                   >  1:154135/1‑139 (MQ=255)
                                       gTTGTGGCGCACGTCTTCGGTGAACGCACACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCa                                                                                  <  1:661695/139‑1 (MQ=255)
                                        ttgttgCGCACGTATTCGGTGAACG‑‑CACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGta                                                                              >  1:741769/1‑138 (MQ=14)
                                        ttGTGGCGCACGTCTTCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGt                                                                               >  1:37946/1‑139 (MQ=255)
                                           ttgcgcACGTATTCGGTGAACG‑‑CACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATc                                                                           >  1:418152/3‑139 (MQ=14)
                                           tgGCGCACGTCTTCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                  >  2:194572/1‑101 (MQ=255)
                                           tgGCGCACGTCTTCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                  <  1:194572/101‑1 (MQ=255)
                                                cACGTCTTCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCa                                                                       >  1:25919/1‑139 (MQ=255)
                                                cACGTCTTCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCa                                                                       >  1:437600/1‑139 (MQ=255)
                                                  cGTCTTCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAg                                                                     >  1:28407/1‑139 (MQ=255)
                                                     cTTCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGt                                                                  >  1:145577/1‑139 (MQ=255)
                                                     attCGGTGAACG‑‑CACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctgg                                                                                                                     <  2:89071/86‑1 (MQ=17)
                                                     attCGGTGAACG‑‑CACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctgg                                                                                                                     >  1:89071/2‑87 (MQ=17)
                                                      ttCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGtt                                                                 <  2:145577/139‑1 (MQ=255)
                                                             aaCG‑‑CACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                                   <  2:53197/81‑1 (MQ=11)
                                                             aaCG‑‑CACTATGGCGACGCTGGGGCGTCTTATGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                                   >  1:53197/1‑81 (MQ=11)
                                                               cG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAgc                                                        >  1:797651/1‑139 (MQ=255)
                                                               cG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAgc                                                        >  1:172788/1‑139 (MQ=255)
                                                               cG‑‑CACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGATGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAgcg                                                       >  2:660193/1‑139 (MQ=11)
                                                               cG‑‑CACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                                        <  1:148442/122‑1 (MQ=17)
                                                               cG‑‑CACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                                        >  2:148442/1‑122 (MQ=14)
                                                                    aCTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCa                                                     <  1:582690/139‑1 (MQ=255)
                                                                     ctctGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCAt                                                    >  2:44202/1‑139 (MQ=255)
                                                                       atggcgacgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCaa                                                                      >  2:454413/10‑118 (MQ=17)
                                                                       atGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCaa                                                                      <  1:454413/117‑1 (MQ=17)
                                                                              cgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                    <  2:93732/63‑1 (MQ=255)
                                                                              cgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                    >  1:93732/3‑65 (MQ=255)
                                                                              cgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCggtat                                          >  2:760333/3‑135 (MQ=11)
                                                                              cgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCggtat                                          >  1:445121/3‑135 (MQ=11)
                                                                              cacTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                 >  1:322225/1‑133 (MQ=255)
                                                                              cacTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                 <  2:322225/133‑1 (MQ=255)
                                                                              cacTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACa                                           >  1:538850/1‑139 (MQ=255)
                                                                                    aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatct                                     >  1:385186/1‑139 (MQ=255)
                                                                                    aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatct                                     >  2:67257/1‑139 (MQ=255)
                                                                                    aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatct                                     >  1:292612/1‑139 (MQ=255)
                                                                                           tctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                 >  1:37676/1‑119 (MQ=255)
                                                                                           tctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                 <  2:37676/119‑1 (MQ=255)
                                                                                               aGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGac                          <  2:538850/139‑1 (MQ=255)
                                                                                                  ctgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGacaac                       >  1:97814/1‑139 (MQ=255)
                                                                                                   tgctgTCGGCCTTTGAGGTCGTGGTATATATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGacaaca                      >  2:496335/1‑139 (MQ=255)
                                                                                                   tgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGacaaca                      >  2:642680/1‑139 (MQ=255)
                                                                                                   tgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGacaaca                      >  1:400667/1‑139 (MQ=255)
                                                                                                    gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                  <  2:188749/46‑1 (MQ=255)
                                                                                                    gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                  >  1:188749/1‑46 (MQ=255)
                                                                                                    gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACAt                     >  1:31075/1‑139 (MQ=255)
                                                                                                    gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACAt                     >  1:134478/1‑139 (MQ=255)
                                                                                                    gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACAt                     >  2:753687/1‑139 (MQ=255)
                                                                                                      tgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCa                                                                                  <  2:758280/76‑1 (MQ=255)
                                                                                                      tgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCa                                                                                  >  1:758280/1‑76 (MQ=255)
                                                                                                         cGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGc                <  1:67257/139‑1 (MQ=255)
                                                                                                          ggCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAg                                                >  2:333052/1‑106 (MQ=11)
                                                                                                          ggCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAg                                                <  1:333052/106‑1 (MQ=1)
                                                                                                                tGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTg         <  2:25919/139‑1 (MQ=255)
                                                                                                                tGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTg         <  2:437600/139‑1 (MQ=255)
                                                                                                                    gTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGAc     >  1:496598/1‑139 (MQ=255)
                                                                                                                     tCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACg    <  1:44202/139‑1 (MQ=255)
                                                                                                                       gTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGa  <  1:496335/139‑1 (MQ=255)
                                                                                                                               |                                                                                                                                  
CGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACG‑‑CACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGA  >  NC_000913/4518889‑4519144

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