Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,518,721	A→G	G65G (GGA→GGG)	insA‑7 →	KpLE2 phage‑like element; IS1 repressor protein InsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,518,721	0	A	G	100.0%	179.0 / NA	54	G65G (GGA→GGG)	insA‑7	KpLE2 phage‑like element; IS1 repressor protein InsA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/32);  total (22/32)

AAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAA  >  NC_000913/4518591‑4518855                                                                                                                                   |                                                                                                                                       aaaGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGcc                                                                                                                               <  1:73879/140‑1 (MQ=255)       cTGCCGGACATCAGCTCTATCTCTGTTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGCCGCGCCAGTGCa                                                                                                                         <  1:252443/140‑1 (MQ=17)       cTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCa                                                                                                                         <  1:134872/140‑1 (MQ=255)       cTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCa                                                                                                                         <  2:409239/140‑1 (MQ=255)       cTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCa                                                                                                                         <  2:474171/140‑1 (MQ=255)       cTGCCGGACATCAGCACTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCa                                                                                                                         <  1:195614/140‑1 (MQ=255)       cTGCCGGAAATCAGCGCTATCTTTGCTCTCATTTCCGTAAAACATGGCCACTACAGTTCACTTCCACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCa                                                                                                                         <  2:350930/140‑1 (MQ=17)         gCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACg                                                                                                                       <  1:549288/140‑1 (MQ=255)                     cgcTATCTCTGCTCTTCTTGCCGTAGAACTTGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGt                                                                                                           <  1:296543/140‑1 (MQ=255)                               gCTCTCCTTGCCGTAAATCATGGCAACTACAGTTCTCTTCCACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAac                                                                                                 <  2:254441/140‑1 (MQ=255)                               gCTCTCCTTGCCGTAAAACATGGCACCTACAGTCCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAac                                                                                                 <  1:34373/140‑1 (MQ=255)                               gCTCTCCTGGCCGTAAAACATGCCAACTACTTTTCTCTTCCACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAac                                                                                                 <  2:353980/140‑1 (MQ=18)                                ctctCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAaca                                                                                                <  1:486252/140‑1 (MQ=255)                                  ctcCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCAGTATGGGAGTTGGCCTCAACACg                                                                                              >  1:255656/1‑140 (MQ=255)                                    ccTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGt                                                                                            >  1:75649/1‑140 (MQ=255)                                     cTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGTCAGTGCACGCAGTATGGGCGTTGGCCTCAACACGGtt                                                                                           >  1:254441/1‑140 (MQ=255)                                               aaCATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTTATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGACGTAGGACTCAACTCGGTTTTACGTCACt                                                                                 >  1:338512/1‑140 (MQ=18)                                                aCATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACtt                                                                                <  1:610236/140‑1 (MQ=255)                                                            cAGTTCACTTACACCGCCTCTCAGCGCGGAAAGCACCTGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAgg                                                                    <  1:233566/140‑1 (MQ=255)                                                              gTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGcc                                                                  >  2:535186/1‑140 (MQ=255)                                                               ttCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCg                                                                 <  2:75649/140‑1 (MQ=255)                                                               ttCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCg                                                                 <  2:686533/140‑1 (MQ=255)                                                               ttCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCg                                                                 >  1:686533/1‑140 (MQ=255)                                                                 cACTTAAACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCa                                                               <  1:452989/140‑1 (MQ=255)                                                                  aCTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTTGCCTCGACACGGTTTTACGTCACTTAATAAACTCAGGCCGCAg                                                              >  1:309228/1‑140 (MQ=255)                                                                    ttACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCATTGGCCTCCACAAGCTTTTACGTCAATTAAGAAACTCAGGCCGCAGTc                                                            >  2:344419/1‑140 (MQ=18)                                                                       cacCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGt                                                         <  1:531613/140‑1 (MQ=255)                                                                        acCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTa                                                        >  2:570395/1‑140 (MQ=255)                                                                         ccGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTaa                                                       >  1:704973/1‑140 (MQ=255)                                                                             ctctCAGCCCGGTAAGCTCCAGAAAATCATTGATATGGCTATGAATGGGGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTTGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTc                                                   <  2:309228/140‑1 (MQ=17)                                                                             ctctCAGCCCGGTAAGCACCAGTAACTCACTGATATGACCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTc                                                   <  1:270531/140‑1 (MQ=17)                                                                             ctctCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTc                                                   <  1:179562/140‑1 (MQ=255)                                                                             ctctCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTc                                                   >  2:55415/1‑140 (MQ=255)                                                                             ctctCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAAAGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCgg                                                          >  2:468031/1‑133 (MQ=255)                                                                             ctctCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAAAGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCgg                                                          <  1:468031/133‑1 (MQ=255)                                                                              tctcAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTATGTCACTTCAAAAACGCAGGCCCCAGTCGGTAACCTcg                                                  >  1:350407/1‑140 (MQ=18)                                                                              tctcAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGCCGCGTAAAAAACTCAGGCCGCAGTCGGTAACCTcg                                                  >  2:346751/1‑140 (MQ=255)                                                                                       ggTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACaa                                         <  2:170590/140‑1 (MQ=255)                                                                                         tAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTAGGTCACTTACAAATCTCAGGCCGCAGGCGGTAACCTCGCGCATACAAcc                                       >  2:222511/1‑140 (MQ=17)                                                                                         tAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTAGGTCACTTAAAAAACGCCGGCCGCAGTCGGTGCCCTCGCGGATACAAcc                                       >  2:355724/1‑140 (MQ=17)                                                                                         tAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATa                                            <  2:596082/135‑1 (MQ=255)                                                                                         tAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATa                                            >  1:596082/1‑135 (MQ=255)                                                                                            gCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCggg                                    <  1:570395/140‑1 (MQ=255)                                                                                            gCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCATCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCATTCGGTAACCTCGCGCATACAACCggg                                    <  2:180606/140‑1 (MQ=255)                                                                                              aCCAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCa                                  >  2:436570/1‑140 (MQ=255)                                                                                               ccAGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAg                                 <  1:541252/140‑1 (MQ=255)                                                                                                 aGAAAATCGTTTATATAGCCATGAATGGCGTGGGGGGTCGCGCCAGTGTACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAgtg                               <  2:338512/140‑1 (MQ=17)                                                                                                 aGAAAATCATTGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAgtg                               <  2:704973/140‑1 (MQ=255)                                                                                                           tGATATGGCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGt                     <  2:51414/140‑1 (MQ=255)                                                                                                                 ggCCATGAATGGCGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTgcgc               >  2:734625/1‑140 (MQ=255)                                                                                                                      tGAATGGCGTCGGGTGTAGCGCTAGTGCACGGATTTTGGTCGTTGGCGGCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGaaa          <  2:225288/140‑1 (MQ=18)                                                                                                                             cGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGa   >  1:697947/1‑140 (MQ=255)                                                                                                                             cGTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGa   >  2:618255/1‑140 (MQ=255)                                                                                                                              gTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGaa  >  1:429382/1‑140 (MQ=255)                                                                                                                              gTCGGGTGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGaa  <  1:55415/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                                       AAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCCTTGCCGTAAAACATGGCAACTACAGTTCACTTACACCGCCTCTCAGCCCGGTAAGCACCAGAAAATCATTGATATGGCCATGAATGGCGTCGGATGTCGCGCCAGTGCACGCATTATGGGCGTTGGCCTCAACACGGTTTTACGTCACTTAAAAAACTCAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAA  >  NC_000913/4518591‑4518855

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