BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	459,278	Δ69 bp	coding (391‑459/2355 nt)	lon →	Lon protease

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	459278	459346	69	52 [1]	[0] 49	lon	Lon protease

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 459277	1 (0.020)	41 (0.820)	25/110	0.2	99.0%	coding (390/2355 nt)	lon	Lon protease
*	?	NC_000913	459347 =	0 (0.000)	41 (0.820)	25/110	0.2	99.0%	coding (460/2355 nt)	lon	Lon protease

ACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGA  >  NC_000913/459213‑459339

acTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTG                                                                >  1:2622498‑M1/1‑65 (MQ=255)

acTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTG                                                                >  1:3857055‑M1/1‑65 (MQ=255)

     tCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGacgtc                                                           >  1:638535‑M1/1‑60 (MQ=255)

               cGGATTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagc                                                 <  1:2035454‑M1/65‑16 (MQ=255)

               cGGAGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagc                                                 <  1:3631629‑M1/65‑16 (MQ=255)

                  aGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatc                                              >  1:499844‑M1/1‑47 (MQ=255)

                  aGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatc                                              <  1:3392343‑M1/65‑19 (MQ=255)

                  aGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGc                                              >  1:43573/1‑65 (MQ=255)

                   ggatcTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcg                                             <  1:1244098‑M1/63‑20 (MQ=255)

                    tATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcga                                            >  1:2422312‑M1/1‑45 (MQ=255)

                    tATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcga                                            >  1:152000‑M1/1‑45 (MQ=255)

                    tATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcga                                            >  1:695591‑M1/1‑45 (MQ=255)

                       cTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacga                                         >  1:4062261‑M1/1‑42 (MQ=255)

                        tGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgat                                        >  1:3358409‑M1/1‑41 (MQ=255)

                        tGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgat                                        >  1:3118209‑M1/1‑41 (MQ=255)

                          gAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatcc                                      >  1:2419179‑M1/1‑39 (MQ=255)

                           aGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccg                                     >  1:2487223‑M1/1‑38 (MQ=255)

                           aGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccg                                     >  1:2665646‑M1/1‑38 (MQ=255)

                                ccGACCATTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcg                                <  1:2596035‑M1/65‑33 (MQ=255)

                                     cATTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctgg                           <  1:2449576‑M1/65‑38 (MQ=255)

                                     cATTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctgg                           >  1:3647294‑M1/1‑28 (MQ=255)

                                      aTTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggc                          >  1:3924347‑M1/1‑27 (MQ=255)

                                      aTTGATGAGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggc                          <  1:3057613‑M1/65‑39 (MQ=255)

                                        tgatgaGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcgg                        <  1:258068‑M1/65‑41 (MQ=255)

                                        tgatgaGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcgg                        <  1:2819017‑M1/65‑41 (MQ=255)

                                        tgatgaGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcgg                        <  1:1947633‑M1/65‑41 (MQ=255)

                                        ggatgagCGAGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcgg                        <  1:2423222‑M1/64‑41 (MQ=255)

                                           tgaGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggata                     <  1:1160044‑M1/65‑44 (MQ=255)

                                           tgaGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggata                     <  1:3428354‑M1/65‑44 (MQ=255)

                                             aGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggatact                   >  1:4100706‑M1/1‑20 (MQ=255)

                                             aGCGGGAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggatacc                   <  1:1651881‑M1/65‑46 (MQ=255)

                                             aGCGGGAACAGGAAGTGCTGacgtcgctgaatagcaccgacgatccggcgcgtctggcggatacc                   >  1:1245175‑M1/1‑20 (MQ=255)

                                             aGCGGGAACAGAAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggatacc                   <  1:3451249‑M1/65‑46 (MQ=255)

                                                gggAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggataccatt                <  1:1587582‑M1/65‑49 (MQ=255)

                                                gggAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggataccatt                >  1:776173‑M1/1‑17 (MQ=255)

                                                gggAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggataccatt                >  1:2281865‑M1/1‑17 (MQ=255)

                                                  gAACAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggataccattgc              <  1:1337288‑M1/65‑51 (MQ=255)

                                                   aaCAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggataccattgct             <  1:201840‑M1/65‑52 (MQ=255)

                                                    aCAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggatacccttgctg            >  1:1312962‑M1/1‑13 (MQ=255)

                                                    aCAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggataccattgctg            >  1:2516001‑M1/1‑13 (MQ=255)

                                                    aCAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggataccattgctg            >  1:1295960‑M1/1‑13 (MQ=255)

                                                     cAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggataccattgctgc           <  1:2631639‑M1/65‑54 (MQ=255)

                                                     cAGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggataccattgctgc           <  1:1331202‑M1/65‑54 (MQ=255)

                                                      aGGAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggataccattgctgca          <  1:3762976‑M1/65‑55 (MQ=255)

                                                        gAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggataccattgctgcaca        >  1:345252‑M1/1‑9 (MQ=255)

                                                        gAAGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggataccattgctgcaca        >  1:3885750‑M1/1‑9 (MQ=255)

                                                         aaGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggatacgattgctgcacat       >  1:1422123‑M1/1‑8 (MQ=255)

                                                          aGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggataccattgctggaccta      >  1:1484063‑M1/1‑7 (MQ=255)

                                                          aGTGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggataccattgctgcacata      >  1:599011‑M1/1‑7 (MQ=255)

                                                          aGTGCTGacgtcgctgaataccatcgacgatccggcgcgtctggcggatcccattgctgcacata      >  1:564279‑M1/1‑7 (MQ=255)

                                                            tGCTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggaaaccttttctgcacataag    >  1:840986‑M1/1‑5 (MQ=255)

                                                              cTGacgtcgctgaatagcatcgacgatccggcgcgtctggcggataccattgctgcacatatgcc  <  1:3974755‑M1/65‑63 (MQ=255)

ACTTTTCTGCGAAGGCGGAGTATCTGGAGTCGCCGACCATTGATGAGCGGGAACAGGAAGTGCTGGTGCGTACTGCAATCAGCCAGTTCGAAGGCTACATCAAGCTGAACAAAAAAATCCCACCAGA  >  NC_000913/459213‑459339

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