Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	591,428	C→A	G584V (GGC→GTC)	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,428	0	C	A	100.0%	129.5 / NA	40	G584V (GGC→GTC)	nfrB	bacteriophage N4 receptor, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (18/22);  total (18/22)

GTGACCGCGTCGGCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTT  >  NC_000913/591292‑591558                                                                                                                                         |                                                                                                                                   gTGACCGCGTCGGCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGCACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGa                                                                                                                                 <  2:445158/140‑1 (MQ=255)   gACCGCGTCGGCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACa                                                                                                                               >  2:396179/1‑140 (MQ=255)        cgTCGGCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTaa                                                                                                                          <  1:11818/140‑1 (MQ=255)          tCGGCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTGAATCCCATCTTCACTGACGACAATTAACt                                                                                                                        >  1:272257/1‑140 (MQ=255)                           gACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGc                                                                                                       <  2:623588/140‑1 (MQ=255)                            aCGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCa                                                                                                      >  1:520847/1‑140 (MQ=255)                               tAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATCCCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCg                                                                                                   >  1:736332/1‑140 (MQ=255)                               tAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCg                                                                                                   >  2:110452/1‑140 (MQ=255)                                aaCCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCgg                                                                                                  <  2:23504/140‑1 (MQ=255)                                aaCCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCgg                                                                                                  <  1:110452/140‑1 (MQ=255)                                aaCCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCgg                                                                                                  <  2:482209/140‑1 (MQ=255)                                aaCCCTGTGACAATTTGTCCCCCCAGAACAAAGACCTAACGCACTTTGTGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCgg                                                                                                  <  2:272257/140‑1 (MQ=255)                                aaCCCTGTGACAATTTGGCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACCTTACGCGACAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCgg                                                                                                  <  2:136990/140‑1 (MQ=255)                                 aCCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGc                                                                                                 >  1:719279/1‑140 (MQ=255)                                 aCCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGc                                                                                                 >  2:139064/1‑140 (MQ=255)                                     tgtgACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTa                                                                                             <  1:120725/140‑1 (MQ=255)                                     tgtgACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTa                                                                                             <  2:19033/140‑1 (MQ=255)                                     tgtgACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGAAACTGGACTCGTCATTTTTAAGACGCATCGGCAGTa                                                                                             >  2:351442/1‑140 (MQ=255)                                      gtgACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTAc                                                                                            >  2:59557/1‑140 (MQ=255)                                        gACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCg                                                                                          >  2:370391/1‑140 (MQ=255)                                        gACAATTTGTCCCCGAAGAACAATGACGTAACGCACTTAGCGTCCGACTTTACGCGTCAGCGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGGCATTTTCCAGACGCAGCGGCAGTACCg                                                                                          >  1:136990/1‑140 (MQ=255)                                                               tGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCCTTTTCCAGTCGCGGCGGCAGTACCGAATAATGCAGCGCCACGGAGGcc                                                                   >  2:287499/1‑140 (MQ=255)                                                                   gTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCa                                                               <  2:110860/140‑1 (MQ=255)                                                                    tAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCAt                                                              >  2:184623/1‑140 (MQ=255)                                                                      aCGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCAttt                                                            >  2:529963/1‑140 (MQ=255)                                                                             ttGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAAt                                                     >  1:122294/1‑140 (MQ=255)                                                                              tGCGTCCGCCTTTACGTGTCAGCGCCGCCATGGAACCCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATc                                                    <  2:342368/140‑1 (MQ=255)                                                                              tGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGAATAATGCAGCGTCATGGAGGCAGGCATTTCGGCAATc                                                    >  1:342368/1‑140 (MQ=255)                                                                                cGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAg                                                  <  2:122294/140‑1 (MQ=255)                                                                                 gTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGc                                                 <  2:595955/140‑1 (MQ=255)                                                                                  tCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCg                                                <  2:520847/140‑1 (MQ=255)                                                                                    cGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAg                                              >  2:599731/1‑140 (MQ=255)                                                                                    cGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGAACCGCATAATGCGGCGACACGGAGGCCGGCCTTTCGGCAATCAGCGAg                                              >  2:281046/1‑140 (MQ=255)                                                                                       cTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGaa                                           <  1:599499/140‑1 (MQ=255)                                                                                          tACGCGTCAGGGCCGCCAGCGATACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAggg                                        <  1:281046/140‑1 (MQ=255)                                                                                          tACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAggg                                        <  1:599731/140‑1 (MQ=255)                                                                                                tCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGc                                  <  1:396179/140‑1 (MQ=255)                                                                                                       cgccAGCGAAACCGGGTGAATACCATCTTCACTGACGACCATTAACTTGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCAt                           <  2:253921/140‑1 (MQ=255)                                                                                                                           tACCATCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCccacgccacg       <  2:621577/140‑1 (MQ=255)                                                                                                                                tCTTCACTGACGACATTTAACTCGTAATTCTCCAGACTCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGtt  <  1:351442/140‑1 (MQ=255)                                                                                                                                tCTTCACTGACGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGtt  >  1:520605/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                   GTGACCGCGTCGGCGTGCATACCAGTGACGTAACCCTGTGACAATTTGTCCCCGCAGAACAATGACGTAACGCACTTTGCGTCCGACTTTACGCGTCAGGGCCGCCAGCGAAACCGGGTCAATACCATCTTCACTGCCGACAATTAACTCGTCATTTTCCAGACGCAGCGGCAGTACCGCATAATGCAGCGCCACGGAGGCCGGCATTTCGGCAATCAGCGAGGAAGGGATCTGCCAGGCATCGATGGATTCCCACGCCACGCCGTT  >  NC_000913/591292‑591558

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