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

ACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAA  >  NC_000913/2068826‑2069068                                                                                                                               |                                                                                                                     acaaacaTTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAg                                                                                                           <  1:889161/138‑1 (MQ=255) acaaacaTTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTAGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAg                                                                                                           >  2:889161/1‑138 (MQ=255)  caaacaTTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTtat                                                                                                                   >  2:892477/1‑129 (MQ=255)  caaacaTTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTtat                                                                                                                   <  1:892477/129‑1 (MQ=255)               aCTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCAcc                                                                                            <  1:1042996/139‑2 (MQ=255)                cTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACtt                                                                                           <  1:1042239/139‑1 (MQ=255)                  ccaccaATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTaa                                                                                         <  1:1647607/139‑1 (MQ=255)                   caccaATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACa                                                                                                        <  1:413648/123‑1 (MQ=255)                   caccaATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACa                                                                                                        >  2:413648/1‑123 (MQ=255)                        aTGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                            <  2:1585217/130‑1 (MQ=255)                        aTGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                            >  1:1585217/1‑130 (MQ=255)                          gCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTAc                                                                                 >  1:1136900/1‑139 (MQ=255)                           ccTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                 >  2:747450/1‑122 (MQ=255)                           ccTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                 <  1:747450/122‑1 (MQ=255)                             tCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTtat                                                                                                                   <  2:545360/102‑1 (MQ=255)                             tCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTtat                                                                                                                   >  1:545360/1‑102 (MQ=255)                             tCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGa                                                                                                          <  2:1005871/111‑1 (MQ=255)                             tCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGa                                                                                                          >  1:1005871/1‑111 (MQ=255)                               cATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGcc                                                                            <  2:1087747/139‑1 (MQ=255)                                       aTAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                            <  1:504349/115‑1 (MQ=255)                                       aTAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                            >  2:504349/1‑115 (MQ=255)                                                    cATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGt                                                       >  1:1002029/1‑139 (MQ=255)                                                       ccccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                          >  1:110716/1‑117 (MQ=255)                                                       ccccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                          <  2:110716/117‑1 (MQ=255)                                                                 gagaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAAc                                          >  1:857271/1‑139 (MQ=255)                                                                  agaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                   <  1:1069703/129‑1 (MQ=255)                                                                  agaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                   >  2:1069703/1‑129 (MQ=255)                                                                      aaCAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgc                                     <  2:1002029/139‑1 (MQ=255)                                                                        cAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATa                                            <  2:900254/130‑1 (MQ=255)                                                                        cAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATa                                            >  1:900254/1‑130 (MQ=255)                                                                             aaCTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCa                                  <  1:45038/135‑1 (MQ=255)                                                                             aaCTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCa                                  >  2:45038/1‑135 (MQ=255)                                                                              aCTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGATATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATg                             >  1:1394199/1‑139 (MQ=255)                                                                              aCTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                 >  1:439955/1‑71 (MQ=255)                                                                              aCTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                 <  2:439955/71‑1 (MQ=255)                                                                              aCTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATg                             >  1:1608080/1‑139 (MQ=255)                                                                                  cTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGa                         >  1:1148172/1‑139 (MQ=255)                                                                                     ggCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                   <  1:368422/110‑1 (MQ=255)                                                                                     ggCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcg                                                   >  2:368422/1‑110 (MQ=255)                                                                                      gCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                     <  2:1394199/139‑1 (MQ=255)                                                                                      gCAGGACGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTATCCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                     <  2:1039636/139‑1 (MQ=255)                                                                                          gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGtt                                                                                                  <  1:469106/58‑1 (MQ=255)                                                                                          gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGtt                                                                                                  >  2:469106/1‑58 (MQ=255)                                                                                          gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                          <  1:1383096/82‑1 (MQ=255)                                                                                          gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                          >  2:1383096/1‑82 (MQ=255)                                                                                          gCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTc                 >  2:1236123/1‑139 (MQ=255)                                                                                            cGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                          <  2:54072/80‑1 (MQ=255)                                                                                            cGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                          >  1:54072/1‑80 (MQ=255)                                                                                            cGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGt               >  1:8621/1‑139 (MQ=255)                                                                                                aaCAAATCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTtata           <  2:1148172/139‑1 (MQ=255)                                                                                                    aaTCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                          <  1:830909/72‑1 (MQ=255)                                                                                                    aaTCCGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                          >  2:830909/1‑72 (MQ=255)                                                                                                        cGCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGa   >  1:1439021/1‑139 (MQ=255)                                                                                                         gCATTCGTGGCTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGaa  <  2:8621/139‑1 (MQ=255)                                                                                                                               |                                                                                                                     ACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAA  >  NC_000913/2068826‑2069068

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