BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,756,138	Δ1 bp	intergenic (+183/‑21)	ssrA → / → intA	tmRNA/CP4‑57 prophage; integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,756,137	0	C	.	100.0%	250.7 / NA	51	intergenic (+182/‑22)	ssrA/intA	tmRNA/CP4‑57 prophage; integrase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base . (25/26); total (25/26)

TAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTT  >  NC_000913/2756090‑2756181

tAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGc                                              >  1:2452110/1‑48 (MQ=255)

tAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGc                                              <  1:837822/48‑1 (MQ=255)

 aGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGct                                             <  1:961674/48‑2 (MQ=255)

 aGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGct                                             >  1:902418/1‑47 (MQ=255)

 aGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGct                                             >  1:1226917/1‑47 (MQ=255)

    cccATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAAc                                        >  1:727176/1‑49 (MQ=255)

    cccATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAAc                                        >  1:1766514/1‑49 (MQ=255)

     ccATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACg                                       <  1:1879527/49‑1 (MQ=255)

     ccATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACg                                       <  1:188001/49‑1 (MQ=255)

     ccATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACg                                       >  1:1431605/1‑49 (MQ=255)

      cATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACg                                       >  1:2038999/1‑48 (MQ=255)

       aTTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCa                                     <  1:2023583/49‑1 (MQ=255)

        ttGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAg                                    >  1:585431/1‑49 (MQ=255)

          gATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAgg                                   <  1:538268/48‑1 (MQ=255)

          gATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAgg                                   >  1:2628341/1‑48 (MQ=255)

           aTAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAgg                                   >  1:48455/1‑47 (MQ=255)

           aTAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAgg                                   >  1:597034/1‑47 (MQ=255)

           aTAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAgg                                   >  1:2307342/1‑47 (MQ=255)

           aTAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAgg                                   >  1:2069358/1‑47 (MQ=255)

           aTAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAgg                                   >  1:2252999/1‑47 (MQ=255)

           aTAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGga                                 <  1:2052266/49‑1 (MQ=255)

           aTAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGga                                 <  1:2511730/49‑1 (MQ=255)

             aGGCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGgaga                               >  1:855112/1‑49 (MQ=255)

              ggCCCAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGAc                              <  1:25863/49‑1 (MQ=255)

                cccAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTc                            >  1:2461809/1‑49 (MQ=255)

                 ccAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTcc                           <  1:2326918/49‑1 (MQ=255)

                 ccAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTcc                           >  1:1451727/1‑49 (MQ=255)

                  cAACGAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTccc                          >  1:1254622/1‑49 (MQ=255)

                      gAAAAGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcc                      >  1:1545440‑M1/1‑47 (MQ=255)

                       aaaaGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcct                     <  1:1123471‑M1/49‑4 (MQ=255)

                        aaaGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctt                    >  1:489601‑M1/1‑45 (MQ=255)

                         aaGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAccttt                   <  1:1428397‑M1/49‑6 (MQ=255)

                          aGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttg                  <  1:2728860‑M1/49‑7 (MQ=255)

                          aGCTCTATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttg                  <  1:2837932‑M1/49‑7 (MQ=255)

                            ctctATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaa                <  1:96056‑M1/49‑9 (MQ=255)

                             tctATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaa                >  1:679087‑M1/1‑40 (MQ=255)

                             tctATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaa                >  1:1574137‑M1/1‑40 (MQ=255)

                             tctATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaa                >  1:2428496‑M1/1‑40 (MQ=255)

                             tctATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaa                >  1:315216‑M1/1‑40 (MQ=255)

                             tctATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaa                >  1:1166318‑M1/1‑40 (MQ=255)

                               tATTGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaa               <  1:319417‑M1/47‑10 (MQ=255)

                                 ttGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacag           <  1:398348‑M1/49‑14 (MQ=255)

                                 ttGTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaa              <  1:892805‑M1/46‑11 (MQ=255)

                                   gTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacagga         >  1:2714389‑M1/1‑34 (MQ=255)

                                   gTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacagg          <  1:1264129‑M1/48‑15 (MQ=255)

                                   gTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacagg          <  1:918422‑M1/48‑15 (MQ=255)

                                   gTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacagg          <  1:1747073‑M1/48‑15 (MQ=255)

                                   gTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacagg          <  1:1633945‑M1/48‑15 (MQ=255)

                                   gTTTACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacagg          <  1:2440648‑M1/48‑15 (MQ=255)

                                    tttACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggat        <  1:528186‑M1/49‑17 (MQ=255)

                                    tttACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggat        <  1:1906189‑M1/49‑17 (MQ=255)

                                    tttACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggat        >  1:901551‑M1/1‑33 (MQ=255)

                                      tACGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatg       >  1:2611210‑M1/1‑31 (MQ=255)

                                        cGTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtag    <  1:2032052‑M1/49‑21 (MQ=255)

                                         gTTGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc   <  1:2073546‑M1/49‑22 (MQ=255)

                                          ttGGG‑CTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcg  <  1:2568447‑M1/49‑23 (MQ=255)

TAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTT  >  NC_000913/2756090‑2756181

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