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%	114.0 / NA	36	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/16);  total (20/16)

GAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGC  >  NC_000913/382457‑382713                                                                                                                           |                                                                                                                                       gagGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTAc                                                                                                                        >  1:81247/1‑139 (MQ=17)       ttCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCa                                                                                                                  >  1:387197/1‑139 (MQ=37)        tCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAg                                                                                                                 >  1:167769/1‑139 (MQ=38)                  taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTACCAGGCAGTTTgcg                                                                                                       <  2:664624/139‑1 (MQ=255)                       aaTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCTCCAAt                                                                                                  >  1:139440/1‑139 (MQ=255)                        aTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATc                                                                                                 >  2:8147/1‑139 (MQ=255)                          ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCAc                                                                                               >  1:743884/1‑139 (MQ=255)                             aTCCGCATAGTGCGCTGGGTTATCGCTCGACACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTAc                                                                                            >  2:89951/1‑139 (MQ=255)                                 gCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGa                                                                                        >  2:542769/1‑139 (MQ=255)                                         cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAAt                                                                                >  1:430640/1‑139 (MQ=255)                                              ggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACgg                                                                           >  1:609320/1‑139 (MQ=255)                                                     gctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATg                                                                    <  1:8147/139‑1 (MQ=255)                                                                          ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGcc                                               <  2:743884/139‑1 (MQ=255)                                                                          ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGcc                                               >  2:607085/1‑139 (MQ=255)                                                                             aGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGAt                                            <  1:380131/139‑1 (MQ=255)                                                                               cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaa                                          >  1:396523/1‑139 (MQ=255)                                                                                 ggCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATaaca                                        <  2:387197/139‑1 (MQ=255)                                                                                     tGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACtt                                    <  2:139440/139‑1 (MQ=255)                                                                                            ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCt                               >  1:35552/1‑137 (MQ=255)                                                                                            ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCt                               <  2:35552/137‑1 (MQ=255)                                                                                               tAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGt                          <  2:165849/139‑1 (MQ=255)                                                                                                 aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGttt                        >  1:52108/1‑139 (MQ=255)                                                                                                    gATAACAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTtcat                     <  2:52108/139‑1 (MQ=255)                                                                                                       aaCAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCa                                                                         <  1:631355/84‑1 (MQ=255)                                                                                                       aaCAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCa                                                                         >  2:631355/1‑84 (MQ=255)                                                                                                         cAGATGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCaata                <  2:396523/139‑1 (MQ=255)                                                                                                            aTGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGAAAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAAtt             >  1:347296/1‑139 (MQ=255)                                                                                                             tGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTtcat                     >  2:134964/1‑130 (MQ=255)                                                                                                             tGTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTtcat                     <  1:134964/130‑1 (MQ=255)                                                                                                              gTCTGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTAc           >  1:273031/1‑139 (MQ=255)                                                                                                                 tGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAAt                                                                                >  2:11301/1‑67 (MQ=255)                                                                                                                 tGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAAt                                                                                <  1:11301/67‑1 (MQ=255)                                                                                                                 tGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACt                                     <  2:724077/110‑1 (MQ=255)                                                                                                                 tGGAAATATAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACt                                     >  1:724077/1‑110 (MQ=255)                                                                                                                      atatAAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTg   <  2:81247/139‑1 (MQ=255)                                                                                                                       tataAGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGc  <  1:752801/139‑1 (MQ=255)                                                                                                                           |                                                                                                                                       GAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGC  >  NC_000913/382457‑382713

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