Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,768,876	G→A	R71Q (CGG→CAG)	yfjQ →	CP4‑57 prophage, uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,768,876	0	G	A	100.0%	135.2 / NA	39	R71Q (CGG→CAG)	yfjQ	CP4‑57 prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/22);  total (17/22)

GGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTA  >  NC_000913/2768748‑2769002                                                                                                                                 |                                                                                                                               ggtgCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGAATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGAc                                                                                                                      <  1:429282/139‑1 (MQ=255)          cGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCcgtcgt                                                                                                             <  2:343923/139‑1 (MQ=255)          cGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAACCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCcgtcgt                                                                                                             <  2:103555/139‑1 (MQ=255)            tATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCcgtcgtcg                                                                                                           >  1:442484/1‑139 (MQ=255)               tCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGa                                                                                                        <  2:728857/139‑1 (MQ=255)                cAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAg                                                                                                       <  1:403821/139‑1 (MQ=255)                cAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAg                                                                                                       <  2:64786/139‑1 (MQ=255)                cAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAg                                                                                                       <  2:416308/139‑1 (MQ=255)                   tGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTg                                                                                                         <  1:301246/134‑1 (MQ=255)                   tGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTg                                                                                                         >  2:301246/1‑134 (MQ=255)                      ggATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACt                                                                                                 <  2:442484/139‑1 (MQ=255)                      ggATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACt                                                                                                 >  1:162677/1‑139 (MQ=255)                                 aaTCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATAtgct                                                                                      >  1:5488/1‑139 (MQ=255)                                 aaTCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATAtgct                                                                                      >  2:592340/1‑139 (MQ=255)                                     ccGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTgcgc                                                                                  <  1:286002/139‑1 (MQ=255)                                       gCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGa                                                                                                        <  1:437170/115‑1 (MQ=255)                                       gCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGa                                                                                                        >  2:437170/1‑115 (MQ=255)                                                       gTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                                                >  1:801754/1‑139 (MQ=255)                                                         atataCCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAt                                                              <  2:801754/139‑1 (MQ=255)                                                            taCCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAc                                                           >  2:411599/1‑139 (MQ=255)                                                               ccACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGa                                                                                                        >  2:336493/1‑91 (MQ=255)                                                               ccACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGa                                                                                                        <  1:336493/91‑1 (MQ=255)                                                                 acaATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTg                                                                                                         >  2:677201/1‑88 (MQ=255)                                                                 acaATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTg                                                                                                         <  1:677201/88‑1 (MQ=255)                                                                 acaATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGa                                                                                                        <  1:762904/89‑1 (MQ=255)                                                                 acaATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGa                                                                                                        >  2:762904/1‑89 (MQ=255)                                                                 acaATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                      >  1:66647/1‑139 (MQ=255)                                                                  caATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                     >  1:781207/1‑139 (MQ=255)                                                                                 aCAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAt                                      >  1:296486/1‑139 (MQ=255)                                                                                    gCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                   <  2:296486/139‑1 (MQ=255)                                                                                    gCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                   <  2:671530/139‑1 (MQ=255)                                                                                              ggAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATAc                                                                                                  <  2:216041/66‑1 (MQ=255)                                                                                              ggAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATAc                                                                                                  >  1:216041/1‑66 (MQ=255)                                                                                                aaGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAc                                                           <  2:633954/103‑1 (MQ=255)                                                                                                aaGGCTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAc                                                           >  1:633954/1‑103 (MQ=255)                                                                                                    cTTCCAGCCATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                   <  2:162677/139‑1 (MQ=255)                                                                                                            cATTCTTTGCCTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCAc           <  2:556373/139‑1 (MQ=255)                                                                                                                     ccTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa  <  2:5488/139‑1 (MQ=255)                                                                                                                     ccTGTCAGACACAGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa  >  2:341374/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                               GGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTA  >  NC_000913/2768748‑2769002

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