Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,252,023	A→T	D246E (GAT→GAA)	nusA ←	transcription elongation factor NusA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,252,023	0	A	T	100.0%	113.6 / NA	31	D246E (GAT→GAA)	nusA	transcription elongation factor NusA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (13/18);  total (13/18)

CCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTA  >  REL606/3251958‑3252092                                                                  |                                                                      ccACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGTTCGATACGtt                                                             <  1:797795/76‑1 (MQ=255)  cACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGTTCGATACGttt                                                            <  1:730248/76‑1 (MQ=255)   ccccAGTTCAGTAGACCCCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTa                                                           <  2:1186119/75‑1 (MQ=255)   ccccAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTa                                                           <  1:546709/75‑1 (MQ=255)       aGTTCAGTAGAAACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTATCGt                                                       >  2:299488/1‑76 (MQ=255)         ttCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTATCGTTg                                                     >  1:763168/1‑76 (MQ=255)          tCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTATCGTTgg                                                    >  2:932688/1‑76 (MQ=255)            aGTANACACCGCCTGAACACGCGCGCCATGCATACCTACNCAAGCACCTACCGGTTCGATACGTTTATCGTTGGtt                                                  <  1:141747/76‑1 (MQ=255)            aGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTATCGTTGGtt                                                  <  1:706250/76‑1 (MQ=255)             gTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTATCGTTGGttt                                                 >  1:975597/1‑76 (MQ=255)                 acacCGCCTGAACTCGCGCGCCATGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTATCGTTGGTTTTCAc                                             <  1:1038127/76‑1 (MQ=255)                 acacCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTATCGTTGGTTTTCAc                                             <  1:399537/76‑1 (MQ=255)                 acacCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTATCGTTGGTTTTCAc                                             <  1:851689/76‑1 (MQ=255)                   acCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTATCGTTGGTTTTCACcg                                           <  1:815429/76‑1 (MQ=255)                      nccTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTATCGTNGGTTTTCACCGCGa                                        >  2:174611/2‑76 (MQ=255)                        cTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTATCGTTGGTTTTCACCGCGAtt                                      >  1:441371/1‑76 (MQ=255)                          gAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTATCGTTGGTTTTCACCGCGGtttt                                    <  2:1197195/76‑1 (MQ=255)                          gAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTATCGTTGGTTTTCACCGCGAtttt                                    <  2:945403/76‑1 (MQ=255)                          gAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTATCGTTGGTTTTCACCGCGAtttt                                    <  2:385902/76‑1 (MQ=255)                                gcgcgcCACGCATATCTACGCAAGCACCTACCGGTTCGATACGTTTATCGTTGGTTTTGACCGCGATTTTCGCACg                              <  2:176162/76‑1 (MQ=255)                                   cgcCACGCATACCTACGCAAGCACCTACCGTTTCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACgaga                           <  2:108493/76‑1 (MQ=255)                                    gcCACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTATCGTTGGTTTTTACCGCGATTTTTGCCCgagga                          >  2:912692/1‑74 (MQ=255)                                    gcCACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGaa                          >  1:883548/1‑76 (MQ=255)                                     ccACGCATACCTACGCAAGCACCTACCGGTTCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAAc                         <  2:716497/76‑1 (MQ=255)                                               cTACGCAAGCACCTACCGGTTCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATcgcg               >  2:749343/1‑76 (MQ=255)                                                 aCGCAAGCACCTACCGGTTCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAg             >  2:1201212/1‑76 (MQ=255)                                                    cAAGCACCTACCGGTTCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCg          >  1:832576/1‑76 (MQ=255)                                                     aaGCACCTACCGGTTCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCgc         >  1:801356/1‑76 (MQ=255)                                                     aaGCACCTACCGGTTCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCgc         >  1:631413/1‑76 (MQ=255)                                                       gCACCTCCGGGTTCGATACGTTTATCGTGGGTTTTCACCGGGATTTTCGCACGAGACCCCGGATCGCGAGCCgctg       <  2:631005/76‑1 (MQ=255)                                                           ctcccGGTTCGATACGTTTATCGTTGGTTTTCACCGGGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCttt   <  2:621044/73‑1 (MQ=255)                                                            tACCGGTTCGATACGTTCATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTa  <  1:134156/76‑1 (MQ=255)                                                                  |                                                                      CCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTA  >  REL606/3251958‑3252092

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