BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AP009048	4,357,331	T→A	5.3%	F24I (TTC→ATC)	yjdK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AP009048	4,357,331	0	T	A	5.3%	384.6 / 10.2	132	F24I (TTC→ATC)	yjdK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (71/54); new base A (5/2); total (76/56)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

Only 100 of 151 total aligned reads displayed.

TATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATGGGATAC  >  AP009048/4357261‑4357400

tgtGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAgtgt                                                                       >  1:7157440/3‑71 (MQ=255)

 atGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAgtgt                                                                       <  1:532000/70‑1 (MQ=255)

  tGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTc                                                                     >  1:5124297/1‑71 (MQ=255)

         tAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagat                                                              <  1:5314986/71‑1 (MQ=255)

         tAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagat                                                              >  1:6243278/1‑71 (MQ=255)

          aaaaaCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagat                                                              >  1:2385946/1‑70 (MQ=255)

          aaaaaCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagat                                                              >  1:1726273/1‑70 (MQ=255)

          aaaaaCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagat                                                              >  1:7678018/1‑70 (MQ=255)

           aaaaCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagatta                                                            <  1:5116879/71‑1 (MQ=255)

           aaaaCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagatta                                                            <  1:4168918/71‑1 (MQ=255)

           aaaaCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagatt                                                             <  1:350967/70‑1 (MQ=255)

            aaaCAAGTTCAATAGTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagatt                                                             >  1:2964946/1‑69 (MQ=255)

            aaaCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagatt                                                             >  1:1577165/1‑69 (MQ=255)

            aaaCAAGTTCAATACTTATGTTGTTTCTTTTAATTATCCATCATCTTATTCCTCAGTGTTCttaagatt                                                             >  1:3433550/1‑69 (MQ=255)

             aaCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCtt                                                                   <  1:451043/62‑1 (MQ=255)

              aCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagatta                                                            <  1:5936349/68‑1 (MQ=255)

              aCAAGTTCAATACTTATGTTGTTTCTTTTGATTATCAATCATCTTATTCCTCAGTGTTCttaagatt                                                             >  1:2442710/1‑67 (MQ=255)

               cAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagatt                                                             >  1:1008557/1‑66 (MQ=255)

               cAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCttaagat                                                              >  1:3604005/1‑65 (MQ=255)

               cAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTtaagattaagat                                                        <  1:2739540/71‑1 (MQ=255)

               cAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATATCATACATCAATGATCttaatcttaagaa                                                        <  1:4302477/71‑2 (MQ=17)

                 aGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCa                                                      >  1:5653615/1‑71 (MQ=255)

                  gTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCAt                                                     >  1:4873046/1‑71 (MQ=255)

                  gTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCAt                                                     >  1:7438407/1‑71 (MQ=255)

                  gTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCAt                                                     >  1:7447532/1‑71 (MQ=255)

                     cAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATc                                                       <  1:835483/66‑1 (MQ=255)

                      aaTACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGa                                                  >  1:1814116/1‑70 (MQ=255)

                      aaTACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGAt                                                 >  1:4211756/1‑71 (MQ=255)

                             aTGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtatga                                           <  1:2054130/70‑1 (MQ=255)

                             aTGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtatga                                           <  1:574675/70‑1 (MQ=255)

                              tgttgtTTCTTTTGATTATCCATCATCTTAATCCTCAGTGTTCTTAAGATTAAGATCATTGATGtatgat                                          <  1:6026784/70‑1 (MQ=255)

                               gttgttTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTGTTAAGATTAAGATCa                                                      >  1:4364414/1‑57 (MQ=255)

                                ttgttTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtatgatgt                                        <  1:7076172/70‑3 (MQ=255)

                                ttgttTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtatgatat                                        <  1:2866252/70‑1 (MQ=255)

                                ttgttTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtatgatat                                        <  1:5499642/70‑1 (MQ=255)

                                  gttTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtatgat                                          <  1:6679447/66‑1 (MQ=255)

                                  gttTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGaa                                     >  1:7557179/1‑71 (MQ=255)

                                  gttTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGaa                                     >  1:7554384/1‑71 (MQ=255)

                                    ttCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAAtt                                   >  1:3957689/1‑71 (MQ=255)

                                    ttCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAAt                                    <  1:2330545/70‑1 (MQ=255)

                                      cTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtat                                             <  1:7379275/59‑1 (MQ=255)

                                        tttGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGAt                                                 >  1:6130025/1‑53 (MQ=255)

                                        tttGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTtctc                               >  1:6785309/1‑71 (MQ=255)

                                         ttGATTATCCATCATCTTATTCCTCAGTGTTCTtaagattaagat                                                        >  1:7197851/1‑45 (MQ=255)

                                         ttGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGtat                                             >  1:2876694/1‑56 (MQ=255)

                                         ttGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTtc                                 <  1:3926448/68‑1 (MQ=255)

                                          ttaTTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTtct                                <  1:5387204/66‑1 (MQ=255)

                                             ttATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAAttt                                  >  1:7874406/1‑63 (MQ=255)

                                             ttATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTa                          >  1:2451424/1‑71 (MQ=255)

                                               aTCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTAt                         <  1:770864/70‑1 (MQ=255)

                                               aTCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTAt                         <  1:6327128/70‑1 (MQ=255)

                                               aTCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATc                        >  1:4003922/1‑71 (MQ=255)

                                               aTCCATCATCTTATTCCTCAATGATCTTAATCTTAAGAACac                                                     >  1:3457465/1‑41 (MQ=2)

                                               aTCCATCATCTTATTCCTCAATGATCTTAATCTTAAGAACac                                                     >  1:5480597/1‑41 (MQ=2)

                                               aTCCATCATCTTATTCCTCAATGATCTTAATCTTAAGAACac                                                     >  1:3747075/1‑41 (MQ=2)

                                                tCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATc                        <  1:6803998/70‑1 (MQ=255)

                                                tCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATc                        >  1:1230252/1‑70 (MQ=255)

                                                tCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCg                       <  1:2443517/71‑1 (MQ=255)

                                                    tcatcTTATTCCTCAGTGTTCTTAAGATTAAGATCa                                                      <  1:5980845/36‑1 (MQ=255)

                                                    tcatcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAAt                                    <  1:5413538/54‑1 (MQ=255)

                                                    tcatcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTtctc                               <  1:36111/59‑1 (MQ=255)

                                                    tcatcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTtctc                               <  1:7923413/59‑1 (MQ=255)

                                                     catcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCt                  >  1:5813902/1‑71 (MQ=255)

                                                      atcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTctcctc                            >  1:97069/1‑60 (MQ=255)

                                                      atcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTctcctc                            >  1:3554774/1‑60 (MQ=255)

                                                      atcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCtg                 <  1:6261208/71‑1 (MQ=255)

                                                      atcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCt                  >  1:2059021/1‑70 (MQ=255)

                                                      atcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCt                  <  1:6738139/70‑1 (MQ=255)

                                                       tcTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTg                     >  1:6322325/1‑66 (MQ=255)

                                                        cTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGt                      <  1:6889891/64‑1 (MQ=255)

                                                         ttATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCtgat               >  1:7842091/1‑70 (MQ=255)

                                                         ttATACCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCtgat               >  1:3909232/1‑70 (MQ=255)

                                                          tATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATTTGAAt                                    >  1:6611989/1‑48 (MQ=255)

                                                          tATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAAt                                    >  1:1919660/1‑48 (MQ=255)

                                                          tATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAAt                                    >  1:4505927/1‑48 (MQ=255)

                                                          tATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTgg                    >  1:275983/1‑64 (MQ=255)

                                                             tCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAata          >  1:5578614/1‑71 (MQ=255)

                                                              ccTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAata          <  1:2408991/70‑1 (MQ=255)

                                                              ccTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAata          <  1:4846516/70‑1 (MQ=255)

                                                              ccTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAata          <  1:1002387/70‑1 (MQ=255)

                                                              ccTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAata          <  1:5629002/70‑1 (MQ=255)

                                                              ccTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTACTCTATCGt                      <  1:257423/58‑1 (MQ=255)

                                                                tCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCtgatga             <  1:5465597/65‑1 (MQ=255)

                                                                 cAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTtct                                >  1:2135179/1‑45 (MQ=255)

                                                                 cAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCtga                >  1:1816914/1‑61 (MQ=255)

                                                                 cAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCtga                >  1:4169149/1‑61 (MQ=255)

                                                                 cAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCtga                >  1:4225418/1‑61 (MQ=255)

                                                                 cAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATgg       >  1:5916984/1‑70 (MQ=255)

                                                                 cAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATg        >  1:1854131/1‑69 (MQ=255)

                                                                  aGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATgg       <  1:3108018/69‑1 (MQ=255)

                                                                  aGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATgg       <  1:70560/69‑1 (MQ=255)

                                                                   gtgtTCTTAAGATTAAGATCATTGATGTATGATATGAATTtct                                <  1:3645466/43‑1 (MQ=255)

                                                                   gtgtTCTTAAGATTAAGATCATTGATGTATGATATGAATTtct                                <  1:5032439/43‑1 (MQ=255)

                                                                      ttCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCg                       >  1:6926135/1‑49 (MQ=255)

                                                                      ttCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCg                       >  1:6140624/1‑49 (MQ=255)

                                                                      ttCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCg                       >  1:1768227/1‑49 (MQ=255)

                                                                      ttCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCg                       >  1:1740162/1‑49 (MQ=255)

                                                                      ttCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCg                       >  1:1262918/1‑49 (MQ=255)

                                                                      ttCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATGGGATAc  >  1:7355053/1‑70 (MQ=255)

                                                                      ttCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATGGGATAc  <  1:6910833/70‑1 (MQ=255)

TATGGAAGGTAAAAACAAGTTCAATACTTATGTTGTTTCTTTTGATTATCCATCATCTTATTCCTCAGTGTTCTTAAGATTAAGATCATTGATGTATGATATGAATTTCTCCTCTATCGTGGCTGATGAATATGGGATAC  >  AP009048/4357261‑4357400

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