Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	4,400,432	C→A	N48K (AAC→AAA)	hfq →	global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	4,400,432	0	C	A	100.0%	168.2 / NA	48	N48K (AAC→AAA)	hfq	global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (24/24);  total (24/24)

AGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACA  >  NC_000913_3_bme_pgi/4400286‑4400575                                                                                                                                                   |                                                                                                                                                agaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACg                                                                                                                                              <  3:305025/150‑1 (MQ=255)  gaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGGACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACgg                                                                                                                                             <  2:523571/150‑1 (MQ=255)    attgCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTc                                                                                                                                           <  1:494089/147‑1 (MQ=255)         tAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCa                                                                                                                                      >  3:184509/1‑150 (MQ=255)         tAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCa                                                                                                                                      >  3:184475/1‑150 (MQ=255)                            gATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCAcgcg                                                                                                                   >  1:540781/1‑150 (MQ=255)                                                  gTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCAc                                                                                                                      >  2:97833/1‑125 (MQ=255)                                                  gTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCAc                                                                                                                      <  1:97833/125‑1 (MQ=255)                                                   tCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCg                                                                                            >  1:102205/1‑150 (MQ=255)                                                         aCGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGt                                                                                      <  1:2922/150‑1 (MQ=255)                                                              ttCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGttt                                                                                    >  2:454302/1‑147 (MQ=255)                                                              ttCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGttt                                                                                    <  1:454302/147‑1 (MQ=255)                                                              ttCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtctc                                                                                 >  3:336127/1‑150 (MQ=255)                                                                 cAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACg                                                                                                                                              <  4:75737/86‑1 (MQ=255)                                                                 cAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACg                                                                                                                                              >  3:75737/1‑86 (MQ=255)                                                                    tttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAaaac                                                                                                                                               <  1:511738/82‑1 (MQ=255)                                                                           ttatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtc                                                                                   >  3:212873/1‑135 (MQ=255)                                                                           ttatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtc                                                                                   <  4:212873/135‑1 (MQ=255)                                                                                   tGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                                                           <  2:475157/135‑1 (MQ=255)                                                                                   tGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                                                           >  1:475157/1‑135 (MQ=255)                                                                                     aaTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTa                                                                                                            <  1:340137/100‑1 (MQ=255)                                                                                     aaTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTa                                                                                                            >  2:340137/1‑100 (MQ=255)                                                                                     aaTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCg                                                                                        >  2:300322/1‑120 (MQ=255)                                                                                     aaTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCg                                                                                        <  1:300322/120‑1 (MQ=255)                                                                                     aaTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGt                                                             >  3:19980/1‑147 (MQ=255)                                                                                     aaTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGt                                                             <  4:19980/147‑1 (MQ=255)                                                                                      aTGGTATTAAGCTGCAAGGGCAAATAGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCg                                                         <  2:284978/150‑1 (MQ=255)                                                                                        ggTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGt                                                       >  1:483393/1‑150 (MQ=255)                                                                                             tAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcag                                                  <  2:102205/150‑1 (MQ=255)                                                                                                        ggCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTAcc                                       >  3:120651/1‑150 (MQ=255)                                                                                                         gCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaac                                                                      >  2:64387/1‑118 (MQ=255)                                                                                                         gCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaac                                                                      <  1:64387/118‑1 (MQ=255)                                                                                                           aaaTCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcag                                                  <  2:254070/136‑1 (MQ=255)                                                                                                           aaaTCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcag                                                  >  1:254070/1‑136 (MQ=255)                                                                                                               cGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGc                                                                                                                        >  4:207344/1‑62 (MQ=255)                                                                                                               cGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGc                                                                                                                        <  3:207344/62‑1 (MQ=255)                                                                                                                gAGTCTTTTGATAAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGt                               <  4:34250/150‑1 (MQ=255)                                                                                                                  gTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTag                             <  2:540781/150‑1 (MQ=255)                                                                                                                    cTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTagca                           >  1:505372/1‑150 (MQ=255)                                                                                                                      tttGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTagcagc                         >  2:114676/1‑150 (MQ=255)                                                                                                                               ttCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAAt                <  4:184475/150‑1 (MQ=255)                                                                                                                               ttCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAAt                <  4:184509/150‑1 (MQ=255)                                                                                                                                 cGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGAt                                                                                                                 <  2:332850/51‑1 (MQ=255)                                                                                                                                 cGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGAt                                                                                                                 >  1:332850/1‑51 (MQ=255)                                                                                                                                 cGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAAc                                           >  2:226299/1‑121 (MQ=255)                                                                                                                                 cGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAAc                                           <  1:226299/121‑1 (MQ=255)                                                                                                                                       ccTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCcgc        <  4:120651/150‑1 (MQ=255)                                                                                                                                           ttGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCaa    >  2:222418/1‑150 (MQ=255)                                                                                                                                             gAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACa  >  2:454604/1‑150 (MQ=255)                                                                                                                                              aaaaaaACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCa                    >  2:442706/1‑131 (MQ=255)                                                                                                                                              aaaaaaACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCa                    <  1:442706/131‑1 (MQ=255)                                                                                                                                               aaaaaaCGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACcatcat                                  <  2:155618/116‑1 (MQ=255)                                                                                                                                               aaaaaaCGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACcatcat                                  >  1:155618/1‑116 (MQ=255)                                                                                                                                                   |                                                                                                                                                AGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACA  >  NC_000913_3_bme_pgi/4400286‑4400575

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