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% 129.6 / NA 43*302* (TAG→TAA) insD1IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/22);  total (21/22)

CGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACC  >  NC_000913/382461‑382705
                                                                                                                      |                                                                                                                              
cGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTAt                                                                                                            >  1:247522/1‑139 (MQ=25)
        cATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGt                                                                                                    <  2:260104/139‑1 (MQ=255)
             taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCACGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTgcg                                                                                               >  1:35665/1‑139 (MQ=255)
                    tGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAAt                                                                                          >  1:1146256/1‑137 (MQ=255)
                    tGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAAt                                                                                          <  2:1146256/137‑1 (MQ=255)
                    tGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTCTGCGCCAATCa                                                                                        <  1:682987/139‑1 (MQ=255)
                     ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTgcg                                                                                               <  1:1423820/131‑1 (MQ=255)
                     ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTgcg                                                                                               >  2:1423820/1‑131 (MQ=255)
                     ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCAc                                                                                       >  1:532507/1‑139 (MQ=255)
                     ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCAc                                                                                       >  2:1181073/1‑139 (MQ=255)
                         tCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAAGGCAAATCCAATTACCTATCAGACAGTTTGCGCCAATCACTAcc                                                                                   <  2:35665/139‑1 (MQ=255)
                            gCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAAtt                                                                                                                  <  1:1645789/105‑1 (MQ=11)
                            gCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAAtt                                                                                                                  >  2:1645789/1‑105 (MQ=11)
                              aTAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAcc                                                                              <  2:247522/139‑1 (MQ=255)
                              aTAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAcc                                                                              <  1:1181073/139‑1 (MQ=255)
                                 gTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTg                                                                           <  1:1629363/139‑1 (MQ=255)
                                         ggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACgg                                                                   >  2:223839/1‑139 (MQ=255)
                                            tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAAtt                                                                                                                  <  2:89558/89‑1 (MQ=14)
                                            tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAAtt                                                                                                                  >  1:89558/1‑89 (MQ=14)
                                                   cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTAcc                                                                                   <  2:317601/113‑1 (MQ=255)
                                                   cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTAcc                                                                                   >  1:317601/1‑113 (MQ=255)
                                                   cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGc                                                                  >  2:364585/1‑130 (MQ=255)
                                                   cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGc                                                                  <  1:364585/130‑1 (MQ=255)
                                                   cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACa                                                         >  1:255667/1‑139 (MQ=255)
                                                                tCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGAGTATTc                                            >  1:349621/1‑139 (MQ=255)
                                                                   gCGGCAGCGGACTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTcgcg                                         >  1:616685/1‑139 (MQ=255)
                                                                     ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATaa                                                                      >  1:884401/1‑108 (MQ=255)
                                                                     ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATaa                                                                      <  2:884401/108‑1 (MQ=255)
                                                                         gCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATa                                   <  2:255667/139‑1 (MQ=255)
                                                                         gCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATa                                   <  2:349621/139‑1 (MQ=255)
                                                                                       ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCg                                                                                  <  2:1244453/78‑1 (MQ=255)
                                                                                       ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCg                                                                                  >  1:1244453/1‑78 (MQ=255)
                                                                                       ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGc                     >  2:445924/1‑139 (MQ=255)
                                                                                            aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATg                                                            <  1:409139/95‑1 (MQ=255)
                                                                                            aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATg                                                            >  2:409139/1‑95 (MQ=255)
                                                                                            aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGa                                                           <  2:1566790/96‑1 (MQ=255)
                                                                                            aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGa                                                           >  1:1566790/1‑96 (MQ=255)
                                                                                                    cAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCaata        <  2:616685/139‑1 (MQ=255)
                                                                                                          tCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTAcc  <  1:445924/139‑1 (MQ=255)
                                                                                                            tGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCa                                                                                        >  2:1590298/1‑51 (MQ=255)
                                                                                                            tGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCa                                                                                        <  1:1590298/51‑1 (MQ=255)
                                                                                                            tGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAcc                                                                              >  2:475941/1‑61 (MQ=255)
                                                                                                            tGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAcc                                                                              <  1:475941/61‑1 (MQ=255)
                                                                                                                      |                                                                                                                              
CGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACC  >  NC_000913/382461‑382705

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