BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,347,288	Δ1 bp	intergenic (+138/‑96)	nrdA → / → nrdB	ribonucleoside‑diphosphate reductase 1 subunit alpha/ribonucleoside‑diphosphate reductase 1 subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,347,288	0	A	.	100.0%	95.8 / NA	26	intergenic (+138/‑96)	nrdA/nrdB	ribonucleoside‑diphosphate reductase 1 subunit alpha/ribonucleoside‑diphosphate reductase 1 subunit beta
Reads supporting (aligned to +/- strand): ref base A (0/0); new base . (15/11); total (15/11)

GGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCATATACCACCTTTTCACAGACGAAAAATGAT  >  NC_000913/2347243‑2347419

agatGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAAc                                                 <  1:1155221/126‑1 (MQ=18)

agatGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAAc                                                 >  2:1155221/4‑129 (MQ=17)

agatGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTc                                      >  2:579694/4‑140 (MQ=18)

agatGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTc                                      >  2:409666/4‑140 (MQ=18)

 gatGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCa                                     <  1:297930/138‑1 (MQ=255)

  atGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGaca                                           <  1:404027/132‑1 (MQ=255)

  atGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGaca                                           >  2:404027/2‑133 (MQ=18)

    gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATgg                                  >  1:1099376/1‑140 (MQ=255)

    gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATgg                                  >  1:449483/1‑140 (MQ=255)

    gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATgg                                  <  2:364504/140‑1 (MQ=255)

      cGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCa                                <  1:9608/140‑1 (MQ=255)

            cggCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCATATacc                          >  1:963760/3‑140 (MQ=255)

             ggCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCATATacca                         >  2:379073/2‑140 (MQ=255)

              gCGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCATATaccac                        >  2:890305/1‑140 (MQ=255)

               cGTAAACGCCTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCATATaccacc                       >  2:994543/1‑140 (MQ=255)

                        cTTATCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCATATACCACCTTTTCACAg              >  1:72109/1‑140 (MQ=255)

                           aTCCGGCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCATATACCACCTTTTCACAGACg           >  2:756552/1‑140 (MQ=255)

                               ggCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCATATACCACCTTTTCACAGACGaaaa       >  1:932386/1‑140 (MQ=255)

                               ggCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCATATACCACCTTTTCACAGACGaaaa       >  1:140184/1‑140 (MQ=255)

                                gCCTACGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCATATACCACCTTTTCACAGACGaaaaa      >  1:644772/1‑140 (MQ=255)

                                    aCGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCATATACCACCTTTTCACAGACGAAAAATGAt  <  2:1153129/140‑1 (MQ=255)

                                    aCGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCATATACCACCTTTTCACAGACGAAAAATGAt  <  2:168314/140‑1 (MQ=255)

                                    aCGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCATATACCACCTTTTCACAGACGAAAAATGAt  <  1:994543/140‑1 (MQ=255)

                                    aCGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCATATACCACCTTTTCACAGACGAAAAATGAt  <  1:53618/140‑1 (MQ=255)

                                    aCGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCATATACCACCTTTTCACAGACGAAAAATGAt  <  2:963760/140‑1 (MQ=255)

                                    aCGGCTCGG‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCATATACCACCTTTTCACAGACGAAAAATGAt  <  1:1012221/140‑1 (MQ=255)

GGGTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGGCTCGGATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAGGACACACTCATGGCATATACCACCTTTTCACAGACGAAAAATGAT  >  NC_000913/2347243‑2347419

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