Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,569,309 G→A pseudogene (312/921 nt) yjiP → pseudogene, transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,569,3090GA100.0% 166.9 / NA 52pseudogene (312/921 nt)yjiPpseudogene, transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/24);  total (28/24)

CCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCAC  >  NC_000913/4569206‑4569439
                                                                                                       |                                                                                                                                  
cccGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGc                                                                                                     <  2:1313884/135‑1 (MQ=255)
cccGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGc                                                                                                     >  1:1313884/1‑135 (MQ=255)
                         tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATTGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATg                                                                                 <  2:730538/130‑1 (MQ=255)
                         tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATTGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATg                                                                                 >  1:730538/1‑130 (MQ=255)
                         tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                        >  1:1346281/1‑139 (MQ=255)
                           tggtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCAtata                                                                                                                           <  1:985762/86‑1 (MQ=255)
                           tggtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCAtata                                                                                                                           >  2:985762/1‑86 (MQ=255)
                              tgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACgg                                                                   >  1:1077988/1‑139 (MQ=255)
                                           gagCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGcc                                                                                                                                <  1:1446556/65‑1 (MQ=255)
                                           gagCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGcc                                                                                                                                >  2:1446556/1‑65 (MQ=255)
                                              ccGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCt                                                                                                    >  2:36455/1‑90 (MQ=255)
                                              ccGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCt                                                                                                    <  1:36455/90‑1 (MQ=255)
                                                 cgAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGTTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGAt                                                                                  <  2:116691/105‑1 (MQ=255)
                                                 cgAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATa                                                                                                                 >  1:767048/1‑74 (MQ=255)
                                                 cgAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATa                                                                                                                 <  2:767048/74‑1 (MQ=255)
                                                 cgAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCGGCTACCGTTGGTCATCCCGAt                                                                                  >  1:116691/1‑105 (MQ=255)
                                                   aGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGcc                                                                                                                                >  2:700516/1‑57 (MQ=255)
                                                   aGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGcc                                                                                                                                <  1:700516/57‑1 (MQ=255)
                                                   aGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCa                                                                                                                      <  1:525199/67‑1 (MQ=255)
                                                   aGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCa                                                                                                                      >  2:525199/1‑67 (MQ=255)
                                                   aGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTg                                              >  1:383001/1‑139 (MQ=255)
                                                      aCATTCATATGGCCTTTCGCCTGCTGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCg                                                             <  2:902786/121‑1 (MQ=255)
                                                      aCATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATGACGGTAGCCg                                                             >  1:902786/1‑121 (MQ=255)
                                                                ggCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATcc                                                                                      >  2:1089813/1‑86 (MQ=255)
                                                                ggCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATcc                                                                                      <  1:1089813/86‑1 (MQ=255)
                                                                  ccTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAacg                                                                                                                                  >  2:181065/1‑40 (MQ=255)
                                                                  ccTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAacg                                                                                                                                  <  1:181065/40‑1 (MQ=255)
                                                                      tCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGcc                                                                                                                                >  1:667301/1‑38 (MQ=38)
                                                                      tCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGcc                                                                                                                                <  2:667301/38‑1 (MQ=255)
                                                                      tCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTAcc                                                                                                 >  2:475616/1‑69 (MQ=255)
                                                                      tCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTAcc                                                                                                 <  1:475616/69‑1 (MQ=255)
                                                                              tGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGc                   >  1:915923/1‑139 (MQ=255)
                                                                                    aTTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGccgaccc             >  2:435345/1‑139 (MQ=255)
                                                                                          tGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATg                                                                                 >  2:1234499/1‑65 (MQ=255)
                                                                                          tGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATg                                                                                 <  1:1234499/65‑1 (MQ=255)
                                                                                           ggCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCg                                                                                                <  1:399018/49‑1 (MQ=255)
                                                                                           ggCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCg                                                                                                >  2:399018/1‑49 (MQ=255)
                                                                                           ggCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGccgacccgac          <  2:274927/135‑1 (MQ=255)
                                                                                           ggCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGccgacccgac          >  1:274927/1‑135 (MQ=255)
                                                                                             cGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGc                                                                                                     >  2:315266/1‑42 (MQ=255)
                                                                                             cGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGc                                                                                                     <  1:315266/42‑1 (MQ=255)
                                                                                             cGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCg                                                             <  1:528581/82‑1 (MQ=255)
                                                                                             cGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCg                                                             >  2:528581/1‑82 (MQ=255)
                                                                                             cGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCt                                                      <  1:365228/89‑1 (MQ=255)
                                                                                             cGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCt                                                      >  2:365228/1‑89 (MQ=255)
                                                                                               gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTa                                                                              <  2:1203370/63‑1 (MQ=255)
                                                                                               gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTa                                                                              >  1:1203370/1‑63 (MQ=255)
                                                                                               gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                        <  2:158884/69‑1 (MQ=255)
                                                                                               gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                        >  1:158884/1‑69 (MQ=255)
                                                                                               gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCt                                               >  2:242330/1‑94 (MQ=255)
                                                                                               gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCt                                               <  1:242330/94‑1 (MQ=255)
                                                                                               gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTAc       <  2:393833/134‑1 (MQ=255)
                                                                                               gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTAc       >  1:393833/1‑134 (MQ=255)
                                                                                               gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCAc  <  2:383001/139‑1 (MQ=255)
                                                                                                       |                                                                                                                                  
CCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCAC  >  NC_000913/4569206‑4569439

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