Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,449,228	C→T	*274* (TAG→TAA)	yfcO ←	DUF2544 family putative outer membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,449,228	0	C	T	100.0%	119.0 / NA	38	*274* (TAG→TAA)	yfcO	DUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/19);  total (19/19)

TGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTGACCGAGACTTCC  >  NC_000913/2449100‑2449363                                                                                                                                 |                                                                                                                                        tGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCAt                                                                                                                               <  2:983588/139‑1 (MQ=255)     gcgTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAg                                                                                                                           >  1:1374538/1‑139 (MQ=255)                                      cTGCCGCCGGAGTTGCCCTGAGGAGTTGAGTGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTgg                                                                                          <  1:1131397/139‑1 (MQ=255)                                        gccgccGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCt                                                                                        <  2:1374538/139‑1 (MQ=255)                                           gccgGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCa                                                                                     <  1:134708/139‑1 (MQ=255)                                                    gCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCTTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGAt                                                                                             <  1:548349/122‑1 (MQ=255)                                                    gCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCTTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGAt                                                                                             >  2:548349/1‑122 (MQ=255)                                                          aggagTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAAACGATTGGCTGCATACTGCTGGTATCaa                                                                      >  2:279065/1‑139 (MQ=255)                                                              gTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAAACGATTGGCTGCATACTGCTGGTATCAAcacg                                                                  <  1:279065/139‑1 (MQ=255)                                                                tGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgca                                                                <  2:256509/139‑1 (MQ=255)                                                                       gTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCat                                                         >  1:822702/1‑139 (MQ=255)                                                                                ggTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTaa                                                >  1:777483/1‑139 (MQ=255)                                                                                     gggTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATActg                                                                                <  2:370514/102‑1 (MQ=255)                                                                                     gggTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATActg                                                                                >  1:370514/1‑102 (MQ=255)                                                                                      ggTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGc                                                                                                          <  1:1071958/75‑1 (MQ=255)                                                                                      ggTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGc                                                                                                          >  2:1071958/1‑75 (MQ=255)                                                                                          aTGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAg                                                                                                    <  1:609838/77‑1 (MQ=255)                                                                                          aTGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAg                                                                                                    >  2:609838/1‑77 (MQ=255)                                                                                           tGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgca                                                                >  2:508728/1‑112 (MQ=255)                                                                                           tGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgca                                                                <  1:508728/112‑1 (MQ=255)                                                                                               aGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGAt                                                                                             >  2:570472/1‑79 (MQ=255)                                                                                               aGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGAt                                                                                             <  1:570472/79‑1 (MQ=255)                                                                                                        tGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCAt                                                                                                                               >  1:248422/1‑36 (MQ=255)                                                                                                        tGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCAt                                                                                                                               <  2:248422/36‑1 (MQ=255)                                                                                                                cctgGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATAc                                                                                  <  2:176796/73‑1 (MQ=255)                                                                                                                cctgGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATAc                                                                                  >  1:176796/1‑73 (MQ=255)                                                                                                                    gtgtATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCgtgt                                                  <  2:1384609/101‑1 (MQ=255)                                                                                                                    gtgtATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCgtgt                                                  >  1:1384609/1‑101 (MQ=255)                                                                                                                          cTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGAt                                                                                             <  2:863782/52‑1 (MQ=255)                                                                                                                          cTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGAt                                                                                             >  1:863782/1‑52 (MQ=255)                                                                                                                          cTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATActgc                                                                               <  1:808640/66‑1 (MQ=255)                                                                                                                          cTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATActgc                                                                               >  2:808640/1‑66 (MQ=255)                                                                                                                          cTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTg              <  1:845956/131‑1 (MQ=255)                                                                                                                          cTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTg              >  2:845956/1‑131 (MQ=255)                                                                                                                          cTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTGACCgaga       <  2:886230/138‑1 (MQ=255)                                                                                                                          cTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTGACCgaga       >  1:886230/1‑138 (MQ=255)                                                                                                                          cTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTGACCGAGAc      >  1:565915/1‑139 (MQ=255)                                                                                                                              gggTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTGACCGAGACTTcc  >  1:1227590/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                        TGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTGACCGAGACTTCC  >  NC_000913/2449100‑2449363

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