Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	333,708	(A)6→7	coding (726/933 nt)	yahB ←	putative LysR‑type DNA‑binding transcriptional regulator YahB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	333,702	1	.	A	100.0%	106.6 / NA	30	L244F (TTG→TTT)	yahB	putative LysR‑type DNA‑binding transcriptional regulator YahB
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (17/13);  total (17/13)

GCGAATCCTGGCGCGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGC‑AAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCAGAATTGACTCCTGC  >  NC_000913/333623‑333777                                                                                 |                                                                            gcgAATCCTGGCGCGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGc                                                                              >  2:74114/1‑80 (MQ=255) gcgAATCCTGGCGCGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGc                                                                              >  1:74154/1‑80 (MQ=255) gcgAATCCTGGCGCGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGC‑a                                                                            <  2:74078/81‑1 (MQ=255) gcgAATCCTGGCGCGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGC‑a                                                                            <  2:74077/81‑1 (MQ=255)   gAATCCTGGCGCGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGC‑aaa                                                                          <  1:74079/81‑1 (MQ=255)   gAATCCTGGCGCGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGC‑aa                                                                           >  1:74112/1‑80 (MQ=255)   gAATCCTGGCGCGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGC‑aa                                                                           >  1:74118/1‑80 (MQ=255)       ccTGGCGCGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGC‑aaaaaaa                                                                      >  2:74121/1‑80 (MQ=255)           gcgcGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGC‑aaaaaaacc                                                                    <  1:74080/79‑4 (MQ=255)              cgGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAAACCAATTc                                                                <  2:74081/81‑1 (MQ=255)               ggATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAAACCAATTcc                                                               >  2:74133/1‑81 (MQ=255)               ggATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAAACCAATTcc                                                               >  2:74115/1‑81 (MQ=255)                gATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCAGCAAAAAAACCAATTCCt                                                              <  2:74082/81‑1 (MQ=255)                 aTTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAAACCAATTCCtt                                                             >  2:74127/1‑81 (MQ=255)                           tGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTaaaa                                                    <  2:74083/80‑1 (MQ=255)                           tGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAAt                                                   >  1:74117/1‑81 (MQ=255)                            gTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATc                                                  >  2:74120/1‑81 (MQ=255)                                    aCCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTg                                          <  1:74084/81‑1 (MQ=255)                                    aCCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAAACCAATCCCTTCACCTAAAATCTGGCACTg                                          <  2:74085/81‑1 (MQ=255)                                      cagcagGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGa                                         >  2:74131/1‑80 (MQ=255)                                       agcagGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACa                                       >  2:74135/1‑81 (MQ=255)                                              aTTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATTTGGCACTGACATTAGgtg                                >  1:74128/1‑81 (MQ=255)                                              aTTGCGTCATCGCCTCACGGACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGgtg                                >  1:74134/1‑81 (MQ=255)                                                   gTCATCGCCTCACGGACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGaa                           <  1:74086/81‑1 (MQ=255)                                                      aTCGCCTCACGGACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATc                        <  1:74087/81‑1 (MQ=255)                                                         gCCTCACGGACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTgg                     >  2:74146/1‑81 (MQ=255)                                                             cACGGACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCAc                  >  2:74141/1‑80 (MQ=255)                                                              aCGGACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCAcc                 <  1:74088/80‑1 (MQ=255)                                                              aCGGACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCAcc                 <  1:74089/80‑1 (MQ=255)                                                              aCGGACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCa                <  2:74090/81‑1 (MQ=255)                                                               cGGACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCAg               <  2:74091/81‑1 (MQ=255)                                                                ggACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCAg               >  2:74155/1‑80 (MQ=255)                                                                 gACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCAGaa             >  2:74148/1‑81 (MQ=255)                                                                 gACCATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCAGaa             >  1:74139/1‑81 (MQ=255)                                                                    cATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCAGAAtt           >  2:74124/1‑80 (MQ=255)                                                                    cATGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCAGAATTg          <  1:74092/81‑1 (MQ=255)                                                                      tGTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCAGAATTGAc        <  2:74093/81‑1 (MQ=255)                                                                       gTAATCCGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCAGAATTGACt       >  2:74142/1‑81 (MQ=255)                                                                            ccGGCAAAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCAGAATTGACTCCTGc  >  1:74129/1‑81 (MQ=255)                                                                                 |                                                                            GCGAATCCTGGCGCGGATTATGGATTTGTCGGGTAACCAGCAGGGATTGCGTCATCGCCTCACGGACCATGTAATCCGGC‑AAAAAACCAATTCCTTCACCTAAAATCTGGCACTGACATTTGGTGTTGAAATCTGGCACCAGAATTGACTCCTGC  >  NC_000913/333623‑333777

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