BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	66,614	T→A	intergenic (‑64/+221)	araD ← / ← araA	L‑ribulose‑5‑phosphate 4‑epimerase/L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	66,614	0	T	A	100.0%	59.6 / NA	22	intergenic (‑64/+221)	araD/araA	L‑ribulose‑5‑phosphate 4‑epimerase/L‑arabinose isomerase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base A (13/9); total (13/9)

GCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTAT  >  NC_000913/66481‑66669

gCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATaa                                                    >  2:1701897/1‑139 (MQ=255)

              tgtTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTc                                      >  2:1650829/1‑139 (MQ=255)

                   tGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATc                                 <  1:1650829/139‑1 (MQ=255)

                     gCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGAcgcgt                                             >  1:1565163/1‑124 (MQ=255)

                     gCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGAcgcgt                                             <  2:1565163/125‑2 (MQ=255)

                        gcgcCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCg                            >  1:390643/1‑139 (MQ=255)

                                ttGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGttacatac                    >  2:722937/1‑133 (MQ=255)

                                  ggCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGttacataccg                  <  1:661737/139‑9 (MQ=18)

                                    cTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGttacataccgga                <  1:722937/139‑11 (MQ=17)

                                    cTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGttacataccgga                <  2:940685/139‑11 (MQ=17)

                                        tAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCACATCAGGCGTTACATACCGGATgcg            >  1:818117/1‑139 (MQ=17)

                                             ccTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATgcggctac       >  1:1155668/1‑134 (MQ=14)

                                                 gCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATTGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATgcggctacttag   >  2:1624855/1‑130 (MQ=11)

                                                  cGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAgg                              >  1:190848/1‑111 (MQ=255)

                                                  cGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAgg                              <  2:190848/111‑1 (MQ=255)

                                                  cGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATgcggctacttagc  >  2:231143/1‑129 (MQ=11)

                                                                       ttGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCa                                   >  1:605673/1‑85 (MQ=255)

                                                                       ttGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCa                                   <  2:605673/85‑1 (MQ=255)

                                                                       ttGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAgg                              <  2:1032125/90‑1 (MQ=255)

                                                                       ttGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAgg                              >  1:1032125/1‑90 (MQ=255)

                                                                         gACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCAt                                  <  1:985883/84‑1 (MQ=255)

                                                                         gACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCAt                                  >  2:985883/1‑84 (MQ=255)

GCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTAT  >  NC_000913/66481‑66669

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