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%	144.0 / NA	45	*274* (TAG→TAA)	yfcO	DUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/24);  total (21/24)

CGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTGACCGAGACTT  >  NC_000913/2449094‑2449361                                                                                                                                       |                                                                                                                                      cGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACa                                                                                                                                   >  2:1411510/1‑139 (MQ=255) cGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACa                                                                                                                                   >  2:1502586/1‑139 (MQ=255)    tGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACaa                                                                                                                                  <  2:447168/137‑1 (MQ=255)    tGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACaa                                                                                                                                  >  1:447168/1‑137 (MQ=255)    tGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACt                                                                                                                                >  1:1419816/1‑139 (MQ=255)        gcggcgTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCAtt                                                                                                                            <  2:1419816/139‑1 (MQ=255)            cgTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGt                                                                                                                        <  1:1411510/139‑1 (MQ=255)            cgTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGt                                                                                                                        <  1:218033/139‑1 (MQ=255)             gttgttGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTa                                                                                                                       <  2:5997/139‑1 (MQ=255)             gttgttGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTa                                                                                                                       <  2:317324/139‑1 (MQ=255)                   ggTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACt                                                                                                                 >  1:1269661/1‑139 (MQ=255)                      gtTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGa                                                                                                              >  2:1006770/1‑139 (MQ=255)                       ttgattgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGaa                                                                                                             <  2:1269661/139‑1 (MQ=255)                       ttgattgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGaa                                                                                                             <  1:1502586/139‑1 (MQ=255)                          attgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCg                                                                                                          <  2:1047874/139‑1 (MQ=255)                                         tGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGAt                                                                                           <  1:1006770/139‑1 (MQ=255)                                           gCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCa                                                                                                                              >  1:1381619/1‑102 (MQ=255)                                           gCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCa                                                                                                                              <  2:1381619/102‑1 (MQ=255)                                                                  gagTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAaca                                                                  <  1:570904/139‑1 (MQ=255)                                                                     ttGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGAt                                                                                           >  1:505372/1‑111 (MQ=255)                                                                     ttGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGAt                                                                                           <  2:505372/111‑1 (MQ=255)                                                                      tGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCg                                                                                                          <  1:793317/95‑1 (MQ=255)                                                                      tGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCg                                                                                                          >  2:793317/1‑95 (MQ=255)                                                                        aGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCg                                                                                                          <  1:1283512/93‑1 (MQ=255)                                                                        aGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCg                                                                                                          >  2:1283512/1‑93 (MQ=255)                                                                          cgcgTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCAtt                                                                                                                            >  1:392383/1‑73 (MQ=255)                                                                          cgcgTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCAtt                                                                                                                            <  2:392383/73‑1 (MQ=255)                                                                          cgcgTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGAc                                                          >  2:1603488/1‑139 (MQ=255)                                                                              ttGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCata                                                      <  1:1603488/139‑1 (MQ=255)                                                                                            ggTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgca                                                              >  1:1424849/1‑117 (MQ=255)                                                                                            ggTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgca                                                              <  2:1424849/117‑1 (MQ=255)                                                                                                aTGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAAc                                                                                                                                 >  2:874506/1‑46 (MQ=255)                                                                                                aTGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAAc                                                                                                                                 <  1:874506/46‑1 (MQ=255)                                                                                                aTGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGt                                      <  2:126445/137‑1 (MQ=255)                                                                                                aTGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGt                                      >  1:126445/1‑137 (MQ=255)                                                                                                 tGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGt                                                                                                                        <  2:1085870/54‑1 (MQ=255)                                                                                                 tGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGt                                                                                                                        >  1:1085870/1‑54 (MQ=255)                                                                                                       gggCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGaa                                                                                                             <  2:1249026/59‑1 (MQ=255)                                                                                                       gggCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGaa                                                                                                             >  1:1249026/1‑59 (MQ=255)                                                                                                       gggCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacg                                                                >  1:763217/1‑104 (MQ=255)                                                                                                       gggCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacg                                                                <  2:763217/104‑1 (MQ=255)                                                                                                       gggCGAATGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTAGTCCGTTATAg                             >  2:1404162/1‑139 (MQ=255)                                                                                                            aaTGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacg                                                                <  2:316428/99‑1 (MQ=255)                                                                                                            aaTGCCCCTGCCTGGTGTATCTAAGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacg                                                                >  1:316428/1‑99 (MQ=255)                                                                                                                                  aaGGGTTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTGACCGAGACtt  >  2:1150059/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                      CGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTGACCGAGACTT  >  NC_000913/2449094‑2449361

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