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%	159.9 / NA	38	coding (226/684 nt)	yjiH	Gate family protein YjiH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (24/14);  total (24/14)

CCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGCCACCAGCTGTCTGTTCCGGAAAAAACACCGGAAAA  >  NC_000913/4560316‑4560587                                                                                                                                            |                                                                                                                                     cccAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCGGCTCCAGCTCAAACAGAAAGCCATCTTTAGc                                                                                                                                      >  2:252617/1‑140 (MQ=255) cccAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGc                                                                                                                                      >  2:438370/1‑140 (MQ=255) cccAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGc                                                                                                                                      >  1:356062/1‑140 (MQ=255)          ggTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCa                                                                                                                            >  1:102964/1‑140 (MQ=255)                     aGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACccgcgccgcg                                                                                                                 <  2:356062/140‑1 (MQ=255)                      ggTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTACAGCGCAACGAGTAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGa                                                                                                                >  2:288750/1‑140 (MQ=255)                            aTTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAg                                                                                                          >  1:151471/1‑140 (MQ=255)                               aaCTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAggtg                                                                                                       <  1:344435/140‑1 (MQ=255)                                 ctgctgCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAggtggt                                                                                                     >  2:383076/1‑140 (MQ=255)                                 ctgctgCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCAGCTTCAG‑CGCCGCGCCAACCGCGCCGCGAAATGAggtggt                                                                                                     >  2:177417/1‑140 (MQ=255)                                    ctgcGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACccgcgcc                                                                                                                    <  1:110368/122‑1 (MQ=255)                                    ctgcGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACccgcgcc                                                                                                                    >  2:110368/1‑122 (MQ=255)                                     tgcGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGGGCCGCGAAATGAGGTGGTGGCt                                                                                                 >  1:215042/1‑140 (MQ=255)                                          ggcgcgTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGtt                                                                                            <  1:487986/140‑1 (MQ=255)                                            cgcgTAATCCACCAAGAACATCCGTAATTGATATATTCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCg                                                                                          <  2:215042/140‑1 (MQ=255)                                                      accaAGACCATCGGTCATTGAAGTACTCCCCAGAGAGAGAATGACTTACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAg                                                                                <  1:315720/140‑1 (MQ=255)                                                        caAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCa                                                                                     >  1:409355/1‑133 (MQ=255)                                                        caAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCa                                                                                     <  2:409355/133‑1 (MQ=255)                                                           gACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCAACCGCGCCGCGAATTAAGGTGGTGGCTCCGTTAGCACCAGCCAGTTGcc                                                                           >  2:314145/1‑140 (MQ=255)                                                            aCCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGccc                                                                          >  1:435037/1‑140 (MQ=255)                                                              cATCGGTAATTGAAATAATCCCGAGCGAGAGAATGACAGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCaa                                                                        <  1:282678/140‑1 (MQ=255)                                                              cATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCaa                                                                        <  1:428709/140‑1 (MQ=255)                                                                  ggTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGa                                                                    >  2:523839/1‑140 (MQ=255)                                                                      aTTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCg                                                                >  1:588791/1‑140 (MQ=255)                                                                      aTTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCg                                                                >  2:339050/1‑140 (MQ=255)                                                                            ataatCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGa                                                          >  2:518554/1‑140 (MQ=255)                                                                              aatCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCACACCAGGCAGTTGCCCAAACGATCCGTTCAGAAc                                                        >  1:23788/1‑140 (MQ=255)                                                                                tCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGTCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACgg                                                      >  2:525944/1‑140 (MQ=255)                                                                                  cccAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTCGCTCCGTTCGCATCAGGCAGTTGCCCAAACTTTCCGTTCAGAACGGaa                                                    >  2:312589/1‑140 (MQ=255)                                                                                      gagagagaATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAAt                                                >  2:539518/1‑140 (MQ=255)                                                                                         agagaATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCaa                                             >  1:581797/1‑140 (MQ=255)                                                                                                gactgacGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAAcgcg                                      >  1:490936/1‑140 (MQ=255)                                                                                                      cGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGccacca                                <  2:102964/140‑1 (MQ=255)                                                                                                           ccAGCTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGCCACCAGctgt                           <  1:339050/140‑1 (MQ=255)                                                                                                               cTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGccac                                  <  2:600628/129‑1 (MQ=255)                                                                                                               cTCCAGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGccac                                  >  1:600628/1‑129 (MQ=255)                                                                                                                   aGCGCAAAGAGAAAGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGCCACCAGCTGTCTGTTCCg                   >  2:41512/1‑140 (MQ=255)                                                                                                                              aaaGCCATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGCCACCAGCTGTCTGTTCCGGAAAAAACAcc        <  2:92700/140‑1 (MQ=255)                                                                                                                                  ccATCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGCCACCAGCTGTCTGTTCCGGAAAAAACACCGGaa    >  2:484516/1‑140 (MQ=255)                                                                                                                                    aTCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGCCACCAGCTGTCTGTTCCGGAAAAAACACCGGaaaa  >  2:75142/1‑140 (MQ=255)                                                                                                                                    aTCTTTAG‑CCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGCCACCAGCTGTCTGTTCCGGAAAAAACACCGGaaaa  <  2:151471/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                     CCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCATCGGTAATTGAAATAATCCCCAGAGAGAGAATGACTGACGGTGCCAGCTCCAGCGCAAAGAGAAAGCCATCTTTAGCCCCCGCGCCACCCGCGCCGCGAAATGAGGTGGTGGCTCCGTTCGCACCAGGCAGTTGCCCAAACGATCCGTTCAGAACGGAAAAATCAAAAACGCGCCACCAGCTGTCTGTTCCGGAAAAAACACCGGAAAA  >  NC_000913/4560316‑4560587

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