BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,762,741	A→T	K662I (AAA→ATA)	spoT →	bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,762,741	0	A	T	100.0%	93.1 / NA	27	K662I (AAA→ATA)	spoT	bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base T (14/13); total (14/13)

TGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATG  >  REL606/3762670‑3762811

tGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGagatagat                                                                    >  1:3443/1‑76 (MQ=255)

 gcgcTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATg                                                                   >  2:309068/1‑76 (MQ=255)

   gcTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGt                                                                 <  2:551879/76‑1 (MQ=255)

    cTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTc                                                                <  2:274812/76‑1 (MQ=255)

           aCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa                                                         <  1:501545/76‑1 (MQ=255)

             cTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa                                                       <  1:901106/76‑1 (MQ=255)

               gACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCNACAgc                                                     >  2:775864/1‑76 (MQ=255)

                aCGGCGGCAATTAACACCACGACTTCGACTATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAgcg                                                    <  2:437450/76‑1 (MQ=255)

                    cggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCtt                                                >  2:445316/1‑76 (MQ=255)

                    cggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCtt                                                <  2:437055/76‑1 (MQ=255)

                         aTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTc                                           >  2:278346/1‑76 (MQ=255)

                          taaaCACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTTGCGTCTACAGCGCCTTTATTCg                                          >  1:694087/3‑76 (MQ=255)

                                   ctaCTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGc                                 >  2:535453/3‑76 (MQ=255)

                                    gACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCt                                >  2:722314/1‑76 (MQ=255)

                                       ttCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGt                             <  2:15157/76‑1 (MQ=255)

                                         cGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGa                           <  1:868864/76‑1 (MQ=255)

                                          gAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGAc                          <  1:232941/76‑1 (MQ=255)

                                              atTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtg                      >  2:248298/1‑76 (MQ=255)

                                              atTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtg                      <  1:914720/76‑1 (MQ=255)

                                                     gTTTGAATACGGCATAGATAGATTNTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCNCGTGACCGTGTGCATCt               <  2:797337/76‑1 (MQ=255)

                                                        tGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGc            >  2:384745/1‑76 (MQ=255)

                                                         gAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGCCCGTGTGCATCTGGCg           >  2:612242/1‑76 (MQ=255)

                                                         gAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCg           >  1:374394/1‑76 (MQ=255)

                                                          aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGa          >  2:369320/1‑76 (MQ=255)

                                                           aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGa                           <  2:226882/58‑1 (MQ=255)

                                                           aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGa                           >  1:226882/1‑58 (MQ=255)

                                                                 aaGAGATAGATGGTCGCGTCTACAGCCCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGACTATCAt   >  1:859490/1‑76 (MQ=255)

                                                                  agagATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATg  <  1:281988/76‑1 (MQ=255)

TGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATG  >  REL606/3762670‑3762811

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