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

CTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGC  >  NC_000913/1467793‑1468040                                                                                                                                 |                                                                                                                        ctctGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGTGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCa                                                                                                               <  2:970199/139‑1 (MQ=255)      gggTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACAt                                                                                                          <  2:1369603/139‑1 (MQ=255)      gggTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACAt                                                                                                          <  1:245317/139‑1 (MQ=255)            gCCGGGGTTATTCGGGTACAGGATAAAGCGTTGCTGGAAAATGGCAATGGTGCTACGCTAACAAGGCAGGATATTGATGGTTTTGAACAACATGGGACAAGAGCTCTGGATTTGCCCTTATAATTCTAGCCATCTGTTa                                                                                                    <  2:1451892/139‑1 (MQ=255)                          ggTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCAtt                                                                                      >  2:1148312/1‑139 (MQ=255)                                ggATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGc                                                                                >  2:1280033/1‑139 (MQ=255)                                 gATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGcc                                                                               >  1:66547/1‑139 (MQ=255)                                     aaGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCt                                                                           >  1:518438/1‑139 (MQ=255)                                      aGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGAc                                                                                                            <  2:435455/105‑1 (MQ=255)                                      aGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGAc                                                                                                            >  1:435455/1‑105 (MQ=255)                                           tGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGa                                                                                                             <  2:803037/99‑1 (MQ=255)                                           tGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGa                                                                                                             >  1:803037/1‑99 (MQ=255)                                           tGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTAt                                                                                                   <  1:939987/109‑1 (MQ=255)                                           tGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTAt                                                                                                   >  2:939987/1‑109 (MQ=255)                                           tGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                               <  2:697736/113‑1 (MQ=255)                                           tGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                               >  1:697736/1‑113 (MQ=255)                                           tGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCa                                                                     >  1:1331258/1‑139 (MQ=255)                                                          gTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGTACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                      <  1:1148312/139‑1 (MQ=255)                                                          gTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTc                                                                                                                 >  2:1320924/1‑80 (MQ=255)                                                          gTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTc                                                                                                                 <  1:1320924/80‑1 (MQ=255)                                                          gTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                      <  2:66547/139‑1 (MQ=255)                                                                   cgctaACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAAc                                             >  2:1533174/1‑139 (MQ=255)                                                                        aCAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTta                                                                                                                       <  1:1336411/60‑1 (MQ=255)                                                                         cAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcg                                       <  1:1533174/139‑1 (MQ=255)                                                                                   tatTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                          >  1:1450450/1‑78 (MQ=255)                                                                                   tatTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                          <  2:1450450/78‑1 (MQ=255)                                                                                    atTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                    >  2:723193/1‑67 (MQ=255)                                                                                    atTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                    <  1:723193/67‑1 (MQ=255)                                                                                      tGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATg                          >  1:577253/1‑139 (MQ=255)                                                                                        aCAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                        <  2:1331258/139‑1 (MQ=255)                                                                                        aCAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                        <  2:324483/139‑1 (MQ=255)                                                                                         cAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCa                       <  2:577253/139‑1 (MQ=255)                                                                                             tttGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCg                   >  2:1138374/1‑139 (MQ=255)                                                                                               tGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGtt                 >  1:1419796/1‑139 (MQ=255)                                                                                                        tGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCa                                     >  2:1124735/1‑110 (MQ=255)                                                                                                        tGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCa                                     <  1:1124735/110‑1 (MQ=255)                                                                                                           gACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGaa     <  1:1138374/139‑1 (MQ=38)                                                                                                              aaGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGc  <  1:1280033/139‑1 (MQ=25)                                                                                                                gAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTa                                  >  2:529207/1‑105 (MQ=25)                                                                                                                gAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTa                                  <  1:529207/105‑1 (MQ=25)                                                                                                                 aaCCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                        <  2:1531824/114‑1 (MQ=18)                                                                                                                 aaCCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                        >  1:1531824/1‑114 (MQ=18)                                                                                                                    ccTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCt                                                                           <  2:40615/60‑1 (MQ=9)                                                                                                                    ccTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCt                                                                           >  1:40615/1‑60 (MQ=9)                                                                                                                                 |                                                                                                                        CTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGC  >  NC_000913/1467793‑1468040

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