Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,239,166	C→T	N511N (AAC→AAT)	yjbH →	YjbH family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,239,166	0	C	T	100.0%	108.7 / NA	31	N511N (AAC→AAT)	yjbH	YjbH family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/16);  total (15/16)

TCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGG  >  NC_000913/4239089‑4239243                                                                              |                                                                              tCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAattac                                                                            >  1:866570/1‑81 (MQ=255)  cGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAattact                                                                           >  1:866586/1‑81 (MQ=255)   gCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAATTACtt                                                                          >  2:866589/1‑81 (MQ=255)     aCTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAATTACTTc                                                                         >  1:866598/1‑80 (MQ=255)        tACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAATTACTTCCAg                                                                      >  2:866579/1‑80 (MQ=255)            gcgcgTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAATTACTTCCAGCATCt                                                                 <  2:866532/81‑1 (MQ=255)              gcgtgcgtACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAATTACTTCCAGCATCTgg                                                               <  2:866533/81‑1 (MQ=255)               cgtgcgtACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAATTACTTCCAGCATCTggg                                                              <  1:866534/81‑1 (MQ=255)                gtgcgtACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAATTACTTCCAGCATCTGGGc                                                             <  2:866535/81‑1 (MQ=255)                    gtACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAATTACTTCCAGCATCTGGGCAACg                                                         >  2:866582/1‑81 (MQ=255)                      aCCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAATTACTTCCAGCATCTGGGCAACGGc                                                       >  2:866592/1‑81 (MQ=255)                      aCCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAATTACTTCCAGCATCTGGGCAACGGc                                                       >  2:866591/1‑81 (MQ=255)                        ccATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAATTACTTCCAGCATCTGGGCAACGGCtt                                                     >  2:866583/1‑81 (MQ=255)                           tgtgCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTa                                                  <  2:866536/81‑1 (MQ=255)                             tgCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTAc                                                 <  1:866537/80‑1 (MQ=255)                               cgcgAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTACGGt                                              <  2:866538/81‑1 (MQ=255)                               cgcgAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTACGGt                                              >  1:866608/1‑81 (MQ=255)                                      aTGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTc                                        >  1:866571/1‑80 (MQ=255)                                        gtgCAGAACGATGTCTATGTGAATAACCTGCAAGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTAc                                     >  2:866585/1‑81 (MQ=255)                                                cGATGTCTATGTGAATAACCTGCAAGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTa                             >  2:866563/1‑81 (MQ=255)                                                     tCTATGTGAATAACCTGCAAGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTAtctc                         >  1:866604/1‑80 (MQ=255)                                                         tgtgAATAACCTGCAAGCCAATTACCTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGa                       <  2:866539/78‑1 (MQ=255)                                                           tgAATAACCTGCAAGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAAcc                   <  2:866542/80‑1 (MQ=255)                                                           tgAATAACCTGCAAGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAAcc                   <  2:866540/80‑1 (MQ=255)                                                           tgAATAACCTGCAAGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCa                  >  1:866590/1‑81 (MQ=255)                                                           tgAATAACCTGCAAGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCa                  <  1:866541/81‑1 (MQ=255)                                                              aTAACCTGCAAGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCATGt               <  2:866543/81‑1 (MQ=255)                                                                  ccTGCAAGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCATGTTTgg           <  2:866544/81‑1 (MQ=255)                                                                     gCAAGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCATGTTTGGCgg        >  1:866600/1‑81 (MQ=255)                                                                      cAAGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCATGTTTGGCGGt       <  1:866545/81‑1 (MQ=255)                                                                       aaGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCATGTTTGGCGGTg      <  2:866546/81‑1 (MQ=255)                                                                       aaGCCAATTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCATGTTTGGCGGTg      >  2:866610/1‑81 (MQ=255)                                                                           cAATTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCATGTTTGGCGGTGCggg  <  2:866547/81‑1 (MQ=255)                                                                              |                                                                              TCGCACTTACCGCGCGTGCGTACCCATGTGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGG  >  NC_000913/4239089‑4239243

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