Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,560,459	(C)5→4	coding (222/684 nt)	yjiH ←	Gate family protein YjiH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,560,455	0	C	.	100.0%	113.8 / NA	28	coding (226/684 nt)	yjiH	Gate family protein YjiH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (14/14);  total (14/14)

CCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGCCACCAGCTGTCTGTTCCGGAAAAAACACCGGAAAAAA  >  NC_000913/4560316‑4560589                                                                                                                                            |                                                                                                                                       cccAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGGGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGc                                                                                                                                        >  1:265598/1‑140 (MQ=255) cccAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGc                                                                                                                                        >  2:259011/1‑140 (MQ=255)     acaacGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGccccg                                                                                                                                    <  1:496074/140‑2 (MQ=255)               tAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACccgc                                                                                                                         <  2:624197/140‑1 (MQ=255)                aGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACccgcg                                                                                                                        <  2:265598/140‑1 (MQ=255)                       gtgtCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGaa                                                                                                                 <  2:150007/140‑1 (MQ=255)                             ttAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAgg                                                                                                           <  2:288348/140‑1 (MQ=255)                               aaCTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAggtg                                                                                                         >  2:505277/1‑140 (MQ=255)                                 ctgctgCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAggtggt                                                                                                       >  1:156644/1‑140 (MQ=255)                                      gcgcggcgcgTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACccgcgccgc                                                                                                                    <  2:69586/122‑1 (MQ=255)                                      gcgcggcgcgTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACccgcgccgc                                                                                                                    >  1:69586/1‑122 (MQ=255)                                         cggcgcgTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGt                                                                                               <  1:259011/140‑1 (MQ=255)                                                      accaAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAg                                                                                  >  1:160295/1‑140 (MQ=255)                                                      accaAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAg                                                                                  >  1:334133/1‑140 (MQ=255)                                                        caAGACCATCGGAAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGtt                                                                                >  1:178589/1‑140 (MQ=255)                                                           gACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGcc                                                                             >  2:331121/1‑140 (MQ=255)                                                                              aatCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAAc                                                          >  1:374113/1‑140 (MQ=255)                                                                                             aaTGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAAc                                           >  1:401453/1‑140 (MQ=255)                                                                                               tgactgacGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAAcgc                                         <  1:331121/140‑1 (MQ=255)                                                                                                 actgacGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACgcgc                                       >  1:446625/1‑140 (MQ=255)                                                                                                  ctgacGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGcc                                      <  2:334133/140‑1 (MQ=255)                                                                                                   tgacGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGcca                                     >  2:622896/1‑140 (MQ=255)                                                                                                   tgacGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGcca                                     >  1:434666/1‑140 (MQ=255)                                                                                                    gacGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGccac                                    <  2:160295/140‑1 (MQ=255)                                                                                                    gacGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGccac                                    <  1:622896/140‑1 (MQ=255)                                                                                                      cGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGccacca                                  <  2:401453/140‑1 (MQ=255)                                                                                                           cccgCTCCAGTGCATAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGCCACCAGctgt                             <  2:446625/137‑1 (MQ=255)                                                                                                              gCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGCCACCAGctgtctg                          <  2:156644/140‑1 (MQ=255)                                                                                                                         aaGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGCCACCAGCTGTCTGTTCCGGaaaaa               >  2:566894/1‑140 (MQ=255)                                                                                                                                    aTCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGCCACCAGCTGTCTGTTCCGGAAAAAACACCGGaaaa    <  2:178589/140‑1 (MQ=255)                                                                                                                                      cTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGCCACCAGCTGTCTGTTCCGGAAAAAACACCGGaaaaaa  >  1:371518/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                       CCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGCCACCAGCTGTCTGTTCCGGAAAAAACACCGGAAAAAA  >  NC_000913/4560316‑4560589

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