BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	765,208	Δ5,405 bp		[mngR]–mngB	[mngR], mngA, mngB

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	765208	770612	5405	38 [0]	[0] 38	[mngR]–mngB	[mngR],mngA,mngB

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 765207	0 (0.000)	38 (1.280)	25/96	0.0	100%	coding (669/723 nt)	mngR	DNA‑binding transcriptional repressor MngR
*	?	NC_000913	770613 =	0 (0.000)	38 (1.280)	25/96	0.0	100%	intergenic (+2/‑845)	mngB/cydA	alpha‑mannosidase/cytochrome bd‑I ubiquinol oxidase subunit I

CGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCC  >  NC_000913/765159‑765254

cgATGATTTTCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAta                                               <  1:311374‑M1/50‑3 (MQ=255)

 gatgatTTTCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtag                                              >  1:167373‑M1/1‑47 (MQ=255)

  atgatTTTCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagc                                             <  1:2140299‑M1/50‑5 (MQ=255)

     atTTTCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaat                                          >  1:1844422‑M1/1‑43 (MQ=255)

      ttttCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatc                                         <  1:1263530‑M1/50‑9 (MQ=255)

      ttttCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatc                                         <  1:410685‑M1/50‑9 (MQ=255)

      ttttCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaat                                          >  1:1771179‑M1/1‑42 (MQ=255)

       tttCGCTGCGCTATCCGGGTAAA‑TTATAGTCATCGGTATTAtagcaatc                                         >  1:2340897‑M1/1‑41 (MQ=255)

       tttCGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatc                                         >  1:578162‑M1/1‑41 (MQ=255)

          cGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaac                                     >  1:963074‑M1/1‑38 (MQ=255)

          cGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaac                                     >  1:2659194‑M1/1‑38 (MQ=255)

          cGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaac                                     >  1:1706871‑M1/1‑38 (MQ=255)

          cGCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaa                                      <  1:678574‑M1/49‑12 (MQ=255)

           gCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaacc                                    <  1:1651447‑M1/50‑14 (MQ=255)

           gCTGCGCTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaacc                                    <  1:2512072‑M1/50‑14 (MQ=255)

              gcgcTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccg                                   >  1:1603340‑M1/1‑34 (MQ=255)

              gcgcTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccg                                   >  1:1415707‑M1/1‑34 (MQ=255)

               cgcTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaag                                <  1:1738243‑M1/50‑18 (MQ=255)

               cgcTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaag                                >  1:349574‑M1/1‑33 (MQ=255)

               cgcTATCAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaag                                >  1:1394880‑M1/1‑33 (MQ=255)

                     cAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagccaca                           >  1:248087‑M1/1‑27 (MQ=255)

                     cAGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagccaca                           >  1:2629528‑M1/1‑27 (MQ=255)

                      aGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagccacata                         >  1:1987153‑M1/1‑26 (MQ=255)

                      aGGGTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagccacata                         >  1:1045023‑M1/1‑26 (MQ=255)

                       gggTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagccacatat                        >  1:2382705‑M1/1‑25 (MQ=255)

                       gggTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagccacatat                        >  1:232973‑M1/1‑25 (MQ=255)

                        ggTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagccacatatg                       >  1:599784‑M1/1‑24 (MQ=255)

                         gTAAA‑TTATAGTCATCGGTATTAtagcaatcaaaccgaagccacatatgc                      >  1:1242367‑M1/1‑23 (MQ=255)

                             aattatAGTCATCGGTATTAtagcaatcaaaccgaagccacatatgcgc                    <  1:240991‑M1/47‑30 (MQ=255)

                              attaTAGTCATCGGTATTAtagcaatcaaaccgaagccacatatgcgcgg                  >  1:2429232‑M1/2‑19 (MQ=255)

                                tataGTCATCGGTATTAtagcaatcaaaccgaagccacatatgcgcggc                 >  1:1661357‑M1/1‑17 (MQ=255)

                                 ataGTCATCGGTATTAtagcaatcaaaccgaagccacatatgcgcggcc                <  1:1327627‑M1/49‑34 (MQ=255)

                                  taGTCATCGGTATTAtagcaatcaaaccgaagccacatatgcgcggc                 >  1:799354‑M1/1‑15 (MQ=255)

                                   aGTCATCGGTATTAtagcaatcaaaccgaagccacatatgcgcggc                 >  1:1470190‑M1/1‑14 (MQ=255)

                                   aGTCATCGGTATTAtagcaatcaaaccgaagccacatatgcgcggc                 >  1:140194‑M1/1‑14 (MQ=255)

                                      cATCGGTATTAtagcaatcaaaccgaagccacatatgcgcggccagattg          >  1:179561‑M1/1‑11 (MQ=255)

                                           gTATTAtagcaatcaaaccgaagccacatatgcgcggccagattgttgac     <  1:1948985‑M1/50‑45 (MQ=255)

                                              ttatagcaatcaaaccgaagccacatatgcgcggccagattgttgacaaa  <  1:623503‑M1/50‑48 (MQ=255)

CGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCC  >  NC_000913/765159‑765254

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