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

GAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGC  >  NC_000913/1467786‑1468032                                                                                                                                        |                                                                                                                gaATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTtat                                                                                                              <  2:329060/139‑1 (MQ=255)            ggggTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACa                                                                                                   <  1:159981/139‑1 (MQ=255)            ggggTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACa                                                                                                   >  2:545738/1‑139 (MQ=255)             gggTTTGCCGGGGTTATTCGGGTACAGGATAAAGCGTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACAt                                                                                                  <  1:561177/139‑1 (MQ=255)                                  gTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACAt                                                                                                  <  1:169200/118‑1 (MQ=255)                                  gTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACAt                                                                                                  >  2:169200/1‑118 (MQ=255)                                  gTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTa                                                                             >  1:147164/1‑139 (MQ=255)                                      aGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAg                                                                         >  2:232617/1‑139 (MQ=255)                                          tAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                     <  1:232617/139‑1 (MQ=255)                                          tAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                     <  1:68205/139‑1 (MQ=255)                                          tAAAGCGTTGCTGGAAATTCGAAGTGTCGCAACGCGAACAATGCAGGATATTTACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCGATTACAAGCCCg                                                                     <  1:753664/139‑1 (MQ=25)                                            aaGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACa                                                                                                   <  1:439530/107‑1 (MQ=255)                                            aaGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACa                                                                                                   >  2:439530/1‑107 (MQ=255)                                             aGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                     >  2:383859/1‑136 (MQ=255)                                             aGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                     <  1:383859/136‑1 (MQ=255)                                             aGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTg                                                                  <  1:668333/139‑1 (MQ=255)                                                  tGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCAc                                                                                        <  1:745081/112‑1 (MQ=255)                                                  tGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCAc                                                                                        >  2:745081/1‑112 (MQ=255)                                                  tGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGc                                                                 <  1:517185/135‑1 (MQ=255)                                                  tGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCa                                                             >  1:838627/1‑139 (MQ=255)                                                  tGCTGGAAATTGGCAGTGACGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGc                                                                 >  2:517185/1‑135 (MQ=255)                                                      ggAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGata                                                         <  2:147164/139‑1 (MQ=255)                                                                 gTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                              <  2:685228/139‑1 (MQ=255)                                                                      gctacgctaACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcgagcga                                         <  1:688776/139‑1 (MQ=255)                                                                            ctaACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGtt                                                                                             <  1:42097/81‑1 (MQ=255)                                                                            ctaACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGtt                                                                                             >  2:42097/1‑81 (MQ=255)                                                                               aCAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgc                                >  1:385280/1‑139 (MQ=255)                                                                               aCAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgc                                >  2:403317/1‑139 (MQ=255)                                                                                cAATGCAGGATATTGACAGTTTTGAACATCATGGTACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcg                               <  2:758415/139‑1 (MQ=255)                                                                                cAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcg                               <  2:385280/139‑1 (MQ=255)                                                                                  aTGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                     >  2:577262/1‑99 (MQ=255)                                                                                  aTGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                     <  1:577262/99‑1 (MQ=255)                                                                                    gcgggATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgc                                <  1:879606/131‑1 (MQ=255)                                                                                    gcgggATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgc                                >  2:879606/4‑134 (MQ=255)                                                                                    gCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                            >  1:139881/1‑74 (MQ=255)                                                                                    gCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                            <  2:139881/74‑1 (MQ=255)                                                                                         atatTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCg                      >  1:635733/1‑139 (MQ=255)                                                                                             tGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                              >  2:700664/1‑111 (MQ=255)                                                                                             tGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                              <  1:700664/111‑1 (MQ=255)                                                                                             tGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATg                  <  1:55172/139‑1 (MQ=255)                                                                                               gcAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                            >  2:581118/2‑63 (MQ=255)                                                                                               aCAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATTTTTCCAGACATCTGTTa                                                                                            <  1:581118/63‑1 (MQ=255)                                                                                               aCAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                <  2:635733/139‑1 (MQ=255)                                                                                               aCAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                <  1:495440/139‑1 (MQ=255)                                                                                                  gTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcgagcga                                         <  2:284497/111‑1 (MQ=255)                                                                                                  gTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcgagcga                                         >  1:284497/1‑111 (MQ=255)                                                                                                  gTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                >  2:200530/1‑136 (MQ=255)                                                                                                  gTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                <  1:200530/136‑1 (MQ=255)                                                                                                  gTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAtt             <  1:545738/139‑1 (MQ=255)                                                                                                      tGAACATCATGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGtt         >  1:78493/1‑139 (MQ=255)                                                                                                             catGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                  <  2:522097/59‑1 (MQ=255)                                                                                                             catGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                  >  1:522097/1‑59 (MQ=255)                                                                                                             catGGGACAAGAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGc  >  2:473697/1‑139 (MQ=255)                                                                                                                       gAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                       <  1:931382/44‑1 (MQ=255)                                                                                                                       gAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                       >  2:931382/1‑44 (MQ=255)                                                                                                                       gAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                              >  1:734770/1‑85 (MQ=255)                                                                                                                       gAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                              <  2:734770/85‑1 (MQ=37)                                                                                                                       gAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                <  2:504878/115‑1 (MQ=21)                                                                                                                       gAACCCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                >  1:504878/1‑115 (MQ=21)                                                                                                                          cccTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                  <  2:443967/46‑1 (MQ=255)                                                                                                                          cccTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                  >  1:443967/1‑46 (MQ=255)                                                                                                                                        |                                                                                                                GAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGC  >  NC_000913/1467786‑1468032

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