Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,121,883	G→T	R135R (CGC→CGA)	hslV ←	peptidase component of the HslVU protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,121,883	0	G	T	100.0%	132.4 / NA	38	R135R (CGC→CGA)	hslV	peptidase component of the HslVU protease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (23/15);  total (23/15)

GGGAGATCCTTACGCTTTGTAGCTTAATTCTTCGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTC  >  NC_000913/4121748‑4122011                                                                                                                                        |                                                                                                                                 gggAGATCCTTACGCTTTGTAGCTTAATTCTTCGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgctcgcg                                                                                                                              >  1:681965/1‑140 (MQ=255) gggAGATCCTTACGCTTTGTAGCTTAATTCTTCGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgctcgcg                                                                                                                              >  1:430670/1‑140 (MQ=255)  ggAGATCCTTACGCTTTGTAGCTTAATTCTTCGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgctcgcgc                                                                                                                             >  2:240958/1‑140 (MQ=255)  ggAGATCCTTACGCTTTGTAGCTTAATTCTTCGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgctcgcgc                                                                                                                             >  1:672024/1‑140 (MQ=255)  ggAGATCCTTACGCTTTGTAGCTTAATTCTTCGATGGGGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGATCAGTGTTTTCTAACAgcgctcgcgc                                                                                                                             >  2:645243/1‑140 (MQ=255)    agaTCCTTACGCTTTGTAGCTTAATTCTTCGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGcgccg                                                                                                                           <  1:670333/140‑1 (MQ=255)          ttACGCTTTGTAGCTTAATTCTTCGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCt                                                                                                                     <  1:645243/140‑1 (MQ=255)          ttACGCTTTGTAGCTTAATTCTTCGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCt                                                                                                                     <  1:545647/140‑1 (MQ=255)          ttACGCTTTGTAGCTTAATTCTTCGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGCTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCt                                                                                                                     <  2:324928/140‑1 (MQ=255)                  tGTAGCTTAATTCTTCGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTaa                                                                                                             >  2:552173/1‑140 (MQ=255)                  tGTAGCTTAATTCTTCGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTaa                                                                                                             >  2:53476/1‑140 (MQ=255)                  tGTAGCTTAATTCTTCGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTaa                                                                                                             >  1:880063/1‑140 (MQ=255)                  tGTAGCTTAATTCTTCGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTaa                                                                                                             >  1:747703/1‑140 (MQ=255)                   gTAGCTTAATTCTTCGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAg                                                                                                            >  1:617589/1‑140 (MQ=255)                          aaTTCTTCGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGgccgcc                                                                                                     >  1:151627/1‑140 (MQ=255)                                tcGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGcc                                                                                               <  2:672024/140‑1 (MQ=255)                                  gATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGa                                                                                             <  1:288526/140‑1 (MQ=255)                                   aTGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGAt                                                                                            <  2:604595/140‑1 (MQ=255)                                       tgtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCa                                                                                        >  2:591531/1‑140 (MQ=255)                                          ggAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATa                                                                                     <  2:617589/140‑1 (MQ=255)                                             aaTGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGa                                                                                  >  1:526906/1‑140 (MQ=255)                                              aTGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGAt                                                                                 <  1:26263/140‑1 (MQ=255)                                                ggttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCg                                                                               >  2:42140/1‑140 (MQ=255)                                                 gttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGt                                                                              >  2:951211/1‑140 (MQ=255)                                                   tggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGttt                                                                            >  1:341253/1‑140 (MQ=255)                                                    ggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGtttt                                                                           <  1:648681/140‑1 (MQ=255)                                                           aTGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTg                                                                    >  2:419807/1‑140 (MQ=255)                                                                aaTGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcacca                                                               >  2:205308/1‑140 (MQ=255)                                                                 aTGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcaccac                                                              <  2:681965/140‑1 (MQ=255)                                                                    gcGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTc                                                           <  1:34165/139‑1 (MQ=255)                                                                      gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCAc                                                         >  2:128407/1‑140 (MQ=255)                                                                           cAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTa                                                    <  1:552173/140‑1 (MQ=255)                                                                           cAATATCCAACACCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTa                                                    >  1:733946/1‑140 (MQ=255)                                                                            aaTATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTAc                                                   <  2:151627/140‑1 (MQ=255)                                                                                     cGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACTGGtgatga                                          >  1:548015/1‑140 (MQ=255)                                                                                     cGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatga                                          >  1:60181/1‑140 (MQ=255)                                                                                     cGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatga                                          >  1:454138/1‑140 (MQ=255)                                                                                     cGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatga                                          >  1:195011/1‑140 (MQ=255)                                                                                      gCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTgatgat                                         >  2:651529/1‑140 (MQ=255)                                                                                          tttCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGc                                     >  1:198951/1‑140 (MQ=255)                                                                                                             gcTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGa                  >  2:207268/1‑140 (MQ=255)                                                                                                              cTAAGTTCAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGAc                 >  2:764222/1‑140 (MQ=255)                                                                                                                     cAGTGTTTTCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCcagc          <  1:240958/140‑1 (MQ=255)                                                                                                                             tCTAACAGCGCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTc  >  1:565322/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                 GGGAGATCCTTACGCTTTGTAGCTTAATTCTTCGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTC  >  NC_000913/4121748‑4122011

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