BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,451,734	T→C	intergenic (+172/‑496)	unknown → / → unknown	FIG00698722: hypothetical protein/Bacterial non‑heme ferritin (EC 1.16.3.2)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,451,734	0	T	C	100.0%	44.0 / NA	18	intergenic (+172/‑496)	unknown/unknown	FIG00698722: hypothetical protein/Bacterial non‑heme ferritin (EC 1.16.3.2)
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (8/10); total (8/10)

CAGCTATCCAATTTTATTAAGACGGATAAAGGCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACACTTCTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTTCATCACGTTTAACTAGTAATAGTGACCTGTGAATAACTCTGTGTTTTCATTAAATCGTAGTGTTCATTTACTTTGTCACTAATTTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGGTTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTC  >  1/1451524‑1451932

aaatttcccgcccctgcaaatacccctaaagccacgtaataactgctgactagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATAtatta                                                                                                                                                                <  1:358946‑M2/164‑1 (MQ=255)

   tttcccgcccctgcaaatacccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGt                                                                                                                                                             <  1:641107‑M2/167‑1 (MQ=255)

                   atacccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATgg                                                                                                                                             <  1:334688‑M2/183‑1 (MQ=255)

                      cccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATAtatt                                                                                                                                                                 <  1:332040‑M2/163‑1 (MQ=255)

                      cccctaaagccacgtaataactgctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTaa                                                                                                                                          >  1:814500‑M2/66‑251 (MQ=255)

                                            gctgattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATAtatt                                                                                                                                                                 >  1:591573‑M2/44‑206 (MQ=255)

                                               gattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACatatat                                                                                                                                                                  <  1:873395‑M2/162‑1 (MQ=255)

                                               gattagataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATAtatta                                                                                                                                                                <  1:882711‑M2/164‑1 (MQ=255)

                                                   agataatgcgatacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACatatat                                                                                                                                                                  >  1:413914‑M2/37‑198 (MQ=255)

                                                             atacgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTaaaaa                                                                                                    <  1:732398‑M2/224‑1 (MQ=255)

                                                                cgtgcttgaatttatgcttacatTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGa                                                                                                 >  1:684168‑M2/24‑250 (MQ=255)

                                                                                                               atatTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGt                                                                                                                                                             >  1:796235/1‑143 (MQ=255)

                                                                                                                            ttACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATAtatt                                                                                                                                                                 <  1:659191/126‑1 (MQ=255)

                                                                                                                            ttACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTaaaaa                                                                                                    >  1:661888/1‑187 (MQ=255)

                                                                                                                                   gCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACat                                                                                                                                                                      >  1:676356/1‑114 (MQ=255)

                                                                                                                                           tcttcAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTag                     <  1:352855/251‑1 (MQ=255)

                                                                                                                                            cttcAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTaga                    >  1:313791/1‑251 (MQ=255)

                                                                                                                                                              ttAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTtc  <  1:807161/251‑1 (MQ=255)

CAGCTATCCAATTTTATTAAGACGGATAAAGGCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACACTTCTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTTCATCACGTTTAACTAGTAATAGTGACCTGTGAATAACTCTGTGTTTTCATTAAATCGTAGTGTTCATTTACTTTGTCACTAATTTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGGTTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTC  >  1/1451524‑1451932

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