Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	588,928	A→G	L676P (CTT→CCT)	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	588,928	0	A	G	100.0%	112.4 / NA	33	L676P (CTT→CCT)	nfrA	bacteriophage N4 receptor, outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/14);  total (19/14)

AGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAAGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCCACAAGGCGTAACCAAGCGCTGCCTGGGTGTTGCTATTATTCGGTTCCAGTTCCAGCGCGGCGCGCAAATCACTCAC  >  NC_000913/588799‑589061                                                                                                                                  |                                                                                                                                      aGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATg                                                                                                                             >  1:594739/1‑140 (MQ=255)    ggggTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATgcgc                                                                                                                          <  2:594739/140‑1 (MQ=255)             aGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAg                                                                                                                 >  1:159944/1‑140 (MQ=255)             aGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAg                                                                                                                 >  2:184415/1‑140 (MQ=255)             aGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCATCCTGACGCTCTTGATGCGCCGGTTCGAg                                                                                                                 >  1:338239/1‑140 (MQ=255)                gcCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCa                                                                                                               >  2:430204/1‑139 (MQ=255)                gcCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCa                                                                                                               <  1:430204/139‑1 (MQ=255)                     aTTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAg                                                                                                                 >  1:745685/1‑132 (MQ=255)                     aTTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAg                                                                                                                 <  2:745685/132‑1 (MQ=255)                            aTGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACtg                                                                                                  <  2:159944/140‑1 (MQ=255)                               tcatcaATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGc                                                                                               >  1:608461/1‑140 (MQ=255)                                catcaATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGa                                                                                                     <  2:784135/133‑1 (MQ=255)                                catcaATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGa                                                                                                     >  1:784135/1‑133 (MQ=255)                                catcaATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCg                                                                                              >  1:11230/1‑140 (MQ=255)                                         ccAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCg                                                                                     >  1:171997/1‑140 (MQ=255)                                         ccAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCg                                                                                     <  1:184415/140‑1 (MQ=255)                                          cAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCAGGTTCGAGCATCTCCCGCGACAGTGCGATATCACCGc                                                                                    >  1:272348/1‑140 (MQ=255)                                            gCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCAGGTGCGAACATTTCCCCCGACGGTGCGATATCACCGCTa                                                                                  >  2:246968/1‑140 (MQ=255)                                                 gCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCca                                                                             >  1:467187/1‑140 (MQ=255)                                                  cGTAGTGCTGCGTCGCAGGCATCTCCTCCAGACGCTGATTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCcac                                                                            <  1:246968/140‑1 (MQ=255)                                                   gTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCcaca                                                                           <  2:608461/140‑1 (MQ=255)                                                       tgctgcGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCCACAAGGc                                                                       >  2:45814/1‑140 (MQ=255)                                                                        gTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCCACAAGGCGTAACCAAGCGCTGCCt                                                      <  2:11230/140‑1 (MQ=255)                                                                              cAGACGCTGGTTCACGTGGGCCTGTTGTCGGCTCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCCACAAGGCGTAACCAAGCGCTGCCTGGgtgt                                                <  2:272348/140‑1 (MQ=255)                                                                              cAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCCACAAGGCGTAACCAAGCGCTGCCTGGgtgt                                                <  1:743188/140‑1 (MQ=255)                                                                                       gTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCCACAAGGCGTAACCAAGCGCTGCCTGGGTGTTGCtattat                                       >  1:682396/1‑140 (MQ=255)                                                                                        ttCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCCACAAGGCGTAACCAAGCGCTGCCTGGGTGTTGCTattatt                                      >  2:692636/1‑140 (MQ=255)                                                                                            cGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCCACAAGGCGTAACCAAGCGCTGCCTGGGTGTTGCTATTATTCGgt                                  <  1:426051/140‑1 (MQ=255)                                                                                                         tCGGATCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCCACAAGGCGTAACCAAGCGCTGCCTGGGTGTTGCTATTATTCGGTTCCAGTTCCAgcg                     >  1:700947/1‑140 (MQ=255)                                                                                                             aTCAGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCCACAAGGCGTAACCAAGCGCTGCCTGGGTGTTGCTATTATTCGGTTCCAGTTCCAgcgcggc                 <  2:682396/140‑1 (MQ=255)                                                                                                                aGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCCACAAGGCGTAACCAAGCGCTGCCTGGGTGTTGCTATTATTCGGTTCCAGTTCCAgcgcggcgc               <  2:185134/139‑1 (MQ=255)                                                                                                                aGTGCCGGATCGTCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCCACAAGGCGTAACCAAGCGCTGCCTGGGTGTTGCTATTATTCGGTTCCAGTTCCAgcgcggcgc               >  1:185134/1‑139 (MQ=255)                                                                                                                            tCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCCACAAGGCGTAACCAAGCGCTGCCTGGGTGTTGCTATTATTCGGTTCCAGTTCCAGCGCGGCGCGCAAAtcactcac  >  2:3976/1‑140 (MQ=255)                                                                                                                            tCCGGAGGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCCACAAGGCGTAACCAAGCGCTGCCTGGGTGTTGCTATTATTCGGTTCCAGTTCCAGCGCGGCGCGCAAAtcactcac  >  1:558233/1‑140 (MQ=255)                                                                                                                                  |                                                                                                                                      AGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACGCTGGTTCACGTAGGCCAGTTGTCGGATCAGTGCCGGATCGTCCGGAAGCCCTTTATGCGCCGGTTCGAGCATTTCCCGCGACTGTGCGATATCACCGCTATCCCACAAGGCGTAACCAAGCGCTGCCTGGGTGTTGCTATTATTCGGTTCCAGTTCCAGCGCGGCGCGCAAATCACTCAC  >  NC_000913/588799‑589061

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