Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	589,509	G→A	A482A (GCC→GCT)	nfrA ←	bacteriophage N4 receptor, outer membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	589,509	0	G	A	100.0%	113.6 / NA	33	A482A (GCC→GCT)	nfrA	bacteriophage N4 receptor, outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/16);  total (17/16)

GTTGTTCGGCCTGAAGCCATGCATACAACGCCACACCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAGGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCGGGCAATTATCTGCAATACCCGGC  >  NC_000913/589433‑589584                                                                             |                                                                            gttgttCGGCCTGAAGCCATGCATACAACGCCACACCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAagcgg                                                                         <  1:111861/81‑1 (MQ=255) gttgttCGGCCTGAAGCCATGCATACAACGCCACACCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAagcgg                                                                         >  2:111907/1‑81 (MQ=255) gttgttCGGCCTGAAGCCATGCATACAACGCCACACCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAagcg                                                                          <  2:111862/80‑1 (MQ=255)    gttCGGCCTGAAGCCATGCATACAACGCCACACCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAagcggcgg                                                                      >  1:111931/1‑81 (MQ=255)        ggCCTGAAGCCATGCATACAACGCCACACCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCAt                                                                   <  1:111864/80‑1 (MQ=255)        ggCCTGAAGCCATGCATACAACGCCACACCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATc                                                                  <  1:111863/81‑1 (MQ=255)        ggCCTGAAGCCATGCATACAACGCCACACCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATc                                                                  >  2:111924/1‑81 (MQ=255)         gCCTGAAGCCATGCATACAACGCCACACCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCg                                                                 >  1:111920/1‑81 (MQ=255)          ccTGAAGCCATGCATACAACGCCACACCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCGt                                                                <  1:111866/81‑1 (MQ=255)          ccTGAAGCCATGCATACAACGCCACACCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCGt                                                                <  2:111865/81‑1 (MQ=255)                    tGCATACAACGCCACACCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGa                                                      >  2:111922/1‑81 (MQ=255)                      cATACAACGCCACACCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGAc                                                     <  1:111867/80‑1 (MQ=255)                       aTACAACGCCACACCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACat                                                   <  1:111868/81‑1 (MQ=255)                        tACAACGCCACACCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACata                                                  >  2:111914/1‑81 (MQ=255)                              gCCACACCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGcc                                             >  1:111913/1‑80 (MQ=255)                                 acacCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAg                                          >  2:111906/1‑80 (MQ=255)                                 acacCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGc                                         <  1:111869/81‑1 (MQ=255)                                  cacCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCa                                        <  2:111870/81‑1 (MQ=255)                                     cGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGc                                     <  1:111872/81‑1 (MQ=255)                                     cGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGc                                     <  2:111871/81‑1 (MQ=255)                                             gtgtCCCGATAACACTTTGCCAGACGGATCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTAtt                             >  1:111938/1‑81 (MQ=255)                                                  ccGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCgg                        >  1:111934/1‑81 (MQ=255)                                                   cGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCggg                       <  2:111873/81‑1 (MQ=255)                                                   cGATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCggg                       >  2:117393/1‑81 (MQ=255)                                                    gATAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCggg                       >  1:111921/1‑80 (MQ=255)                                                      tAACACTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCGGGCaa                    <  1:111874/81‑1 (MQ=255)                                                       aaCACTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCGGGCAAt                   >  1:111927/1‑81 (MQ=255)                                                        acacTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCGGGCAAtt                  <  2:111875/81‑1 (MQ=255)                                                        acacTTTGCCAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCGGGCAAt                   <  1:111876/80‑1 (MQ=255)                                                                ccAGACGGTTCCAAGCGGTGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCGGGCAATTATCTGCaa          >  2:111942/1‑81 (MQ=255)                                                                ccAGACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCGGGCAATTATCTGCaa          >  2:111952/1‑81 (MQ=255)                                                                 cAGACGGTTCCAAACGGCAGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCGGGCAATTATCTGCAAt         <  2:111877/81‑1 (MQ=255)                                                                   gACGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCGGGCAATTATCTGCAATAc       >  1:111915/1‑81 (MQ=255)                                                                     cGGTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCGGGCAATTATCTGCAATAccc     >  2:111929/1‑81 (MQ=255)                                                                      ggTTCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCGGGCAATTATCTGCAATACCCg    <  1:111878/81‑1 (MQ=255)                                                                        ttCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCGGGCAATTATCTGCAATACCCgg   >  1:111930/1‑80 (MQ=255)                                                                         tCCAAGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCGGGCAATTATCTGCAATACCCGGc  >  1:111941/1‑80 (MQ=255)                                                                             |                                                                            GTTGTTCGGCCTGAAGCCATGCATACAACGCCACACCGGGTAGCGTGTCCCGATAACACTTTGCCAGACGGTTCCAGGCGGCGGCATCGTAGGAAGGCGACATATCGCCCAGCAAGCGAACTATTGCCGGGCAATTATCTGCAATACCCGGC  >  NC_000913/589433‑589584

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