Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,274,166	A→G	K601K (AAA→AAG)	yejA →	microcin C ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,274,166	0	A	G	100.0%	122.5 / NA	35	K601K (AAA→AAG)	yejA	microcin C ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/16);  total (19/16)

GTACATGGCGGAAGACCGTCTCGCCTGGTGGGATAAATTCTCCCAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAACAGGGAGAGTAGATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCAACTTTTTCATCGTGCAAATTGCGCCTGGCGGTCCGG  >  NC_000913/2274034‑2274281                                                                                                                                     |                                                                                                                    gTACATGGCGGAAGACCGTCTCGCCTGGTGGGATAAATTCTCCCAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAgg                                                                                                                 >  2:4607/1‑137 (MQ=255) gTACATGGCGGAAGACCGTCTCGCCTGGTGGGATAAATTCTCCCAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAgg                                                                                                                 <  1:4607/137‑1 (MQ=255)   aCATGGCGGAAGACCGTCTCGCCTGGTGGGATAAATTCTCCCAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGAAAGCa                                                                                                                   >  2:522881/1‑133 (MQ=255)       ggCGGAAGACCGTCTCGCCTGGTGGGATAAATTCTCCCAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTaa                                                                                                         >  2:27816/1‑138 (MQ=255)       ggCGGAAGACCGTCTCGCCTGGTGGGATAAATTCTCCCAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTaa                                                                                                         >  1:361852/1‑138 (MQ=255)       ggCGGAAGACCGTCTCGCCTGGTGGGATAAATTCTCCCAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTaa                                                                                                         >  2:582351/1‑138 (MQ=255)           gAAGACCGTCTCGCCTGGTGGGATAAATTCTCCCAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATgg                                                                                                     <  2:313090/139‑1 (MQ=255)              gACCGTCTCGCCTGGTGGGATAAATTCTCCCAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGgcg                                                                                                  <  2:361852/139‑1 (MQ=255)                                                  gccgTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCaa                                                              >  1:786281/1‑139 (MQ=255)                                                           cccATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATggg                                                     <  1:504452/139‑1 (MQ=255)                                                              aTCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAgg                                                                                                                 <  2:80453/76‑1 (MQ=255)                                                              aTCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAgg                                                                                                                 >  1:80453/1‑76 (MQ=255)                                                                         cGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTc                                                                                <  1:301166/98‑1 (MQ=255)                                                                         cGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTc                                                                                >  2:301166/1‑98 (MQ=255)                                                                         cGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTAt                                              <  1:730096/132‑1 (MQ=255)                                                                         cGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTAt                                              <  1:689404/132‑1 (MQ=255)                                                                         cGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTAt                                              >  2:689404/1‑132 (MQ=255)                                                                         cGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTAt                                              >  2:730096/1‑132 (MQ=255)                                                                         cGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATc                                       <  1:722278/139‑1 (MQ=255)                                                                          ggTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAg                                                                                                                  <  2:279703/63‑1 (MQ=255)                                                                          ggTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAg                                                                                                                  >  1:279703/1‑63 (MQ=255)                                                                          ggTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTaaat                                                                                                       <  2:593788/74‑4 (MQ=255)                                                                          ggTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTaaat                                                                                                       >  1:593788/1‑71 (MQ=255)                                                                          ggTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCa                                      >  1:26548/1‑139 (MQ=255)                                                                                 aTACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGga                                                                                                               <  2:293016/59‑1 (MQ=255)                                                                                 aTACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGga                                                                                                               >  1:293016/1‑59 (MQ=255)                                                                                      tggtggTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGga                                                                                                               >  2:192027/1‑54 (MQ=255)                                                                                      tggtggTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGga                                                                                                               <  1:192027/54‑1 (MQ=255)                                                                                            tATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCAACTTTTTCATCGTGCaaa                    >  2:472946/1‑139 (MQ=255)                                                                                              tGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATggg                                                     >  1:768057/1‑104 (MQ=255)                                                                                              tGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATggg                                                     <  2:768057/104‑1 (MQ=255)                                                                                                      aTAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCAACTTTTTCATCGTGCAAATTGCGCCTgg          <  2:735624/139‑1 (MQ=255)                                                                                                            cGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCAACTTTTTCATCTTGCAAATTGCGCCTGGcggtcc    >  2:693450/1‑139 (MQ=255)                                                                                                             ggCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCAACTTTTTCATCGTGCAAATTGCGCCTGGcggtccg   >  1:5411/1‑139 (MQ=255)                                                                                                              gCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCAACTTTTTCATCGTGCAAATTGCGCCTGGcggtccgg  <  2:26548/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                    GTACATGGCGGAAGACCGTCTCGCCTGGTGGGATAAATTCTCCCAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAACAGGGAGAGTAGATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCAACTTTTTCATCGTGCAAATTGCGCCTGGCGGTCCGG  >  NC_000913/2274034‑2274281

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