BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,905,025	G→A	17.2%	N299N (AAC→AAT)	ptr ←	protease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,905,025	0	G	A	17.2%	180.7 / 21.0	109	N299N (AAC→AAT)	ptr	protease III
Reads supporting (aligned to +/- strand): ref base G (69/21); new base A (4/15); total (73/36)
Fisher's exact test for biased strand distribution p-value = 8.18e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 114 total aligned reads displayed.

CCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGAT  >  minE/1904960‑1905090

ccTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAgtt                                                                 >  1:3273970/1‑68 (MQ=255)

ccTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAgtt                                                                 >  1:3184857/1‑68 (MQ=255)

ccTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAgtt                                                                 >  1:500736/1‑68 (MQ=255)

ccTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAgtt                                                                 >  1:2250669/1‑68 (MQ=255)

           aTTACCAATCAGATAGGTAATCAATTCATCCGTTTTACTACGGAACTTCGCTGAATTGTTATCGATGCg                                                     <  1:2603580/69‑1 (MQ=255)

           aTTACCAATCAGATAGGTAATCAATTCATCCGTTTTACTACGGAACTTCGCTGAATTGTTATCGATGCg                                                     <  1:841012/69‑1 (MQ=255)

           aTTACCAATCAGATAGGTAATCAATTCATCCGTTTTACTACGGAACTTCGCTGAATTGTTATCGATGCg                                                     <  1:1945239/69‑1 (MQ=255)

           aTTACCAATCAGATAGGTAATCAATTCATCCGTTTTACTACGGAACTTCGCTGAATTGTTATCGATGCg                                                     <  1:683/69‑1 (MQ=255)

           aTTACCAATCAGATAGGTAATCAATTCATCCGTTTTACTACGGAACTTCGCTGAATTGTTATCGATGCg                                                     <  1:587440/69‑1 (MQ=255)

           aTTACCAATCAGATAGGTAATCAATTCATCCGTTTTACTACGGAACTTCGCTGAATTGTTATCGATGCg                                                     <  1:477385/69‑1 (MQ=255)

           aTTACCAATCAGATAGGTAATCAATTCATCCGTTTTACTACGGAACTTCGCTGAATTGTTATCGATGCg                                                     <  1:2347433/69‑1 (MQ=255)

           aTTACCAATCAGATAGGTAATCAATTCATCCGTTTTACTACGGAACTTCGCTGAATTGTTATCGATGCg                                                     <  1:38100/69‑1 (MQ=255)

           aTTACCAATCAGATAGGTAATCAATTCATCCGTTTTACTACGGAACTTCGCTGAATTGTTATCGATGCg                                                     <  1:2638479/69‑1 (MQ=255)

           aTTACCAATCAGATAGGTAATCAATTCATCCGTTTTACTACGGAACTTCGCTGAATTGTTATCGATGCg                                                     <  1:293163/69‑1 (MQ=255)

           aTTACCAATCAGATAGGTAATCAATTCATCCGTTTTACTACGGAACTTCGCTGAATTGTTATCGATGCg                                                     <  1:1147604/69‑1 (MQ=255)

           aTTACCAATCAGATAGGTAATCAATTCATCCGTTTTACTACGGAACTTCGCTGAATTGTTATCGATGCg                                                     <  1:3217929/69‑1 (MQ=255)

           aTTACCAATCAGATAGGTAATCAATTCATCCGTTTTACTACGGAACTTCGCTGAATTGTTATCGATGCg                                                     <  1:3326966/69‑1 (MQ=255)

            ttGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGa                                                    <  1:2118196/69‑1 (MQ=255)

            ttGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGa                                                    <  1:2342397/69‑1 (MQ=255)

            ttGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGa                                                    <  1:363257/69‑1 (MQ=255)

            ttGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGa                                                    <  1:3108365/69‑1 (MQ=255)

            ttGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGa                                                    <  1:3331924/69‑1 (MQ=255)

            ttGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGa                                                    <  1:1320329/69‑1 (MQ=255)

            ttGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGa                                                    <  1:955960/69‑1 (MQ=255)

             tGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGaa                                                   >  1:1907797/1‑69 (MQ=255)

             tGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGaa                                                   >  1:1351235/1‑69 (MQ=255)

             tGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGaa                                                   >  1:2874322/1‑69 (MQ=255)

             tGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGaa                                                   >  1:1462871/1‑69 (MQ=255)

             tGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGa                                                    <  1:3190507/68‑1 (MQ=255)

             tGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGa                                                    >  1:1554914/1‑68 (MQ=255)

             tGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGa                                                    <  1:1302998/68‑1 (MQ=255)

             tGCCAATCAGATAGGTAATCAATTCATAGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGaa                                                   >  1:1397027/1‑69 (MQ=255)

              gCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGaaa                                                  >  1:16234/1‑69 (MQ=255)

              gCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGaaa                                                  >  1:852226/1‑69 (MQ=255)

              gCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGaaa                                                  >  1:3382673/1‑69 (MQ=255)

              gCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGaaa                                                  >  1:2058668/1‑69 (MQ=255)

              gCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGaaa                                                  >  1:419839/1‑69 (MQ=255)

                cAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAAc                                                 >  1:1554021/1‑68 (MQ=255)

                    cAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCaa                                             >  1:2697117/1‑68 (MQ=255)

                    cAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAAc                                            >  1:1034774/1‑69 (MQ=255)

                    cAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAAc                                            >  1:2304321/1‑69 (MQ=255)

                    cAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAAc                                            >  1:140993/1‑69 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAgttgt                                                               >  1:1319289/1‑49 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTc                                               >  1:2167078/1‑65 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTc                                               >  1:104841/1‑65 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTc                                               >  1:3241680/1‑65 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTc                                               >  1:1181991/1‑65 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTc                                               >  1:3125390/1‑65 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTc                                               >  1:3017374/1‑65 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTc                                               >  1:3386967/1‑65 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTc                                               >  1:2721222/1‑65 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTc                                               >  1:1287713/1‑65 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTc                                               >  1:259028/1‑65 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTc                                               >  1:2195258/1‑65 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTc                                               >  1:1348254/1‑65 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTc                                               >  1:879877/1‑65 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAAcg                                           >  1:2101194/1‑69 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAAcg                                           >  1:2526913/1‑69 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAAcg                                           >  1:882272/1‑69 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAAcg                                           >  1:2899403/1‑69 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAAcg                                           >  1:2453141/1‑69 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAAcg                                           >  1:1297405/1‑69 (MQ=255)

                     aGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAAcg                                           >  1:907491/1‑69 (MQ=255)

                             aaTCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAacac                                   >  1:1291685/1‑69 (MQ=255)

                             aaTCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAaca                                    >  1:13955/1‑68 (MQ=255)

                             aaTCAATTCATCGGTTTTACTACGGAACTTCGCAGAGTTGTTATCGATGCGAAACTCAACGCGCAaca                                    >  1:1287376/1‑68 (MQ=255)

                              aTCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAacac                                   >  1:3375261/1‑68 (MQ=255)

                                 aaTTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTa                               >  1:1194603/1‑69 (MQ=255)

                                 aaTTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTa                               >  1:1618056/1‑69 (MQ=255)

                                 aaTTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTa                               >  1:201163/1‑69 (MQ=255)

                                 aaTTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTa                               >  1:1561375/1‑69 (MQ=255)

                                 aaTTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTa                               >  1:2844826/1‑69 (MQ=255)

                                 aaTTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTa                               >  1:871409/1‑69 (MQ=255)

                                   ttCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTa                               >  1:973482/1‑67 (MQ=255)

                                   ttCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTAcg                             >  1:1002237/1‑69 (MQ=255)

                                   ttCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTAcg                             >  1:2742160/1‑69 (MQ=255)

                                    tCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTAcgc                            >  1:793029/1‑69 (MQ=255)

                                    tCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTAcgc                            >  1:2421786/1‑69 (MQ=255)

                                         ggTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCa                                       >  1:914864/1‑53 (MQ=255)

                                         ggTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAg                       >  1:2332662/1‑69 (MQ=255)

                                         ggTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAg                       >  1:2992950/1‑69 (MQ=255)

                                             ttACTACGGAACTTCGCTGAATTGTTATCGATGCGAAACTCAACGCGCAGAACTTTa                               >  1:758657/1‑57 (MQ=255)

                                             ttACTACGGAACTTCGCTGAATTGTTATCGATGCGAAACTCAACGCGCAGAACTTTAcg                             >  1:13918/1‑59 (MQ=255)

                                             ttACTACGGAACTTCGCTGAATTGTTATCGATGCGAAACTCAACGCGCAGAACTTTAc                              >  1:2570486/1‑58 (MQ=255)

                                                     gAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAgg                 >  1:573019/1‑63 (MQ=255)

                                                     gAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTa           >  1:2125854/1‑69 (MQ=255)

                                                        cTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGAc              >  1:1358206/1‑63 (MQ=255)

                                                        cTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTaatg        >  1:2983082/1‑69 (MQ=255)

                                                        cTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTaatg        >  1:1654422/1‑69 (MQ=255)

                                                        cTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTaatg        <  1:17457/69‑1 (MQ=255)

                                                        cTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTaat         >  1:1073868/1‑68 (MQ=255)

                                                              tGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGAt  <  1:823138/69‑1 (MQ=255)

                                                              tGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGAt  <  1:797627/69‑1 (MQ=255)

                                                              tGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGAt  <  1:214529/69‑1 (MQ=255)

                                                              tGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGAt  <  1:1417353/69‑1 (MQ=255)

                                                              tGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGAt  <  1:756331/69‑1 (MQ=255)

                                                              tGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGAt  <  1:2053668/69‑1 (MQ=255)

                                                              tGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGAt  <  1:1108012/69‑1 (MQ=255)

                                                              tGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGAt  <  1:2651421/69‑1 (MQ=255)

                                                              tGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGAt  <  1:59085/69‑1 (MQ=255)

CCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAATGAT  >  minE/1904960‑1905090

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