Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 382,579 G→A *302* (TAG→TAA insD1 → IS2 transposase TnpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913382,5790GA100.0% 89.0 / NA 29*302* (TAG→TAAinsD1IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/13);  total (16/13)

ATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTC  >  NC_000913/382480‑382710
                                                                                                   |                                                                                                                                   
aTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTgcg                                                                                                    <  1:587600/133‑1 (MQ=255)
aTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTgcg                                                                                                    >  2:587600/1‑133 (MQ=255)
  ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCAc                                                                                            >  2:731296/1‑139 (MQ=255)
  ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCAc                                                                                            >  1:592278/1‑139 (MQ=255)
      tCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTAcc                                                                                        <  2:592278/139‑1 (MQ=255)
        cGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCgg                                                                                      <  1:731296/139‑1 (MQ=255)
                 cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAAt                                                                             >  2:777252/1‑139 (MQ=255)
                                                         ggCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGa                                                                >  1:242078/1‑112 (MQ=255)
                                                         ggCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGa                                                                <  2:242078/112‑1 (MQ=255)
                                                                    ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGc                          >  1:15177/1‑139 (MQ=255)
                                                                    ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGc                          >  2:147447/1‑139 (MQ=255)
                                                                         aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAAt                                                                                               <  1:513750/65‑1 (MQ=255)
                                                                         aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAAt                                                                                               >  2:513750/1‑65 (MQ=255)
                                                                         aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGttt                     >  1:437415/1‑139 (MQ=255)
                                                                            gATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTtcat                  <  2:437415/139‑1 (MQ=255)
                                                                              tAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTAcc                                                                                        >  1:220809/1‑67 (MQ=255)
                                                                              tAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTAcc                                                                                        <  2:220809/67‑1 (MQ=255)
                                                                              tAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACa                                                              >  2:296870/1‑93 (MQ=255)
                                                                              tAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACa                                                              <  1:296870/93‑1 (MQ=255)
                                                                                   gATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAgg                                                                                                             >  1:681043/1‑41 (MQ=255)
                                                                                   gATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAgg                                                                                                             <  2:681043/41‑1 (MQ=255)
                                                                                    aTGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCa                                                                                             <  1:717033/56‑1 (MQ=255)
                                                                                    aTGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCa                                                                                             >  2:717033/1‑56 (MQ=255)
                                                                                      gTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTAcc                                                                                        >  1:784578/1‑59 (MQ=255)
                                                                                      gTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTAcc                                                                                        <  2:784578/59‑1 (MQ=255)
                                                                                       tCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCAt                                                                  >  2:229327/1‑80 (MQ=255)
                                                                                       tCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCAt                                                                  <  1:229327/80‑1 (MQ=255)
                                                                                       tCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTAcc       >  1:786412/1‑139 (MQ=255)
                                                                                            aaaTATAAGGGCCAATCCAATTACCTATCCGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACTGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCGGCAGTTTCATCAATAATTACCTGGTc  <  2:719318/139‑1 (MQ=255)
                                                                                                   |                                                                                                                                   
ATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTC  >  NC_000913/382480‑382710

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: