Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,187,415	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	3,187,415	0	G	A	100.0%	122.3 / NA	39	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/18);  total (21/18)

AGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGG  >  NC_000913/3187292‑3187548                                                                                                                            |                                                                                                                                      agagGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTa                                                                                                                        <  2:628635/139‑1 (MQ=17)      cGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTatcatc                                                                                                                   >  2:1109614/1‑139 (MQ=37)                    aaCGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTc                                                                                                     <  1:618473/139‑1 (MQ=255)                             cATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTAt                                                                                              <  1:1037368/137‑1 (MQ=255)                             cATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTAt                                                                                              >  2:1037368/1‑137 (MQ=255)                             cATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGc                                                                                            >  2:206106/1‑139 (MQ=255)                               tCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCg                                                                                          >  1:207857/1‑139 (MQ=255)                                 cGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGtt                                                                                        <  1:780850/139‑1 (MQ=255)                                         gcgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTa                                                                                                                        >  1:106403/1‑99 (MQ=18)                                         gcgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTa                                                                                                                        <  2:106403/99‑1 (MQ=18)                                          cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTcaac                                                                               >  1:985838/1‑139 (MQ=255)                                          cgcTGGGTTATCGCTCGCCACGGGAATATCTGCCGCAGCGGGCTTGTAATGGGTTAAGCGATAACAGATGTCTGGAAATATAACGGCAAATCCGGTCACCATCATTGGTTGTTCGTCACCTTAAGCCGTTGAACTcaac                                                                               >  2:391242/1‑139 (MQ=17)                                           gcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTcaaca                                                                              <  2:300151/139‑1 (MQ=255)                                            cTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTcaacaa                                                                             >  1:954615/1‑139 (MQ=255)                                                  tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTcaac                                                                               >  1:83262/1‑131 (MQ=255)                                                  tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTcaac                                                                               <  2:83262/131‑1 (MQ=255)                                                       ctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCg                                                                  >  2:440427/1‑139 (MQ=255)                                                       ctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCg                                                                  <  2:207857/139‑1 (MQ=255)                                                                ggAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAg                                                         <  2:323590/139‑1 (MQ=255)                                                                         gCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGaaa                                                >  1:618783/1‑139 (MQ=255)                                                                                cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGaa                                                 <  2:1095793/131‑1 (MQ=255)                                                                                cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGaa                                                 >  1:1095793/1‑131 (MQ=255)                                                                                cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTaa                                         <  2:985838/139‑1 (MQ=255)                                                                                cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTaa                                         >  1:944565/1‑139 (MQ=255)                                                                                       gTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTcc                                  >  1:857970/1‑139 (MQ=255)                                                                                             ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTcc                                  >  2:909741/1‑133 (MQ=255)                                                                                             ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTcc                                  <  1:909741/133‑1 (MQ=255)                                                                                                  aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGtt                                                                                                                         >  2:903609/1‑41 (MQ=255)                                                                                                  aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGtt                                                                                                                         <  1:903609/41‑1 (MQ=255)                                                                                                  aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAAc                       >  1:478817/1‑139 (MQ=255)                                                                                                  aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAAc                       >  2:295710/1‑139 (MQ=255)                                                                                                  aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAAc                       >  1:375058/1‑139 (MQ=255)                                                                                                      aTAACAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGAt                   <  2:857970/139‑1 (MQ=255)                                                                                                        aaCAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTg                     >  2:24543/1‑135 (MQ=255)                                                                                                        aaCAGATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTg                     <  1:24543/135‑1 (MQ=255)                                                                                                            gATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCg                                                                  >  2:333688/1‑86 (MQ=255)                                                                                                            gATGTCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCg                                                                  <  1:333688/86‑1 (MQ=255)                                                                                                                tCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCg                                                                                                       >  2:375875/1‑45 (MQ=255)                                                                                                                tCTGGAAATATAAGGGCAAATCCAGTTATCATCATTGGTTGTGCg                                                                                                       <  1:375875/45‑1 (MQ=255)                                                                                                                       atatAAGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCgg  <  2:305226/139‑1 (MQ=255)                                                                                                                            |                                                                                                                                      AGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAGTTATCATCATTGGTTGTGCGTCAGCTTATGCCGTTGAATTCAACAAAGATTTAATCGAAGCCGAAGATCGTGAAAACGTTAACCTTTCCCAATTTGAAACTGATGGCCAATTACCCGTCGG  >  NC_000913/3187292‑3187548

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