BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,863,452	A→G	29.6%	P373P (CCA→CCG)	hflX →	predicted GTPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,863,452	0	A	G	29.6%	58.1 / 27.7	64	P373P (CCA→CCG)	hflX	predicted GTPase
Reads supporting (aligned to +/- strand): ref base A (21/24); new base G (11/8); total (32/32)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

GCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGA  >  minE/2863389‑2863513

gCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCa                                                          <  1:899222/69‑1 (MQ=255)

gCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCa                                                          <  1:734596/69‑1 (MQ=255)

gCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCa                                                          <  1:514537/69‑1 (MQ=255)

gCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCa                                                          <  1:3166002/69‑1 (MQ=255)

gCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCa                                                          <  1:187785/69‑1 (MQ=255)

 cTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCATCTGCCACCGCa                                                          <  1:1055125/68‑1 (MQ=255)

  tttttCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAgg                                                        <  1:2205923/69‑1 (MQ=255)

  tttttCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAgg                                                        <  1:2172608/69‑1 (MQ=255)

  tttttCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAgg                                                        <  1:176796/69‑1 (MQ=255)

   ttttCAGGCTTTGACGGGGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAgg                                                        <  1:746617/68‑1 (MQ=255)

    tttCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCa                                                          <  1:2385807/65‑1 (MQ=255)

          gCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGt                                                >  1:2417473/1‑69 (MQ=255)

          gCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGt                                                >  1:1862517/1‑69 (MQ=255)

          gCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGt                                                >  1:1851809/1‑69 (MQ=255)

          gCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGt                                                >  1:1351995/1‑69 (MQ=255)

                   gAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       <  1:2690075/69‑1 (MQ=255)

                   gAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       <  1:1172859/69‑1 (MQ=255)

                   gAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       <  1:1955753/69‑1 (MQ=255)

                   gAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       <  1:1345896/69‑1 (MQ=255)

                   gAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       <  1:2317449/69‑1 (MQ=255)

                                gCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCa                                                          <  1:2298353/37‑1 (MQ=255)

                                  gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAg                                                         >  1:480932/1‑36 (MQ=255)

                                  gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       >  1:161298/1‑54 (MQ=255)

                                  gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       >  1:2728995/1‑54 (MQ=255)

                                  gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       >  1:2166537/1‑54 (MQ=255)

                                  gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       >  1:3131461/1‑54 (MQ=255)

                                  gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGt                                       >  1:1370601/1‑54 (MQ=255)

                                  gAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAACGGCGTCTgaga                                          >  1:203482/1‑51 (MQ=255)

                                          gcagcaTACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCaa                 <  1:1430466/68‑1 (MQ=255)

                                          gcagcaTACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCaa                 <  1:2084463/68‑1 (MQ=255)

                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTa                               >  1:350162/1‑51 (MQ=39)

                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:1355996/1‑52 (MQ=39)

                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:1628336/1‑52 (MQ=39)

                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:526767/1‑52 (MQ=39)

                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:2739706/1‑52 (MQ=39)

                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:971609/1‑52 (MQ=39)

                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:2691715/1‑52 (MQ=39)

                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:2437332/1‑52 (MQ=39)

                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:1715876/1‑52 (MQ=39)

                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:2281596/1‑52 (MQ=39)

                                             gcaTACATTGCGTCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTAt                              >  1:1043125/1‑52 (MQ=39)

                                             gcaTACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGa             <  1:2625580/69‑1 (MQ=255)

                                              caTACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGa             <  1:3191018/68‑1 (MQ=255)

                                              caTACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGa             <  1:589439/68‑1 (MQ=255)

                                                tACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGaaaa          <  1:2019209/69‑1 (MQ=255)

                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGaaaaa         >  1:2894471/1‑68 (MQ=255)

                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag        >  1:459042/1‑69 (MQ=255)

                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag        >  1:757783/1‑69 (MQ=255)

                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag        >  1:1312595/1‑69 (MQ=255)

                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag        >  1:2931011/1‑69 (MQ=255)

                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag        >  1:447038/1‑69 (MQ=255)

                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag        >  1:3346106/1‑69 (MQ=255)

                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag        >  1:3317135/1‑69 (MQ=255)

                                                  cATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAag        >  1:318484/1‑69 (MQ=255)

                                                      gCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg    <  1:2848947/69‑1 (MQ=255)

                                                      gCGTCTGCCACAGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtg    >  1:1431696/1‑69 (MQ=255)

                                                         tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga  <  1:2181522/68‑1 (MQ=255)

                                                         tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga  <  1:849140/68‑1 (MQ=255)

                                                         tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga  <  1:803378/68‑1 (MQ=255)

                                                         tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga  <  1:1943021/68‑1 (MQ=255)

                                                         tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga  <  1:634910/68‑1 (MQ=255)

                                                         tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga  <  1:2226142/68‑1 (MQ=255)

                                                         tCTGCCGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga  <  1:3261932/68‑1 (MQ=255)

                                                         tCTGACGCCGCAGGAAGGGCGCCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGtgga  <  1:229727/68‑1 (MQ=255)

GCTTTTTCAGGCTTTGACGGAGCGGCTTTCCGGCGAGGTGGCGCAGCATACATTGCGTCTGCCACCGCAGGAAGGGCGTCTGAGAAGTCGTTTTTATCAGCTTCAGGCAATAGAAAAAGAGTGGA  >  minE/2863389‑2863513

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: `‑`