BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	49,765	C→T	intergenic (+134/‑58)	kefC → / → folA	glutathione‑regulated potassium‑efflux system protein KefC/type 3 dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	49,765	0	C	T	100.0%	100.9 / NA	26	intergenic (+134/‑58)	kefC/folA	glutathione‑regulated potassium‑efflux system protein KefC/type 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (11/15); total (11/15)

CATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCA  >  NZ_CP009273/49681‑49829

cATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTtgacga                                                             <  2:131461/90‑1 (MQ=255)

      aaTGCGGCGAGTCCAGGGAGAGAGCGTGGGCTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGttt                                                       <  2:164742/90‑1 (MQ=255)

      aaTGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGttt                                                       <  1:647884/90‑1 (MQ=255)

        tGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTAc                                                     <  2:81127/90‑1 (MQ=255)

           ggcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCt                                                  <  1:510613/90‑1 (MQ=255)

            gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCtt                                                 >  1:266271/1‑90 (MQ=255)

            gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCtt                                                 >  1:773882/1‑90 (MQ=255)

            gcgAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCtt                                                 <  2:773882/90‑1 (MQ=255)

              gAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTa                                               <  1:508360/90‑1 (MQ=255)

                         gagagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGGTGACGACGGTTTACGCTTTACGTATAGGGGc                                    >  1:329058/1‑90 (MQ=255)

                          agagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCg                                   >  2:160948/1‑90 (MQ=255)

                           gagCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGa                                  >  2:887911/1‑90 (MQ=255)

                               gTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGAc                                 <  2:951035/87‑1 (MQ=255)

                               gTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGAc                                 >  1:951035/1‑87 (MQ=255)

                                      cGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAAtttttttt                       >  2:706497/1‑90 (MQ=255)

                                      cGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAAtttttttt                       >  1:523865/1‑90 (MQ=255)

                                           gcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACg                                                    >  2:88597/1‑56 (MQ=255)

                                           gcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACg                                                    <  1:88597/56‑1 (MQ=255)

                                           gcagAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCgg                  <  1:604836/90‑1 (MQ=255)

                                              gAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaa               >  1:174821/1‑90 (MQ=255)

                                              gAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaa               <  1:770829/90‑1 (MQ=255)

                                                          ttCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATc   <  2:144023/90‑1 (MQ=255)

                                                          ttCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATc   <  2:675910/90‑1 (MQ=255)

                                                          ttCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATc   <  1:792660/90‑1 (MQ=255)

                                                          ttCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATc   >  1:211357/1‑90 (MQ=255)

                                                           tCCTCAACATCATCCTCGCACCAGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCa  <  1:665408/90‑1 (MQ=255)

                                                                                 aGTTGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGaa               <  2:60298/55‑1 (MQ=255)

CATTGTAATGCGGCGAGTCCAGGGAGAGAGCGTGGACTCGCCAGCAGAATATAAAATTTTCCTCAACATCATCCTCGCACCAGTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCA  >  NZ_CP009273/49681‑49829

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