Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,467,921	C→T	*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	1,467,921	0	C	T	100.0%	149.0 / NA	46	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/27);  total (19/27)

GGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGC  >  NC_000913/1467798‑1468040                                                                                                                            |                                                                                                                        gggTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACAt                                                                                                          <  1:1058650/139‑1 (MQ=255)     ttGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGt                                                                                                      <  2:948376/139‑1 (MQ=255)       gCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAAGTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGATCCCTGGATTTGCCCTTATATTTCCAGCCATCTGTTa                                                                                                    <  2:562359/139‑1 (MQ=255)            ggTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                               <  1:383211/139‑1 (MQ=255)            ggTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                               <  1:424852/139‑1 (MQ=255)                 ttCGGGTACAGGATAAAGCTTAGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                          <  2:1085248/139‑1 (MQ=255)                              tAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                             <  2:158947/139‑1 (MQ=255)                              tAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                             >  2:829746/1‑139 (MQ=255)                              tAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                             <  1:469271/139‑1 (MQ=255)                              tAAAGCTTTGCTGGAAATTGACAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                             <  2:76654/139‑1 (MQ=255)                                      tGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCa                                                                                                               >  2:323775/1‑97 (MQ=255)                                      tGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCa                                                                                                               <  1:323775/97‑1 (MQ=255)                                      tGCTGGAAATTGGCAGTGGCGCTACGCGAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGACCGCTGCCGCa                                                                     >  2:788780/1‑139 (MQ=255)                                             aaTTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTc                                                              >  1:606504/1‑139 (MQ=255)                                                     gTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                      <  2:227421/139‑1 (MQ=255)                                                        gcgcTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGc                                                                            <  2:258802/114‑1 (MQ=255)                                                        gcgcTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGc                                                                            >  1:258802/1‑114 (MQ=255)                                                                   aCAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgc                                        >  1:424095/1‑139 (MQ=255)                                                                   aCAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgc                                        <  2:594217/139‑1 (MQ=255)                                                                    cAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcg                                       <  1:829746/139‑1 (MQ=255)                                                                       tGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCAc                                    <  2:333181/139‑1 (MQ=255)                                                                        gCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                      >  2:626561/1‑120 (MQ=255)                                                                        gCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                      <  1:626561/120‑1 (MQ=255)                                                                        gCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACt                                   >  2:724270/1‑139 (MQ=255)                                                                         cAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTa                                  >  2:270476/1‑139 (MQ=255)                                                                         cAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTa                                  >  1:164645/1‑139 (MQ=255)                                                                            gATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGc                               >  2:544831/1‑139 (MQ=255)                                                                               atTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGa                                                                   <  2:736464/100‑1 (MQ=255)                                                                               atTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGa                                                                   >  1:736464/1‑100 (MQ=255)                                                                                ttGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGAt                           <  1:544831/139‑1 (MQ=255)                                                                                 tGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                      >  2:611229/1‑111 (MQ=255)                                                                                 tGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                      <  1:611229/111‑1 (MQ=255)                                                                                   aCAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGATCGATAACCCAGCGCACTATGCGGATGcc                        <  1:788780/139‑1 (MQ=255)                                                                                   aCAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                        <  2:606504/139‑1 (MQ=255)                                                                                        tttGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCg                   >  1:269836/1‑139 (MQ=255)                                                                                              catcatGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATaa             >  1:600839/1‑139 (MQ=255)                                                                                                 catGGGACAAGAACCCTGGATTTGCCCTTATATTTCCa                                                                                                               >  2:117813/1‑38 (MQ=255)                                                                                                 catGGGACAAGAACCCTGGATTTGCCCTTATATTTCCa                                                                                                               <  1:117813/38‑1 (MQ=255)                                                                                                    gggACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCg       <  2:424095/139‑1 (MQ=255)                                                                                                         aaGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                      <  1:878280/87‑1 (MQ=255)                                                                                                         aaGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                      >  2:878280/1‑87 (MQ=255)                                                                                                         aaGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGc  <  2:269836/139‑1 (MQ=25)                                                                                                           gAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGcc                                                                        <  1:80183/67‑1 (MQ=255)                                                                                                           gAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGcc                                                                        >  2:80183/1‑67 (MQ=255)                                                                                                            aaCCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                          >  1:83106/1‑48 (MQ=255)                                                                                                            aaCCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                          <  2:83106/48‑1 (MQ=255)                                                                                                                            |                                                                                                                        GGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGC  >  NC_000913/1467798‑1468040

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