Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,274,178	G→A	*605* (TAG→TAA)	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,178	0	G	A	100.0%	102.6 / NA	33	*605* (TAG→TAA)	yejA	microcin C ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/16);  total (17/16)

GGAAGACCGTCTCGCCTGGTGGGATAAATTCTCCCAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAACAGGGAGAGTAGATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCAACTTTTTCATCGTGCAAATTGCGCCTGGCGGTCCGGTCGACCAGGCCATCGCCGCCA  >  NC_000913/2274043‑2274302                                                                                                                                        |                                                                                                                             ggAAGACCGTCTCGCCTGGTGGGATAAATTCTCCCAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATg                                                                                                                           <  1:925826/139‑1 (MQ=255)          tctcGCCTGGTGGGATAAATTCTCCCAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTAc                                                                                                                  >  1:883127/1‑139 (MQ=255)                     gggATAAATTCTCCCAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCg                                                                                                       <  2:463270/139‑1 (MQ=255)                      ggATAAATTCTCCCAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGt                                                                                                      <  2:883127/139‑1 (MQ=255)                      ggATAAATTCTCCCAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGt                                                                                                      <  2:310940/139‑1 (MQ=255)                                 cccAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGt                                                                                           <  1:103754/139‑1 (MQ=255)                                             tGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGCAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACAtt                                                                               >  2:196929/1‑139 (MQ=255)                                              gcgcCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATT‑GCCGTCTGTTGCTGGTGATCCCAACATTAt                                                                             <  2:574896/139‑1 (MQ=255)                                                 ccccATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATgg                                                                           <  2:708367/139‑1 (MQ=255)                                                            gCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAAGTAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCAc                                                                >  1:589819/1‑139 (MQ=255)                                                              cTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCa                                                              >  2:133901/1‑139 (MQ=255)                                                                cGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATc                                                            >  2:640424/1‑139 (MQ=255)                                                                 ggTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCa                                                           >  1:546855/1‑139 (MQ=255)                                                                 ggTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCa                                                           >  2:98226/1‑139 (MQ=255)                                                                    aTCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCAACt                                                        >  2:965729/1‑139 (MQ=255)                                                                           ccTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCAACTTTTTCAt                                                 <  2:546855/139‑1 (MQ=255)                                                                            cTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGcc                                                                                                        >  2:388312/1‑83 (MQ=255)                                                                            cTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGcc                                                                                                        <  1:388312/83‑1 (MQ=255)                                                                             tgatggtATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCCAGCAGGAAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCAACTTTTTCATCg                                               >  1:584985/4‑139 (MQ=255)                                                                               gtggtATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGt                                                                                                      >  1:767632/1‑82 (MQ=255)                                                                               gtggtATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGt                                                                                                      <  2:767632/82‑1 (MQ=255)                                                                               gtggtATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCAACTTTTTCATCGTg                                             >  2:729328/1‑139 (MQ=255)                                                                                   tATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTAt                                                                             <  1:957003/103‑1 (MQ=255)                                                                                   tATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTAt                                                                             >  2:957003/1‑103 (MQ=255)                                                                                                 aGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTg                                                                                                >  1:923574/1‑70 (MQ=255)                                                                                                 aGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTg                                                                                                <  2:923574/70‑1 (MQ=255)                                                                                                 aGCGGCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCCTCTGTTCATGGTGATCCCAACATTATGGGCGATTATCACCATCAACTTTTTCATCGTGCAAATTGCGCCTGGcggt                           >  1:130990/1‑139 (MQ=255)                                                                                                     gCCAAACTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCAACTTTTTCATCGTGCAAATTGCGCCTGGcggtccgg                       <  2:806185/139‑1 (MQ=255)                                                                                                          aCTGCCGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCAACTTTTTCATCGTGCAAATTGCGCCTGGCGGTCCGGTCGAc                  <  2:584985/139‑1 (MQ=255)                                                                                                               cGTCCGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCAACTTTTTCATCGTGCAAATTGCGCCTGGCGGTCCGGTCGACCAGGc             >  1:369134/1‑139 (MQ=255)                                                                                                                   cGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATccc                                                                                     >  2:598112/1‑63 (MQ=255)                                                                                                                   cGCCAGCAAGCAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATccc                                                                                     <  1:598112/63‑1 (MQ=255)                                                                                                                          aagcAGGGAGAGTAAATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCAACTTTTTCATCGTGCAAATTGCGCCTGGCGGTCCGGTCGACCAGGCCATCGCCGCCa  <  2:369134/136‑1 (MQ=255)                                                                                                                                        |                                                                                                                             GGAAGACCGTCTCGCCTGGTGGGATAAATTCTCCCAGCCGGCCGTGCGCCCCATCTATAGCCTCGGTATCGATACCTGGTGGTATGACGTCAATAAAGCGGCCAAACTGCCGTCCGCCAGCAAACAGGGAGAGTAGATGGGCGCTTACCTGATTCGCCGTCTGTTGCTGGTGATCCCAACATTATGGGCGATTATCACCATCAACTTTTTCATCGTGCAAATTGCGCCTGGCGGTCCGGTCGACCAGGCCATCGCCGCCA  >  NC_000913/2274043‑2274302

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