BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,672,374	G→A	42.1%	G296G (GGC→GGT)	hisS ←	histidyl tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,672,374	0	G	A	42.1%	45.7 / 90.1	121	G296G (GGC→GGT)	hisS	histidyl tRNA synthetase
Reads supporting (aligned to +/- strand): ref base G (27/43); new base A (42/9); total (69/52)
Fisher's exact test for biased strand distribution p-value = 1.23e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 121 total aligned reads displayed.

CTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGT  >  minE/1672311‑1672441

cTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAg                                                                <  1:2255711/69‑1 (MQ=255)

cTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAg                                                                <  1:2139023/69‑1 (MQ=255)

cTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAg                                                                <  1:667023/69‑1 (MQ=255)

cTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAg                                                                <  1:1404452/69‑1 (MQ=255)

cTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAg                                                                <  1:460308/69‑1 (MQ=255)

     cTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGCCCGCCCAg                                                                >  1:684259/1‑64 (MQ=255)

     cTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAgttgt                                                            >  1:1838497/1‑68 (MQ=255)

     cTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAgttgt                                                            >  1:2275247/1‑68 (MQ=255)

      tAACAATACAAGACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTc                                                          <  1:1251622/69‑1 (MQ=255)

      tAACAATACAAGACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTc                                                          <  1:817941/69‑1 (MQ=255)

      tAACAATACAAGACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTc                                                          <  1:864342/69‑1 (MQ=255)

      tAACAATACAAGACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTc                                                          <  1:2311112/69‑1 (MQ=255)

      tAACAATACAAGACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTc                                                          <  1:2935790/69‑1 (MQ=255)

      tAACAATACAAGACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTc                                                          <  1:2707096/69‑1 (MQ=255)

       aaCAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTc                                                          <  1:389621/68‑1 (MQ=255)

       aaCAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTc                                                          >  1:2728569/1‑68 (MQ=255)

       aaCAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTc                                                          <  1:2106995/68‑1 (MQ=255)

       aaCAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTc                                                          <  1:1899059/68‑1 (MQ=255)

       aaCAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTc                                                          >  1:226158/1‑68 (MQ=255)

       aaCAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTc                                                          <  1:1612090/68‑1 (MQ=255)

       aaCAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTc                                                          <  1:1395720/68‑1 (MQ=255)

       aaCAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTc                                                          <  1:1421096/68‑1 (MQ=255)

        aCAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTc                                                          <  1:448684/67‑1 (MQ=255)

           aTACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTT‑CACGACCGCCCAGTTGTTCCACAAg                                                    <  1:891683/69‑1 (MQ=255)

            tACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAg                                                    <  1:834426/69‑1 (MQ=255)

             aCAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAg                                                    <  1:2480970/68‑1 (MQ=255)

             aCAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGa                                                   >  1:3326147/1‑69 (MQ=255)

             aCAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGa                                                   >  1:3297168/1‑69 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:2516845/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:3062105/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:310196/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:101093/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:3151085/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:3210797/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:1040720/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:1090103/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:1091622/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:1165109/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:899628/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:1207152/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:858861/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:690950/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:192743/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:895755/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:2002471/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:2013116/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:425949/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:1266857/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:2212807/1‑67 (MQ=255)

                 gACGTTCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACCACCACCCAGTTGTTCCACCAGACCg                                                >  1:1719740/1‑68 (MQ=255)

                 gACGTTCGAGACCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGAcc                                                 >  1:1547930/1‑67 (MQ=255)

                    gTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAg                                                                <  1:1139214/49‑1 (MQ=255)

                      tggAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGACCGTCATa                                           >  1:2359302/3‑68 (MQ=255)

                      tCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGACCGTCATa                                           >  1:2996130/1‑68 (MQ=255)

                      tCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGACCGTCATa                                           >  1:2833430/1‑68 (MQ=255)

                      tCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGACCGTCATa                                           >  1:942618/1‑68 (MQ=255)

                      tCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGACCGTCATa                                           >  1:1804805/1‑68 (MQ=255)

                      tCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGACCGTCATa                                           >  1:1845543/1‑68 (MQ=255)

                      tCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGACCGTCATa                                           >  1:2701098/1‑68 (MQ=255)

                      tCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGACCGTCATa                                           >  1:2015803/1‑68 (MQ=255)

                      tCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGACCGTCATa                                           >  1:1667968/1‑68 (MQ=255)

                      tCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGACCGTCATa                                           >  1:2346383/1‑68 (MQ=255)

                      tCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGACCGTCATa                                           >  1:1552378/1‑68 (MQ=255)

                      tCGAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGACCGTCATa                                           >  1:2259009/1‑68 (MQ=255)

                      tCGAGGACCATAGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTCCACCAGACCGTCATa                                           >  1:437280/1‑68 (MQ=255)

                       cGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACaa                                                     >  1:259884/1‑57 (MQ=255)

                       cGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCAt                                            >  1:2786561/1‑66 (MQ=255)

                       cGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCAt                                            >  1:609339/1‑66 (MQ=255)

                       cGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCAt                                            >  1:2000975/1‑66 (MQ=255)

                       cGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCAt                                            >  1:2056806/1‑66 (MQ=255)

                       cGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCAt                                            >  1:1150667/1‑66 (MQ=255)

                       cGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAAc                                         >  1:2104048/1‑69 (MQ=255)

                       cGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAAc                                         >  1:782818/1‑69 (MQ=255)

                       cGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAAc                                         >  1:281572/1‑69 (MQ=255)

                       cGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAAc                                         >  1:311717/1‑69 (MQ=255)

                       cGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAAc                                         >  1:3387797/1‑69 (MQ=255)

                        gAGGCCCATCGCAAAACCGACTGCCGGTGTTGCACGACCACCCAGTTGTTc                                                          <  1:2411168/51‑1 (MQ=255)

                            cccATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCg                                    <  1:2637901/69‑1 (MQ=255)

                            cccATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCg                                    <  1:174169/69‑1 (MQ=255)

                            cccATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCg                                    <  1:783147/69‑1 (MQ=255)

                            cccATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCg                                    <  1:1936059/69‑1 (MQ=255)

                                  gCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGca                              >  1:1454035/1‑69 (MQ=255)

                                  gCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGca                              >  1:453685/1‑69 (MQ=255)

                                  gCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGca                              >  1:457896/1‑69 (MQ=255)

                                  gCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGca                              >  1:78831/1‑69 (MQ=255)

                                  gCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGca                              >  1:146236/1‑69 (MQ=255)

                                       aCCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCacacac                         <  1:2333963/69‑1 (MQ=255)

                                                         aCGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCga        <  1:1406717/68‑1 (MQ=255)

                                                         aCGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCga        <  1:2800228/68‑1 (MQ=255)

                                                         aCGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCga        <  1:645832/68‑1 (MQ=255)

                                                         aCGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCga        <  1:578271/68‑1 (MQ=255)

                                                         aCGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCga        <  1:3379105/68‑1 (MQ=255)

                                                         aCGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCga        <  1:2594143/68‑1 (MQ=255)

                                                         aCGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCga        <  1:2883189/68‑1 (MQ=255)

                                                             ccGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTg   <  1:591013/69‑1 (MQ=255)

                                                             ccGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTg   <  1:3306597/69‑1 (MQ=255)

                                                              cGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGt  <  1:202862/69‑1 (MQ=255)

                                                              cGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGt  <  1:84773/69‑1 (MQ=255)

                                                              cGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGt  <  1:3016707/69‑1 (MQ=255)

                                                              cGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGt  <  1:1054772/69‑1 (MQ=255)

CTGTACTAACAATACAAGACGTTCGAGGCCCATAGCAAAACCGACAGCCGGTGTTGCACGACCGCCCAGTTGTTCCACAAGACCGTCATAACGACCGCCTGCACACACGGTGCCCTGGGAGCCGAGACTGT  >  minE/1672311‑1672441

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