BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,571,187	G→A	9.8%	L364L (CTG→CTA)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,571,187	0	G	A	9.8%	338.3 / 20.3	112	L364L (CTG→CTA)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (47/54); new base A (7/4); total (54/58)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 118 total aligned reads displayed.

CATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAA  >  minE/2571121‑2571255

  tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                                                 >  1:1434260/1‑70 (MQ=255)

   cTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                                               >  1:2402462/1‑71 (MQ=255)

   cTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGa                                                                >  1:3994567/1‑70 (MQ=255)

    tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                                                              >  1:111777/1‑71 (MQ=255)

    tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                                                              >  1:133312/1‑71 (MQ=255)

    tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                                                              >  1:3220899/1‑71 (MQ=255)

      cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                                                            <  1:3668258/71‑1 (MQ=255)

            ggTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                                                                  <  1:2728731/59‑1 (MQ=255)

             gTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTg                                                     >  1:1095014/1‑71 (MQ=255)

             gTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                     >  1:3945463/1‑71 (MQ=255)

               ggCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                   >  1:545315/1‑71 (MQ=255)

                gCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                                                              >  1:1617535/1‑59 (MQ=255)

                gCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                   <  1:952344/70‑1 (MQ=255)

                 cGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATg                                                 >  1:2814993/1‑71 (MQ=255)

                  ggTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGc                                                <  1:3009714/71‑1 (MQ=255)

                  ggTTTACTCCGACCATTTCAGCCATGTTACGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                     <  1:3312450/66‑1 (MQ=255)

                   gTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCa                                               <  1:3919056/71‑1 (MQ=255)

                    tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                                                                  <  1:1068004/51‑1 (MQ=255)

                    tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                                                                  <  1:4142337/51‑1 (MQ=255)

                     ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                                                            >  1:3380717/1‑56 (MQ=255)

                     ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                              >  1:3169414/1‑70 (MQ=255)

                     ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                              >  1:2338357/1‑70 (MQ=255)

                      tACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGc                                                       >  1:2947690/1‑60 (MQ=255)

                      tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAtt                                             <  1:230781/70‑1 (MQ=255)

                      tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                              >  1:3680293/1‑69 (MQ=255)

                      tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                            >  1:1855529/1‑71 (MQ=255)

                       aCTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                            <  1:248070/70‑1 (MQ=255)

                           cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                                                 >  1:2545533/1‑45 (MQ=255)

                           cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                                                            >  1:2593160/1‑50 (MQ=255)

                           cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCt                                                      >  1:3403223/1‑56 (MQ=255)

                           cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCCGCTAAGAATGACCAGCTGGATGCATTATCtt                                        >  1:2423530/1‑70 (MQ=255)

                               cATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                                               >  1:251483/1‑43 (MQ=255)

                               cATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCa                                               >  1:930626/1‑59 (MQ=255)

                               cATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCa                                               >  1:3276909/1‑59 (MQ=255)

                                aTTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTa                                      <  1:1667984/67‑1 (MQ=255)

                                 tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                                                 <  1:2863215/39‑1 (MQ=255)

                                 tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATg                                                             <  1:2207404/43‑1 (MQ=255)

                                 tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                     <  1:2049599/51‑1 (MQ=255)

                                 tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTgg                                 <  1:679027/71‑1 (MQ=255)

                                   tCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGa                                <  1:2342399/70‑1 (MQ=255)

                                   tCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                               >  1:552544/1‑71 (MQ=255)

                                   tCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                               >  1:746082/1‑71 (MQ=255)

                                     aGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa                             <  1:2335195/71‑1 (MQ=255)

                                       ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                           <  1:3976246/71‑1 (MQ=255)

                                       ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                           <  1:3093849/71‑1 (MQ=255)

                                       ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                           <  1:925931/71‑1 (MQ=255)

                                       ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                           <  1:795084/71‑1 (MQ=255)

                                       ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                           <  1:630160/71‑1 (MQ=255)

                                       ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                           <  1:423590/71‑1 (MQ=255)

                                       ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                           <  1:44487/71‑1 (MQ=255)

                                        cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCtt                                        >  1:4082641/1‑57 (MQ=255)

                                        cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                          <  1:775427/71‑1 (MQ=255)

                                        cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                          <  1:176674/71‑1 (MQ=255)

                                        cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                          <  1:3087011/71‑1 (MQ=255)

                                         aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                     <  1:2191325/43‑1 (MQ=255)

                                         aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                           >  1:3613364/1‑69 (MQ=255)

                                          tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                     <  1:419792/42‑1 (MQ=255)

                                          tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                     <  1:2223298/42‑1 (MQ=255)

                                          tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa                             <  1:824364/66‑1 (MQ=255)

                                          tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAt                            <  1:2066690/67‑1 (MQ=255)

                                          tGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCttt                                       <  1:2911709/56‑1 (MQ=255)

                                          tGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCttt                                       <  1:2338841/56‑1 (MQ=255)

                                           gTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGc                                                       <  1:623439/39‑1 (MQ=255)

                                            tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCt                                                      >  1:485237/1‑39 (MQ=255)

                                            tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                   >  1:148076/1‑42 (MQ=255)

                                            tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                   >  1:2596392/1‑42 (MQ=255)

                                            tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAttt                                   <  1:2604808/58‑1 (MQ=255)

                                            tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggg                      >  1:4022860/1‑71 (MQ=255)

                                            tctCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAt                                     >  1:1880178/3‑56 (MQ=255)

                                            tctCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAt                                     >  1:2445816/3‑56 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                   >  1:899206/1‑40 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                            >  1:1650799/1‑47 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                            >  1:3041325/1‑47 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                            >  1:3582595/1‑47 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                           >  1:3999649/1‑64 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                     <  1:594621/70‑1 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                     <  1:1567541/70‑1 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                     <  1:3833739/70‑1 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                     <  1:1448316/70‑1 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggg                      >  1:1463137/1‑69 (MQ=255)

                                              tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga                         >  1:1228556/1‑66 (MQ=255)

                                              tCGGAGATGCGGAAAGTGCTGCTAAGAATGACCAg                                                        >  1:2288865/1‑35 (MQ=255)

                                               cGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggg                    <  1:1839896/70‑1 (MQ=255)

                                                ggTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                   >  1:3888278/1‑38 (MQ=255)

                                                  tGATGCGGAAAGTGCTGCTAAGAATGACCATCTGGATGCATTATCTTTATTTTTATCa                             <  1:2245708/58‑1 (MQ=255)

                                                   gATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggggagggg               >  1:1510288/1‑71 (MQ=255)

                                                     tGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa                             <  1:2250476/55‑1 (MQ=255)

                                                     tGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgag                        <  1:409267/60‑1 (MQ=255)

                                                      gCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATc                                          >  1:1592872/1‑41 (MQ=255)

                                                      gCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa                             >  1:3914078/1‑54 (MQ=255)

                                                         gAAAGTGCTGCTAAGAATGACCAGCTGGATGCAtt                                             <  1:3005181/35‑1 (MQ=255)

                                                         gAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAt         >  1:3659699/1‑71 (MQ=255)

                                                         gAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAt         >  1:295657/1‑71 (MQ=255)

                                                          aaaGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAt         >  1:2594483/1‑70 (MQ=255)

                                                              tgctgcTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATc                              <  1:1206929/45‑1 (MQ=255)

                                                              tgctgcTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATg     >  1:263586/1‑70 (MQ=255)

                                                               gctgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggg                    >  1:695984/1‑54 (MQ=255)

                                                                ctgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGaa  >  1:1196068/1‑71 (MQ=255)

                                                                ctgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGaa  >  1:3857652/1‑71 (MQ=255)

                                                                 tgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                          <  1:3705773/46‑1 (MQ=255)

CATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAA  >  minE/2571121‑2571255

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