Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,237,643	(A)7→8	coding (10/2097 nt)	yjbH →	YjbH family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,237,636	1	.	A	100.0%	85.5 / NA	24	M1M (ATG→ATA) †	yjbH	YjbH family protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (13/11);  total (13/11)

CCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATG‑AAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCC  >  NC_000913/4237561‑4237710                                                                             |                                                                           ccGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATG‑aaaa                                                                        >  2:846499/1‑80 (MQ=255)          gTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAAGACATc                                                               <  2:846464/81‑1 (MQ=255)           tGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAAGACATCt                                                              <  2:846465/81‑1 (MQ=255)                cGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAAGACATCTgctta                                                         <  1:846466/81‑1 (MQ=255)                 gCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAAGACATCTgcttag                                                        <  1:846467/81‑1 (MQ=255)                 gCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAAGACATCTgcttag                                                        <  2:846468/81‑1 (MQ=255)                  ccTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAAGACATCTgcttagc                                                       <  2:846469/81‑1 (MQ=255)                   cTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAAGACATCTgcttagct                                                      >  2:846505/1‑81 (MQ=255)                       tGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAAGACATCTGCTTAGCTTACt                                                  <  2:846470/81‑1 (MQ=255)                               aCGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAAGACATCTGCTTAGCTTACTGGCGCTgg                                          >  1:846520/1‑81 (MQ=255)                                cGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAAGACATCTGCTTAGCTTACTGGCGCTggg                                         <  2:846471/81‑1 (MQ=255)                                   cGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCAt                                      >  1:846521/1‑81 (MQ=255)                                    gACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAAGACATCTGCTTAGTTTACTGGCGCTGGGCAt                                      >  1:846532/1‑80 (MQ=255)                                     aCATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCAt                                      >  1:846502/1‑79 (MQ=255)                                         tcttcAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGca                                 >  1:846517/1‑80 (MQ=255)                                           ttcAGACTCTGACGCAGCGGATACCTCAATAATGAAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAg                              >  1:846510/1‑81 (MQ=255)                                                   ctGACGCAGCGGATACCTCAATAATGAAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTAc                      >  2:846511/1‑81 (MQ=255)                                                        gCAGCGGATACCTCAATAATGAAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCg                  >  1:846533/1‑80 (MQ=255)                                                        gCAGCGGATACCTCAATAATGAAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGa                 <  2:846472/81‑1 (MQ=255)                                                        gCAGCGGATACCTCAATAATGAAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGa                 <  2:846473/81‑1 (MQ=255)                                                           gCGGATACCTCAATAATGAAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGaaa               >  2:846529/1‑80 (MQ=255)                                                            cGGATACCTCAATAATGAAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACa             >  2:846513/1‑81 (MQ=255)                                                             ggATACCTCAATAATGAAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACat            <  2:846474/81‑1 (MQ=255)                                                             ggATACCTCAATAATGAAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACa             >  1:846508/1‑80 (MQ=255)                                                                       ataatGAAAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGcc  >  2:846509/1‑81 (MQ=255)                                                                             |                                                                           CCGTCTGGAGTGAGACGCCTGATGCCCTTAACGCCGACATTCTTCAGACTCTGACGCAGCGGATACCTCAATAATG‑AAAAAAAGACATCTGCTTAGCTTACTGGCGCTGGGCATTAGCACAGCTTGCTACGGCGAAACATATCCTGCGCC  >  NC_000913/4237561‑4237710

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