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_000913	382,579	0	G	A	100.0%	179.0 / NA	56	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/26);  total (30/26)

GTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCT  >  NC_000913/382462‑382711                                                                                                                      |                                                                                                                                     gTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATc                                                                                                                 >  2:399078/1‑139 (MQ=37)  ttCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCa                                                                                                                >  2:1099284/1‑139 (MQ=37)             taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTgcg                                                                                                     <  1:1099284/139‑1 (MQ=255)                  aaTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAAt                                                                                                <  2:134593/139‑1 (MQ=255)                   aTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATc                                                                                               >  1:1485287/1‑139 (MQ=255)                     ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCAc                                                                                             >  2:322398/1‑139 (MQ=255)                               tAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAccc                                                                                   <  1:322398/139‑1 (MQ=255)                               tAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGAccc                                                                                   <  1:835997/139‑1 (MQ=255)                                    cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTAcc                                                                                         >  2:35211/1‑128 (MQ=255)                                    cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTAcc                                                                                         <  1:35211/128‑1 (MQ=255)                                      cTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATaa                                                                            <  2:855073/139‑1 (MQ=255)                                      cTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGACAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATaa                                                                            >  2:581471/1‑139 (MQ=255)                                      cTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGACAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATaa                                                                            >  2:867997/1‑139 (MQ=255)                                       tGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATCACAGATGTCTGGAAATATAAGGACTAATCCAATTACCTATCAGGCAGTTTGCGCCGATCACTACCAGACCCTGAATAAc                                                                           >  2:1012680/1‑139 (MQ=18)                                        gggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAAt                                                                                                <  1:596114/117‑1 (MQ=255)                                        gggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAAt                                                                                                >  2:596114/1‑117 (MQ=255)                                                gctcgcCACGGGAATATCTGAGGCAGCGGGCTTGTAATGGGTTAAGTGATACCAGATGTCTGGAAATATAAGGGCAAATCCAATTACCGATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATg                                                                  <  1:934373/139‑1 (MQ=255)                                                   cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACgg                                                                         <  2:494988/129‑1 (MQ=255)                                                   cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACgg                                                                         >  1:494988/1‑129 (MQ=255)                                                   cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACa                                                               >  2:868836/1‑139 (MQ=255)                                                      cACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAAt                                                            <  2:638827/139‑1 (MQ=255)                                                                  tGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTcgc                                                >  1:749401/1‑139 (MQ=255)                                                                     ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAgg                                                                                                              <  1:256573/74‑1 (MQ=255)                                                                     ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAgg                                                                                                              >  2:256573/1‑74 (MQ=255)                                                                      gCAGCGGGCTTGTAATGGGTTAAGCGATAACAGATGTTTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGa                                                                                      <  1:308651/97‑1 (MQ=21)                                                                      gCAGCGGGCTTGTAATGGGTTAAGCGATAACAGATGTTTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGa                                                                                      >  2:308651/1‑97 (MQ=21)                                                                            ggCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaaca                                      >  2:860886/1‑139 (MQ=255)                                                                                      gggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTg                            >  2:56742/1‑139 (MQ=255)                                                                                       ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGt                                                                                                          <  1:7234/60‑1 (MQ=255)                                                                                       ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGt                                                                                                          >  2:7234/1‑60 (MQ=255)                                                                                       ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGc                                                                        <  1:164141/94‑1 (MQ=255)                                                                                       ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGc                                                                        >  2:164141/1‑94 (MQ=255)                                                                                       ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGc                           >  1:1487321/1‑139 (MQ=255)                                                                                       ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGc                           >  2:1329306/1‑139 (MQ=255)                                                                                       ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGc                           >  1:69033/1‑139 (MQ=255)                                                                                         ttAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAg                         >  1:612237/1‑139 (MQ=255)                                                                                            aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGttt                      >  2:773350/1‑139 (MQ=255)                                                                                            aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGttt                      >  2:82110/1‑139 (MQ=255)                                                                                               gATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTtcat                   <  2:612237/139‑1 (MQ=255)                                                                                               gATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTtcat                   <  1:773350/139‑1 (MQ=255)                                                                                               gATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTtcat                   <  1:868836/139‑1 (MQ=255)                                                                                                aTAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTtcatc                  <  2:1487321/139‑1 (MQ=255)                                                                                                  aaCAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAAc                                                                           >  2:1478000/1‑80 (MQ=255)                                                                                                  aaCAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAAc                                                                           <  1:1478000/80‑1 (MQ=255)                                                                                                        tGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTAcc                                                                                         <  2:185663/60‑1 (MQ=255)                                                                                                        tGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTAcc                                                                                         >  1:185663/1‑60 (MQ=255)                                                                                                         gTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACgg                                                                         >  1:1128512/1‑75 (MQ=255)                                                                                                         gTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACgg                                                                         <  2:1128512/75‑1 (MQ=255)                                                                                                           cTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCt       <  2:749401/139‑1 (MQ=255)                                                                                                              gAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCg                                                                                        >  2:1211662/1‑55 (MQ=255)                                                                                                              gAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCg                                                                                        <  1:1211662/55‑1 (MQ=255)                                                                                                                aaTATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCt  <  1:860886/139‑1 (MQ=255)                                                                                                                 atatAAGGGCAAATCCAATTACCTATCAGGCAGTTTg                                                                                                       <  2:802359/37‑1 (MQ=255)                                                                                                                 atatAAGGGCAAATCCAATTACCTATCAGGCAGTTTg                                                                                                       >  1:802359/1‑37 (MQ=255)                                                                                                                  tataAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATa                                                           <  2:240024/80‑1 (MQ=255)                                                                                                                  tataAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATa                                                           >  1:240024/1‑80 (MQ=255)                                                                                                                      |                                                                                                                                     GTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCT  >  NC_000913/382462‑382711

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: ‑