Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,651,537 G→A pseudogene (693/693 nt) insD1 → IS2 transposase TnpB,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_0009131,651,5370GA100.0% 87.6 / NA 29pseudogene (693/693 nt)insD1IS2 transposase TnpB,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/13);  total (16/13)

TTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATA  >  NC_000913/1651421‑1651670
                                                                                                                    |                                                                                                                                     
ttCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTa                                                                                                                 >  2:886406/1‑139 (MQ=255)
 tCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTAc                                                                                                                >  1:90977/1‑139 (MQ=255)
    aGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGct                                                                                                             <  2:748969/139‑1 (MQ=255)
    aGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGct                                                                                                             <  1:886406/139‑1 (MQ=255)
                                         ttATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGa                                                                        <  2:379570/139‑1 (MQ=255)
                                                 cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGt                                                                                                                       >  2:670750/1‑84 (MQ=25)
                                                 cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGt                                                                                                                       <  1:670750/84‑1 (MQ=25)
                                                          aaTATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGTAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTc                                                       <  1:57215/139‑1 (MQ=255)
                                                              tCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGa                                                                        >  2:441746/1‑118 (MQ=255)
                                                              tCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGa                                                                        <  1:441746/118‑1 (MQ=255)
                                                                  cGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACt                                                                      >  1:677857/1‑116 (MQ=255)
                                                                  cGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACt                                                                      <  2:677857/116‑1 (MQ=255)
                                                                       gCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCaa                                          >  2:180854/1‑139 (MQ=255)
                                                                          ggCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCt                                       >  1:97092/1‑139 (MQ=255)
                                                                            cTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGct                                                                                                             >  1:797422/1‑67 (MQ=255)
                                                                            cTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGct                                                                                                             <  2:797422/67‑1 (MQ=255)
                                                                                      gTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTgg                           >  1:276558/1‑139 (MQ=255)
                                                                                         aaGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGTTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTc                                                       <  2:462156/108‑1 (MQ=255)
                                                                                         aaGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTc                                                       >  1:462156/1‑108 (MQ=255)
                                                                                          aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACAGATGAACTACTCTAATCTGTTTGTTACGCAATACAAACTGGCGTTATCGACTAAGCAATACAGTTTGCCTCCTCACGCAATCTGGAAAAATTTgggtag                       >  2:93545/1‑135 (MQ=37)
                                                                                               tAACAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCg                  >  2:132118/1‑139 (MQ=255)
                                                                                                aaCAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCTCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGa                 >  2:778638/1‑139 (MQ=255)
                                                                                                aaCAGATGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGa                 <  1:132118/139‑1 (MQ=255)
                                                                                                     aTGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGct                                                                                                             <  2:431152/42‑1 (MQ=255)
                                                                                                     aTGTCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGct                                                                                                             >  1:431152/1‑42 (MQ=255)
                                                                                                        tCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGctcct                                                                                                          >  2:733014/1‑42 (MQ=255)
                                                                                                        tCTGGAAATATAAGGGCAAATCCACAGGTGAACTACGctcct                                                                                                          <  1:733014/42‑1 (MQ=255)
                                                                                                          tGGAAATATAAGGGCAAATCCACGGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTaa       >  1:646262/1‑139 (MQ=255)
                                                                                                               atatAAGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATa  <  2:276558/139‑1 (MQ=255)
                                                                                                                    |                                                                                                                                     
TTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATA  >  NC_000913/1651421‑1651670

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