Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,068,952	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	2,068,952	0	C	T	100.0%	215.2 / NA	66	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/34);  total (32/34)

TGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAA  >  NC_000913/2068820‑2069068                                                                                                                                     |                                                                                                                     tGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATAtt                                                                                                                >  2:1289994/1‑139 (MQ=255) tGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATAtt                                                                                                                >  2:1331603/1‑139 (MQ=255)    acaacaaAC‑GTGCACAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCa                                                                                                            <  2:1219196/138‑1 (MQ=255)     caacaaACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCa                                                                                                            <  2:332980/139‑1 (MQ=255)      aacaaacaTTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAg                                                                                                           >  1:1627742/1‑139 (MQ=255)                  aGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACAt                                                                                                       >  2:343741/1‑131 (MQ=255)                  aGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACAt                                                                                                       <  1:343741/131‑1 (MQ=255)                     aCTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                            <  2:1647951/139‑1 (MQ=255)                         caccaATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAc                                                                                        >  1:1328929/1‑139 (MQ=255)                          accaATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                       <  1:1289994/139‑1 (MQ=255)                                       tAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                          <  2:282653/139‑1 (MQ=255)                                            gATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATTCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGTCATCTGTTATCACTTAACCCATTGCACGCCCGCTGcc                                                                     <  2:1160689/139‑1 (MQ=37)                                             aTAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTta                                                                                                                    >  2:581970/1‑91 (MQ=255)                                             aTAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACAt                                                                                                       >  2:871831/1‑104 (MQ=255)                                             aTAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACAt                                                                                                       <  1:871831/104‑1 (MQ=255)                                             aTAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCg                                                                    >  1:1672112/1‑139 (MQ=255)                                                    gCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                       >  2:6487/1‑113 (MQ=255)                                                    gCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                       <  1:6487/113‑1 (MQ=255)                                                        gACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTccc                                                         <  2:1627742/139‑1 (MQ=255)                                                        gACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTccc                                                         >  2:852819/1‑139 (MQ=255)                                                            tCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTgg                                                     <  1:471310/139‑1 (MQ=255)                                                             ccccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgc                                                    <  2:1672112/139‑1 (MQ=255)                                                             ccccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgc                                                    >  1:478734/1‑139 (MQ=255)                                                               ccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                          <  2:1075762/115‑1 (MQ=255)                                                               ccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCACTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                          >  1:1075762/1‑115 (MQ=255)                                                                  gAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcgagcg                                               >  1:809542/1‑139 (MQ=255)                                                                       gagaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGa                                                                                                          >  1:803925/1‑75 (MQ=255)                                                                       gagaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGa                                                                                                          <  2:803925/75‑1 (MQ=255)                                                                       gagaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                 >  2:405104/1‑84 (MQ=255)                                                                       gagaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                 <  1:405104/84‑1 (MQ=255)                                                                       gagaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                          >  2:457053/1‑107 (MQ=255)                                                                       gagaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                          <  1:457053/107‑1 (MQ=255)                                                                            aaCAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgc                                     <  1:852819/139‑1 (MQ=255)                                                                             aCAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcg                                    <  1:1331603/139‑1 (MQ=255)                                                                              cAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcgc                                   >  1:609058/1‑139 (MQ=255)                                                                               aGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGTCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCa                                  >  2:466168/1‑139 (MQ=255)                                                                               aGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCa                                  >  2:6891/1‑139 (MQ=255)                                                                                      tCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCa                                       <  1:1485982/127‑1 (MQ=255)                                                                                      tCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCa                                       >  2:1485982/1‑127 (MQ=255)                                                                                           ggCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACaaacacg                                                                          >  2:380276/1‑82 (MQ=255)                                                                                           ggCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                          <  1:380276/87‑1 (MQ=255)                                                                                           ggCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAAcc                                         >  1:470280/1‑120 (MQ=255)                                                                                           ggCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAAcc                                         <  2:470280/120‑1 (MQ=255)                                                                                            gCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGTCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                     <  1:466168/139‑1 (MQ=255)                                                                                            gCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                     <  2:809542/139‑1 (MQ=255)                                                                                              aGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAt                   <  2:114292/139‑1 (MQ=255)                                                                                               ggCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCa                                                                                                            >  2:735990/1‑49 (MQ=255)                                                                                               ggCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCa                                                                                                            <  1:735990/49‑1 (MQ=255)                                                                                                gCCTTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGcc                                                                            <  2:109160/80‑1 (MQ=255)                                                                                                gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGa                                                                                                          <  2:1021282/50‑1 (MQ=255)                                                                                                gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGa                                                                                                          >  1:1021282/1‑50 (MQ=255)                                                                                                gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                       <  2:598341/69‑1 (MQ=255)                                                                                                gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                       >  1:598341/1‑69 (MQ=255)                                                                                                gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGcc                                                                            >  1:109160/1‑80 (MQ=255)                                                                                                gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgc                                     >  1:1323960/1‑119 (MQ=255)                                                                                                gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgc                                     <  2:1323960/119‑1 (MQ=255)                                                                                                gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATg                       >  2:496815/1‑133 (MQ=255)                                                                                                gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATg                       <  1:496815/133‑1 (MQ=255)                                                                                                gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTc                 >  2:110611/1‑139 (MQ=255)                                                                                                  cGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTc                                                                                                              <  1:1537955/44‑1 (MQ=255)                                                                                                  cGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTc                                                                                                              >  2:1537955/1‑44 (MQ=255)                                                                                                               gCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGaa  <  1:110611/139‑1 (MQ=255)                                                                                                                cATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCg    >  2:181880/1‑136 (MQ=37)                                                                                                                cATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCg    <  1:181880/136‑1 (MQ=37)                                                                                                                 aTTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTgg                                                     >  1:1018411/1‑86 (MQ=255)                                                                                                                 aTTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTgg                                                     <  2:1018411/86‑1 (MQ=255)                                                                                                                                     |                                                                                                                     TGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAA  >  NC_000913/2068820‑2069068

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