Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	885,221	T→C	K4E (AAA→GAA)	ybjH ←	uncharacterized protein YbjH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	885,221	0	T	C	100.0%	86.2 / NA	25	K4E (AAA→GAA)	ybjH	uncharacterized protein YbjH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/13);  total (12/13)

ATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGC  >  NC_000913/885147‑885298                                                                           |                                                                              aTCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTCCATAAt                                                                         >  2:167514/1‑81 (MQ=255)    gACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTCCATAATCa                                                                       <  1:167484/80‑1 (MQ=255)       ggTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTCCATAATCATTTc                                                                   >  1:167508/1‑81 (MQ=255)         tGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTCCATAATCATTTCGc                                                                 <  1:167485/81‑1 (MQ=255)                     ccATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATGc                                                     <  2:167486/81‑1 (MQ=255)                      cATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATg                                                      >  1:167507/1‑79 (MQ=255)                           ccACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCggg                                               >  1:167502/1‑81 (MQ=255)                             acacCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGtt                                             <  1:167487/81‑1 (MQ=255)                              cacCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGTTa                                            <  1:167488/81‑1 (MQ=255)                                      aaGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCcg                                     <  2:167489/80‑1 (MQ=255)                                        gCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCcgcg                                   >  1:167510/1‑80 (MQ=255)                                          ccATTAAAAGCGCGCCCAACAGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAAt                                <  1:167490/81‑1 (MQ=255)                                           cATTAAAAGCGCGCCCAACAGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAAtt                               <  2:167491/81‑1 (MQ=255)                                                  aGCGCGCCCAACAGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGc                        >  2:167516/1‑81 (MQ=255)                                                  aGCGCGCCCAACAGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGc                        >  1:167521/1‑81 (MQ=255)                                                   gcgcgcCCAACAGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGcc                       >  2:167512/1‑81 (MQ=255)                                                    cgcgcCCAACAGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCt                      >  1:167511/1‑81 (MQ=255)                                                        cccAACAGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCttcat                  >  1:167523/1‑81 (MQ=255)                                                         ccAACAGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCttcatt                 <  2:167492/81‑1 (MQ=255)                                                          cAACAGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCttcattc                <  2:167493/81‑1 (MQ=255)                                                            aCAGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAg              <  1:406481/81‑1 (MQ=255)                                                              aGTAGACAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAgg             >  1:167519/1‑80 (MQ=255)                                                                tAGACAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCAt          <  1:167494/81‑1 (MQ=255)                                                                   aCAATTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATtgt       >  2:167520/1‑81 (MQ=255)                                                                      aTTTTCCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGtt    <  1:167495/81‑1 (MQ=255)                                                                        tttcCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGc  >  1:167518/1‑81 (MQ=255)                                                                           |                                                                              ATCGACGGTGACACTTTGCGCCATCGCCACACCAGTAAAGCCCATTAAAAGCGCGCCCAACAGTAGACAATTTTTCATAATCATTTCGCCTGAAATAATGCAGCGGGTTAATTTCCGCGAATTATGCAGCCTTCATTCAGGCATTGTGTTGC  >  NC_000913/885147‑885298

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