BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,905,908	G→A	40.8%	T5I (ACC→ATC)	ptr ←	protease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,905,908	0	G	A	40.8%	29.2 / 48.1	71	T5I (ACC→ATC)	ptr	protease III
Reads supporting (aligned to +/- strand): ref base G (25/17); new base A (15/14); total (40/31)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

CCGTTTCTGCCTGACTTAAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCTGTCCAGATTCAAC  >  minE/1905846‑1905964

ccGTTTGTGCCTGACTTAAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGtgctgc                                                    >  1:436976/1‑69 (MQ=255)

ccGTTTCTGCCTGACTTAAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGtgctgc                                                    >  1:202088/1‑69 (MQ=255)

ccGTTTCTGCCTGACTTAAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGtgctgc                                                    >  1:961937/1‑69 (MQ=255)

ccGTTTCTGCCTGACTTAAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGtgctgc                                                    >  1:2203314/1‑69 (MQ=255)

ccGTTTCTGCCTGACTTAAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGtgctgc                                                    >  1:2565975/1‑69 (MQ=255)

ccGTTTCTGCCTGACTTAAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGtgctgc                                                    >  1:564972/1‑69 (MQ=255)

ccGTTTCTGCCTGACTTAAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGtgctgc                                                    >  1:1425806/1‑69 (MQ=255)

ccGTTTCTGCCTGACTTAAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGtgctgc                                                    >  1:3423421/1‑69 (MQ=255)

ccGTTTCTGCCTGACTTAAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGtgctgc                                                    >  1:2691983/1‑69 (MQ=255)

ccGTTTCTGCCTGACTTAAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGt                                                         >  1:1650658/1‑64 (MQ=255)

   tttCTGCCTGACTTAAGGGTGCCCAGAGGGCAACAAACAACAATAATGCTTTGAGCCAGATGCTGCggg                                                 <  1:1168407/69‑1 (MQ=255)

   tttCTGCCTGACTTAAGGGTGCCCAGAGGGCAACAAACAACAATAATGCTTTGAACCAGATGCTGCggg                                                 <  1:397922/69‑1 (MQ=255)

   tttCTGCCTGACTTAAGGGTGCCCAGAGGGCAACAAACAACAATAATGCTTTGAACCAGATGCTGCggg                                                 <  1:2076011/69‑1 (MQ=255)

   tttCTGCCTGACTTAAGGGTGCCCAGAGGGCAACAAACAACAATAATGCTTTGAACCAGATGCTGCggg                                                 <  1:2662695/69‑1 (MQ=255)

   tttCTGCCTGACTTAAGGGTGCCCAGAGGGCAACAAACAACAATAATGCTTTGAACCAGATGCTGCggg                                                 <  1:3325940/69‑1 (MQ=255)

   tttCTGCCTGACTTAAGGGTGCCCAGAGGGCAACAAACAACAATAATGCTTTGAACCAGATGCTGCggg                                                 <  1:9896/69‑1 (MQ=255)

   tttCTGCCTGACTTAAGGGTGCCCAGAGGGCAACAAACAACAATAATGCTTTGAACCAGATGCTGCggg                                                 <  1:589136/69‑1 (MQ=255)

   tttCTGCCTGACTTAAGGGTGCCCAGAGGGCAACAAACAACAATAATGCTTTGAACCAGATGCTGCggg                                                 <  1:595682/69‑1 (MQ=255)

   tttCTGCCTGACTTAAGGGTGCCCAGAGGGCAACAAACAACAATAATGCTTTGAACCAGATGCTGCggg                                                 <  1:837943/69‑1 (MQ=255)

   tttCTGCCTGACTTAAGGGTGCCCAGAGGGCAACAAACAACAATAATGCTTTGAACCAGATGCTGCggg                                                 <  1:973218/69‑1 (MQ=255)

   tttCTGCCTGACTTAAGGGTGCCCAGAGGGCAACAAACAACAATAATGCTTTGAACCAGATGCTGCggg                                                 <  1:1986155/69‑1 (MQ=255)

   tttCTGCCTGACTTAAGGGTGCCCAGAGGGCAACAAACAACAATAATGCTTTGAACCAGATGCTGCggg                                                 <  1:1927547/69‑1 (MQ=255)

   tttCTGCCTGACTTAAGGGTGCCCAGAGGGCAACAAACAACAATAATCCTTTGAACCAGATGCTGCggg                                                 <  1:2171277/69‑1 (MQ=255)

       tGCCTGACTTAAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCa                                              <  1:3136053/68‑1 (MQ=255)

       tGCCTGACTTAAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCa                                              <  1:205179/68‑1 (MQ=255)

       tGCCTGACTTAAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCa                                              <  1:1450344/68‑1 (MQ=255)

       tGCCTGACTTAAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCa                                              <  1:1404594/68‑1 (MQ=255)

         ccTGACTTAAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCa                                              <  1:2774672/66‑1 (MQ=255)

                   gggTGCCCAGAGGGCAACAAACAACAATAATGCTTTGAACCAGATGCTGCggg                                                 <  1:732272/53‑1 (MQ=255)

                     gTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCTTTCACGGACCTCATa                               >  1:507010/1‑69 (MQ=255)

                     gTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCAt                                >  1:2884860/1‑68 (MQ=255)

                     gTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCAt                                >  1:2884283/1‑68 (MQ=255)

                     gTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATa                               >  1:3256562/1‑69 (MQ=255)

                     gTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATa                               >  1:3252822/1‑69 (MQ=255)

                      tGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCTGGTGCTGCGGGGCATTCACGGACCTCATa                               >  1:2502423/1‑68 (MQ=255)

                      tGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATa                               >  1:2786012/1‑68 (MQ=255)

                      tGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATa                               >  1:2680103/1‑68 (MQ=255)

                                 cAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAAtt                   >  1:1469555/1‑69 (MQ=255)

                                 cAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAAt                    >  1:786222/1‑68 (MQ=255)

                                 cAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAAt                    >  1:835738/1‑68 (MQ=255)

                                     tAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCt               <  1:2369260/69‑1 (MQ=255)

                                     tAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCt               <  1:1584239/69‑1 (MQ=255)

                                      aacaacAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCt               <  1:1685998/68‑1 (MQ=255)

                                      aacaacAATAATGCTTTGAACCAGATGCTGCGGGGCATTCACg                                        >  1:2322558/1‑43 (MQ=255)

                                      aacaacAATAATGCTTTGAACCAGATGCTGCGGGGCATTCACGGACCTc                                  >  1:1309796/1‑49 (MQ=255)

                                      aacaacAATAATGCTTTGAACCAGATGCTGCGGGGCATTCACGGACCTc                                  >  1:817701/1‑49 (MQ=255)

                                      aacaacAATAATGCTTTGAACCAGATGCTGCGGGGCATTCACGGACCTc                                  >  1:1599192/1‑49 (MQ=255)

                                      aacaacAATAATGCTTTGAACCAGATGCTGCGGGGCATTCACGGACCTc                                  >  1:3277145/1‑49 (MQ=255)

                                      aacaacAATAATGCTTTGAACCAGATGCTGCGGGGCATTCACGGACCTc                                  >  1:1122520/1‑49 (MQ=255)

                                      aacaacAATAATGCTTTGAACCAGATGCTGCGGGGCATTCACGGACCTc                                  >  1:1051713/1‑49 (MQ=255)

                                      aacaacAATAATGCTTTGAACCAGATGCTGCGGGGCATTCACGGACCTc                                  >  1:2945892/1‑49 (MQ=255)

                                      aacaacAATAATGCTTTGAACCAGATGCTGCGGGGCATTCACGGACCTc                                  >  1:3016937/1‑49 (MQ=255)

                                      aacaacAATAATGCTTTGAACCAGATGCTGCGGGGCATTCACGGACCTc                                  >  1:202608/1‑49 (MQ=255)

                                      aacaacAATAATGCTTTGAACCAGATGCTGCGGGGCATTCACGGACCTc                                  >  1:2048343/1‑49 (MQ=255)

                                      aacaacAATAATGCTTTGAACCAGATGCTGCGGGGCATTCACGGACCTc                                  >  1:2055481/1‑49 (MQ=255)

                                      aacaacAATAATGCTTTGAACCAGATGCTGCGGGGCATTCACGGACCTc                                  >  1:1743863/1‑49 (MQ=255)

                                      aacaacAATAATGCTTTGAACCAGATGCTGCGGGGCATTCACGGACCTc                                  >  1:2629950/1‑49 (MQ=255)

                                      aacaacAATAATGCTTTGAACCAGATGCTGCGGGGCATTCACGGACCTc                                  >  1:1712375/1‑49 (MQ=255)

                                             ataatGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATGAACATAATTTTCTGTCCAGa        >  1:2732410/1‑68 (MQ=255)

                                             ataatGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCTGTCCAGa        >  1:1961491/1‑68 (MQ=255)

                                             ataatGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCTGTCCAGa        >  1:2370816/1‑68 (MQ=255)

                                             ataatGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCTGTCCAGa        >  1:1405844/1‑68 (MQ=255)

                                                  gCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCTGTCCAGATTCAac  <  1:2551674/69‑1 (MQ=255)

                                                  gCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCTGTCCAGATTCAac  <  1:2536346/69‑1 (MQ=255)

                                                  gCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCTGTCCAGATTCAac  <  1:1648753/69‑1 (MQ=255)

                                                  gCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCTGTCCAGATTCAac  <  1:3412770/69‑1 (MQ=255)

                                                  gCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCTGTCCAGATTCAac  <  1:3354454/69‑1 (MQ=255)

                                                  gCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCTGTCCAGATTCAac  <  1:1233197/69‑1 (MQ=255)

                                                  gCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCTGTCCAGATTCAac  <  1:1170183/69‑1 (MQ=255)

                                                  gCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCTGTCCAGATTCAac  <  1:3351331/69‑1 (MQ=255)

                                                  gCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCTGTCCAGATTCAac  <  1:3055945/69‑1 (MQ=255)

CCGTTTCTGCCTGACTTAAGGGTGCCCAAAGGGCAACTAACAACAATAATGCTTTGAACCAGGTGCTGCGGGGCATTCACGGACCTCATAATCAACTTAATTTTCTGTCCAGATTCAAC  >  minE/1905846‑1905964

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