Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	591,410	T→G	D590A (GAC→GCC)	nfrB ←	bacteriophage N4 receptor, inner membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	591,410	0	T	G	100.0%	129.1 / NA	38	D590A (GAC→GCC)	nfrB	bacteriophage N4 receptor, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (17/21);  total (17/21)

AACATTGCCCGCGGATCGTGACCGCGTCGGCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGG  >  NC_000913/591275‑591530                                                                                                                                        |                                                                                                                         aacaTTGCCCGCGGATCGTGACCGCGTCGGCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAAt                                                                                                                      <  2:751358/140‑1 (MQ=255) aacaTTGCCCGCGGATCGTGACCGCGTCGGCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAAt                                                                                                                      <  2:416298/140‑1 (MQ=255)      tGCCCGCGGATCGTGACCGCGTCGGCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCAt                                                                                                                 >  2:659044/1‑140 (MQ=255)        cccGCGGATCGTGACCGCGTCGGCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAAt                                                                                                                      <  1:684737/133‑1 (MQ=255)        cccGCGGATCGTGACCGCGTCGGCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAAt                                                                                                                      >  2:684737/1‑133 (MQ=255)        cccGCGGATCGTGACCGCGTCGGCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCt                                                                                                               >  1:431299/1‑140 (MQ=255)             ggATCGTGACCGCGTCGGCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACt                                                                                                          <  1:795352/140‑1 (MQ=255)                       cgcgTCGGCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGcc                                                                                                       <  2:523456/133‑1 (MQ=255)                       cgcgTCGGCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGcc                                                                                                       >  1:523456/1‑133 (MQ=255)                              gCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGt                                                                                         <  2:431299/140‑1 (MQ=255)                                         aGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGa                                                                              >  2:632511/1‑140 (MQ=255)                                           tGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACg                                                                            >  1:600061/1‑140 (MQ=255)                                                 aaCCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCgg                                                                      <  2:54665/140‑1 (MQ=255)                                                  aCCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGc                                                                     >  1:177007/1‑140 (MQ=255)                                                   cccTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCa                                                                    <  2:177007/140‑1 (MQ=255)                                                      tgtgACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTa                                                                 <  1:659044/140‑1 (MQ=255)                                                      tgtgACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTa                                                                 <  2:755709/140‑1 (MQ=255)                                                         gACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTa                                                                 <  2:666679/137‑1 (MQ=255)                                                         gACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTa                                                                 >  1:666679/1‑137 (MQ=255)                                                              tttGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATaa                                                         >  1:39269/1‑140 (MQ=255)                                                                tGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATg                                                       <  2:525305/140‑1 (MQ=255)                                                                     ccGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAgcg                                                  <  2:309212/140‑1 (MQ=255)                                                                               aTGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGc                                        >  2:208229/1‑140 (MQ=255)                                                                               aTGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGc                                        >  2:331879/1‑140 (MQ=255)                                                                                tGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGcc                                       <  1:331879/140‑1 (MQ=255)                                                                                tGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGcc                                       >  1:212777/1‑140 (MQ=255)                                                                                tGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGcc                                       <  1:208229/140‑1 (MQ=255)                                                                                 gACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCg                                      <  2:212777/140‑1 (MQ=255)                                                                                        cGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCCGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTc                               >  1:520191/1‑140 (MQ=255)                                                                                              ttGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAAt                         <  1:283542/140‑1 (MQ=255)                                                                                                   tCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAACACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCa                       <  1:818081/137‑1 (MQ=255)                                                                                                   tCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAACACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCa                       >  2:818081/1‑137 (MQ=255)                                                                                                       aCTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGa                >  2:759997/1‑140 (MQ=255)                                                                                                        cTTTACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGaa               <  1:759997/140‑1 (MQ=255)                                                                                                           tACGCGTCAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAggg            <  2:346980/140‑1 (MQ=255)                                                                                                                  cAGGGCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGcc     >  1:260574/1‑140 (MQ=255)                                                                                                                     ggCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAgg  >  2:40452/1‑140 (MQ=255)                                                                                                                     ggCCGCCAGCGAAACCGGGGCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAgg  <  2:260574/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                         AACATTGCCCGCGGATCGTGACCGCGTCGGCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGG  >  NC_000913/591275‑591530

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: ‑