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%	81.1 / NA	27	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/14);  total (13/14)

CTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCT  >  NC_000913/1467793‑1468042                                                                                                                                 |                                                                                                                          ctctGGTGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCATGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCa                                                                                                                 <  2:1409143/139‑1 (MQ=255)              cGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATc                                                                                                    >  2:556251/1‑139 (MQ=255)                 ggTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                 <  1:556251/139‑1 (MQ=255)                  ggtATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAAGCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACtt                                                                                                >  2:903554/3‑139 (MQ=255)                                     aaGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                            <  1:677404/124‑1 (MQ=255)                                     aaGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCACTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                            >  2:677404/1‑124 (MQ=255)                                         tttGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCg                                                                         <  2:468266/139‑1 (MQ=255)                                           tGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGt                                                                                                        <  1:1233076/106‑1 (MQ=255)                                           tGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGt                                                                                                        >  2:1233076/1‑106 (MQ=255)                                           tGCTGGAAATTGGAAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCa                                                                       >  2:52143/1‑139 (MQ=255)                                                          gTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                        >  1:200736/1‑139 (MQ=255)                                                                ctacgctaACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGAt                                                  >  2:811275/1‑139 (MQ=255)                                                                        acacTGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTCCAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCa                                            >  2:970293/1‑137 (MQ=255)                                                                        aCAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCa                                            <  1:970293/137‑1 (MQ=255)                                                                         cAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcg                                         <  2:200736/139‑1 (MQ=255)                                                                                ggATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATg                                  >  2:1351563/1‑139 (MQ=255)                                                                                  atatTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTATAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCAc                                      <  2:188338/133‑1 (MQ=255)                                                                                  atatTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTATAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCAc                                      >  1:188338/1‑133 (MQ=255)                                                                                    atTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGa                              <  1:52143/139‑1 (MQ=255)                                                                                        aCAGTTTTGAACATCATGGGACAAGAACTCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                          <  1:903554/139‑1 (MQ=255)                                                                                           gTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAtt                       <  1:811275/139‑1 (MQ=255)                                                                                                      catGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGc                                                                        >  1:810066/1‑79 (MQ=255)                                                                                                      catGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGc                                                                        <  2:810066/79‑1 (MQ=255)                                                                                                           gACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGaa       <  1:1351563/139‑1 (MQ=255)                                                                                                                gAACCCTGGATTTGCCCTTATATTTCCAGTCATCTGTTATCACTTTACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTta                  >  2:869944/1‑123 (MQ=255)                                                                                                                gAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTta                  <  1:869944/123‑1 (MQ=21)                                                                                                                gAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCct  >  1:1350845/1‑139 (MQ=18)                                                                                                                                 |                                                                                                                          CTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCT  >  NC_000913/1467793‑1468042

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