Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260479	724,768	+C	coding (82/669 nt)	lysE →	lysine efflux permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260479	724,766	1	.	C	100.0%	143.0 / NA	37	R27P (CGC→CCC)	lysE	lysine efflux permease
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (21/16);  total (21/16)

GCGTGATCCTGGGTGTAGTGGCGTGAGGTAGTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCG‑CCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTCTGCCTGGTG  >  AM260479/724635‑724878                                                                                                                                     |                                                                                                                 gCGTGATCCTGGGTGTAGTGGCGTGAGGTAGTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACgccg‑cc                                                                                                                <  1:484874/134‑1 (MQ=255)        ccTGGGTGTAGTGGCGTGAGGTAGTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCg                                                                                                          >  1:875965/1‑134 (MQ=255)                       gTGAGGTAGTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTc                                                                                           >  2:998963/1‑134 (MQ=255)                            gTAGTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTcc                                                                                      >  1:502516/1‑134 (MQ=255)                            gTAGTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATACACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTcc                                                                                      >  1:1030189/1‑134 (MQ=255)                              aGTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCt                                                                                    >  1:615064/1‑134 (MQ=255)                               gTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATACACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCt                                                                                    <  2:1030189/133‑1 (MQ=255)                                              tGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTtct                                                                                                    >  2:639033/1‑102 (MQ=255)                                              tGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTtct                                                                                                    <  1:639033/102‑1 (MQ=255)                                                         gcgggcATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGc                                                         >  2:910786/1‑134 (MQ=255)                                                          cgggcATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACgccg‑ccct                                                                                                              <  2:1031823/79‑3 (MQ=255)                                                          cgggcATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACgccg‑ccct                                                                                                              >  1:1031823/1‑77 (MQ=255)                                                             gcATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCa                                                                         >  1:10197/1‑114 (MQ=255)                                                             gcATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCa                                                                         <  2:10197/114‑1 (MQ=255)                                                              cATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCg                                                    <  1:998963/134‑1 (MQ=255)                                                              cATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCg                                                    <  2:262462/134‑1 (MQ=255)                                                                   cggcTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGtgctg                                               <  1:910786/134‑1 (MQ=255)                                                                      cTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGa                                            >  2:625942/1‑134 (MQ=255)                                                                           cccGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAgg                                       >  2:643849/1‑134 (MQ=255)                                                                                cGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTgc                                  <  2:875965/134‑1 (MQ=255)                                                                                                cccAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTg                                                                            <  2:901100/76‑1 (MQ=255)                                                                                                cccAAAAAAACAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTg                                                                            >  1:901100/1‑76 (MQ=255)                                                                                                   aaaaaaaCAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTg                                                                            <  2:212806/73‑1 (MQ=255)                                                                                                   aaaaaaaCAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTg                                                                            >  1:212806/1‑73 (MQ=255)                                                                                                   aaaaaaaCAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGa                           <  2:241994/122‑1 (MQ=255)                                                                                                   aaaaaaaCAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGa                           >  1:241994/1‑122 (MQ=255)                                                                                                   aaaaaaaCAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTc                 >  2:508427/1‑132 (MQ=255)                                                                                                   aaaaaaaCAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTc                 <  1:508427/132‑1 (MQ=255)                                                                                                   aaaaaaaCAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTg               >  2:116115/1‑134 (MQ=255)                                                                                                   aaaaaaaCAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTg               >  1:769405/1‑134 (MQ=255)                                                                                                     aaaaaCAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTc                 <  2:203663/130‑1 (MQ=255)                                                                                                     aaaaaCAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTc                 >  1:203663/1‑130 (MQ=255)                                                                                                     aaaaaCAGCCGTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTggg             >  2:814499/1‑134 (MQ=255)                                                                                                               gTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCg                                                    >  1:876661/1‑85 (MQ=255)                                                                                                               gTTCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCg                                                    <  2:876661/85‑1 (MQ=255)                                                                                                                ttCCACATCCAATGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTCTGCCTGgtg  <  1:625942/134‑1 (MQ=255)                                                                                                                           aTGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGa                           <  2:257495/98‑1 (MQ=255)                                                                                                                           aTGCACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGa                           >  1:257495/1‑98 (MQ=255)                                                                                                                              cACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGt                                                   >  1:274273/1‑71 (MQ=255)                                                                                                                              cACGCCGCCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGt                                                   <  2:274273/71‑1 (MQ=255)                                                                                                                                cgccgcCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCAt                                                                        >  2:790921/1‑48 (MQ=39)                                                                                                                                cgccgcCCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCAt                                                                        <  1:790921/48‑1 (MQ=39)                                                                                                                                     |                                                                                                                 GCGTGATCCTGGGTGTAGTGGCGTGAGGTAGTCCTGCGGCAGGAATGGGAAGATGGGCGGGCATGGCGGCTGCCCCCGGCGTGCAATAATAGCCACCCAAAAAAACAGCCGTTCCACATCCAATGCACGCCG‑CCTGGCCGGCTTCTCCCTCGGTCTTTCCCTGATTCTTGCCATCGGTTCCCAGAATGCTTTCGTGCTGCGACAAGGGTTGCGGCGCGAGCATGTGTTCTGGGTCTGCCTGGTG  >  AM260479/724635‑724878

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