BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,729,373	+C	intergenic (‑195/+247)	rrsG ← / ← clpB	16S ribosomal RNA/ClpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,729,371	1	.	C	98.2%	134.5 / ‑3.2	55	intergenic (‑193/+249)	rrsG/clpB	16S ribosomal RNA/ClpB
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (25/29); total (26/29)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTT‑C‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCCGGCG  >  NC_000913/2729325‑2729441

gCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCt                                            <  2:4692973/76‑1 (MQ=9)

gCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCt                                            >  2:4074038/1‑76 (MQ=9)

  cgcttcgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCttt                                          >  1:825722/1‑76 (MQ=255)

    cttcgcttTTTCTCAGCGTCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTtc                                        <  1:2282755/76‑1 (MQ=11)

    cttcgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTtc                                        >  2:505458/1‑76 (MQ=255)

    cttcgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTT‑CTCGCTACAGAGTCAAGCATTTCTTTGCTTTtc                                        >  1:4475263/1‑76 (MQ=255)

     ttcgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTtct                                       >  2:1286577/1‑76 (MQ=255)

     ttcgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTtct                                       <  1:2377072/76‑1 (MQ=255)

     ttcgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTtct                                       <  1:4046911/76‑1 (MQ=255)

     ctcgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTtct                                       <  1:115030/75‑1 (MQ=255)

      tcgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                                      <  2:346495/76‑1 (MQ=255)

      tcgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                                      >  1:4214385/1‑76 (MQ=255)

      tcgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                                      <  1:637062/76‑1 (MQ=255)

       cgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTctct                                     <  1:1903771/76‑1 (MQ=255)

       cgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTctct                                     <  1:1936257/76‑1 (MQ=255)

        gcttTTTCTCAGAGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTg                                    >  2:1155914/1‑76 (MQ=14)

          tttttCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGtt                                  >  2:4655598/1‑76 (MQ=255)

           ttttCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                                 <  1:4574013/76‑1 (MQ=255)

           ttttCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                                 <  2:2938534/76‑1 (MQ=255)

             ttCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgag                               >  1:4057112/1‑76 (MQ=255)

             ttCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgag                               <  1:4518525/76‑1 (MQ=255)

              tctcAGCGGCGCGGGGCGTGCATAATACGCCTTCC‑CGCTACAGAGCCAAGCATTCCTTTGCTTTTCTCTGTTgaga                              >  2:3426954/1‑76 (MQ=17)

               ctcAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAt                             >  1:4427791/1‑76 (MQ=255)

               ctcAGCGGAGCGGGGTGTACATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAt                             <  1:891082/76‑1 (MQ=17)

                 cAGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtc                           <  1:1123129/76‑1 (MQ=255)

                  aGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtct                          >  2:3285384/1‑76 (MQ=255)

                  aGCGGCGCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtct                          >  1:2618188/1‑76 (MQ=255)

                   gcggcgCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtctc                         <  1:3837197/76‑1 (MQ=255)

                    cggcgCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCa                        <  2:2470777/76‑1 (MQ=255)

                    cggcgCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCa                        <  1:1140523/76‑1 (MQ=255)

                     ggcgCGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg                       >  1:4727476/1‑76 (MQ=255)

                      gcgcGGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTGTCTCTGTTGAGATTCTCAgg                      <  1:1550787/76‑1 (MQ=18)

                       cgcgGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGga                     >  2:3888599/1‑76 (MQ=255)

                        gcgGGGTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGgag                    >  2:2185473/1‑76 (MQ=255)

                          ggggTGTGCATAATACGCCTTCC‑CGCTCCAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGgagac                  >  1:1072251/1‑75 (MQ=21)

                          ggggTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGaa                  >  1:1023950/1‑76 (MQ=255)

                          ggggTGTGCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGaa                  >  2:2632733/1‑76 (MQ=255)

                              tgtgCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAcccc              <  2:528290/76‑1 (MQ=255)

                              tgtgCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAcccc              <  1:1130794/76‑1 (MQ=255)

                               gtgCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgaga                              <  1:4310078/59‑1 (MQ=18)

                               gtgCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgaga                              >  2:4310078/1‑59 (MQ=18)

                               gtgCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCccg             >  1:2718403/1‑76 (MQ=255)

                                 gCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCccgcc           <  2:3828202/76‑1 (MQ=255)

                                 gCATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCccgcc           <  2:1611415/76‑1 (MQ=255)

                                  cATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCccgccg          >  2:320464/1‑76 (MQ=255)

                                  cATAATACGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCccgccg          <  1:1354753/76‑1 (MQ=255)

                                    taataCGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGAc        <  1:1745751/76‑1 (MQ=255)

                                     aataCGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGAcc       >  1:2730721/1‑76 (MQ=255)

                                       tccGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCcg     <  2:4584120/74‑1 (MQ=37)

                                       taCGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCcg     <  1:4027828/76‑1 (MQ=255)

                                         cGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGa         <  1:1857691/70‑1 (MQ=255)

                                         cGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGa         >  2:1857690/1‑70 (MQ=255)

                                         cGCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCcggc   <  2:4427791/76‑1 (MQ=255)

                                          gCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCcggcg  >  1:1459891/1‑76 (MQ=255)

                                          gCCTTCC‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCcggcg  >  2:1520838/1‑76 (MQ=255)

GCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTT‑C‑CGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCCGGCG  >  NC_000913/2729325‑2729441

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