BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,867	G→A	45.3%	noncoding (1228/2905 nt)	rrlA ←	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,359,867	0	G	A	45.3%	68.9 / 122.6	95	noncoding (1228/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base G (22/30); new base A (31/12); total (53/42)
Fisher's exact test for biased strand distribution p-value = 4.11e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 108 total aligned reads displayed.

TTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCA  >  minE/2359821‑2359936

ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                               >  1:3593259/1‑71 (MQ=34)

ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                               >  1:2972836/1‑71 (MQ=34)

ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                               >  1:2879441/1‑71 (MQ=34)

     gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                             >  1:476828/1‑68 (MQ=35)

     gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                           <  1:172415/70‑1 (MQ=35)

     gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                          <  1:1532954/71‑1 (MQ=34)

       tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                        >  1:1221893/1‑71 (MQ=255)

       tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                        >  1:2770200/1‑71 (MQ=255)

       tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                        >  1:5197373/1‑71 (MQ=255)

       tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGa                                         >  1:355594/1‑70 (MQ=255)

        aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTcac                                                    <  1:4116475/58‑1 (MQ=33)

        aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTcac                                                    <  1:3671486/58‑1 (MQ=33)

        aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        <  1:3074979/70‑1 (MQ=35)

        aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        <  1:2224595/70‑1 (MQ=35)

           tATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                             <  1:4822703/62‑1 (MQ=35)

           tATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                             <  1:2568997/62‑1 (MQ=35)

           tATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                             <  1:3633800/62‑1 (MQ=35)

             tGTCAGGATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                  >  1:669856/1‑71 (MQ=21)

             tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                  >  1:2942111/1‑71 (MQ=34)

             tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                  >  1:2695897/1‑71 (MQ=34)

             tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                  >  1:3058184/1‑71 (MQ=34)

                 aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                               <  1:2449715/70‑1 (MQ=255)

                 aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                               <  1:3206644/70‑1 (MQ=255)

                 aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTa                              >  1:3989590/1‑71 (MQ=255)

                   cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                            >  1:967328/1‑71 (MQ=255)

                     ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCca                          <  1:502430/71‑1 (MQ=255)

                       cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                         >  1:2962940/1‑70 (MQ=35)

                           cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                    <  1:661449/71‑1 (MQ=35)

                           cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                    >  1:2943447/1‑70 (MQ=255)

                           cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                    <  1:3460941/71‑2 (MQ=255)

                           cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                    <  1:8550/71‑2 (MQ=255)

                           cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                    >  1:785345/1‑70 (MQ=255)

                           cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                    <  1:130031/71‑2 (MQ=255)

                           cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                    <  1:440670/71‑2 (MQ=255)

                           cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                    <  1:1654902/71‑2 (MQ=255)

                           cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                    <  1:1759132/71‑2 (MQ=255)

                           cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                    <  1:2063386/71‑2 (MQ=255)

                            ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCtt                                              <  1:3720816/44‑1 (MQ=11)

                            ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCtt                                              <  1:400776/44‑1 (MQ=11)

                              cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                  >  1:743714/1‑70 (MQ=40)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:3793280/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:549057/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:3799094/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:5020869/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:3614209/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:3930489/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:4423917/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:88289/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:1168903/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:1129106/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:1814784/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:2601325/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:277183/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:2875357/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:295439/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:1632501/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:3101587/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:126803/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTca                                                     >  1:117889/1‑35 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTa                              >  1:3180386/1‑58 (MQ=11)

                              cTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTa                              >  1:2324336/1‑58 (MQ=11)

                                  tACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATa                <  1:3050231/68‑1 (MQ=40)

                                      tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                          <  1:1611612/38‑1 (MQ=40)

                                      tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                          <  1:1649474/38‑1 (MQ=40)

                                      tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                          <  1:525227/38‑1 (MQ=40)

                                      tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                          <  1:1215444/38‑1 (MQ=40)

                                      tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                          <  1:3843980/38‑1 (MQ=40)

                                      tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                          <  1:3850218/38‑1 (MQ=40)

                                      tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                         <  1:1492212/39‑1 (MQ=40)

                                      tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                         <  1:4220735/39‑1 (MQ=40)

                                      tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                         <  1:4974850/39‑1 (MQ=40)

                                      tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                         <  1:4226893/39‑1 (MQ=40)

                                      tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                         <  1:2368094/39‑1 (MQ=40)

                                      tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                         <  1:3372018/39‑1 (MQ=40)

                                      tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                         <  1:342063/39‑1 (MQ=40)

                                      tCCAGCATGACTCACAGCACACCTTCACAGGCTTACAGa                                         <  1:3443849/39‑1 (MQ=255)

                                       ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                      >  1:2796178/1‑57 (MQ=35)

                                       ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                      >  1:52195/1‑57 (MQ=35)

                                        cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:4646475/1‑38 (MQ=40)

                                        cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:2036929/1‑38 (MQ=40)

                                        cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:727559/1‑38 (MQ=40)

                                        cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:3492344/1‑38 (MQ=40)

                                        cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:4988248/1‑38 (MQ=40)

                                        cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:1953393/1‑38 (MQ=40)

                                        cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                        >  1:4032354/1‑38 (MQ=40)

                                        cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                         >  1:2959700/1‑37 (MQ=40)

                                        cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGc         <  1:2456418/69‑1 (MQ=255)

                                          gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc     >  1:3661244/1‑71 (MQ=255)

                                              gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt        <  1:3410515/64‑1 (MQ=255)

                                              gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt        <  1:3277032/64‑1 (MQ=255)

                                              gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt        <  1:2001148/64‑1 (MQ=255)

                                              gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt        <  1:3976395/64‑1 (MQ=255)

                                              gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt        <  1:4096692/64‑1 (MQ=255)

                                              gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt        <  1:2479820/64‑1 (MQ=255)

                                              gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt        <  1:4918407/64‑1 (MQ=255)

                                              gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt        <  1:214763/64‑1 (MQ=255)

                                              gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt        <  1:4536685/64‑1 (MQ=255)

                                              gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt        <  1:4429718/64‑1 (MQ=255)

                                              gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt        <  1:2244141/64‑1 (MQ=255)

                                              gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt        <  1:3344076/64‑1 (MQ=255)

TTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCA  >  minE/2359821‑2359936

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