BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,906,894	G→A	5.1%	L477L (CTG→TTG)	pyrG ←	CTP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,906,894	0	G	A	5.1%	516.2 / 10.7	156	L477L (CTG→TTG)	pyrG	CTP synthetase
Reads supporting (aligned to +/- strand): ref base G (80/68); new base A (5/3); total (85/71)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

Only 100 of 161 total aligned reads displayed.

GACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAAC  >  W3110S.gb/2906829‑2906963

     aCTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATgttgtt                                                              >  1:5780813/1‑70 (MQ=255)

     aCTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGttgttg                                                             <  1:2615633/71‑1 (MQ=255)

      cTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGa                                                            <  1:666842/71‑1 (MQ=255)

      cTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAACATGCTGTTGa                                                            <  1:78140/71‑1 (MQ=255)

        gATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACt                                                          >  1:2472298/1‑71 (MQ=255)

        gATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACt                                                          >  1:2430654/1‑71 (MQ=255)

         atcatcCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACtt                                                         >  1:7583018/1‑71 (MQ=255)

         atcatcCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACtt                                                         >  1:3250202/1‑71 (MQ=255)

         atcatcCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACt                                                          <  1:755204/70‑1 (MQ=255)

          tcatcCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTc                                                        >  1:5596907/1‑71 (MQ=255)

            atcCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGt                                                      <  1:1905740/71‑1 (MQ=255)

              ccccGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTaa                                                    >  1:1543935/1‑71 (MQ=255)

              ccccGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTaa                                                    >  1:5927107/1‑71 (MQ=255)

                 cGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATgttgtt                                                              <  1:5632715/58‑1 (MQ=255)

                 cGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATgttgtt                                                              <  1:332618/58‑1 (MQ=255)

                 cGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATgttgtt                                                              <  1:88711/58‑1 (MQ=255)

                 cGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACgg                                                 <  1:3420264/71‑1 (MQ=255)

                    aaCGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACgg                                                 >  1:5622246/1‑68 (MQ=255)

                    aaCGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGa                                              >  1:4764422/1‑71 (MQ=255)

                    aaCGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGa                                              <  1:1207151/71‑1 (MQ=255)

                       ggCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTaa                                                    >  1:4507118/1‑62 (MQ=255)

                        gCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGAt                                          <  1:5271328/71‑1 (MQ=255)

                           cGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGAc                                       >  1:6652018/1‑71 (MQ=255)

                            gCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtgttg                                                                >  1:6045819/1‑45 (MQ=255)

                            gCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTa                                                     <  1:5560745/56‑1 (MQ=255)

                             cAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtg                                                                   >  1:2842959/1‑41 (MQ=255)

                             cAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCAtg                                                                   >  1:5551971/1‑41 (MQ=255)

                             cAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAACAtgct                                                                 <  1:2086630/43‑3 (MQ=39)

                              aaCGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCt                                    <  1:4526120/71‑1 (MQ=255)

                               aCGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACg                                      >  1:4103280/1‑68 (MQ=255)

                               aCGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACg                                      >  1:1527242/1‑68 (MQ=255)

                               aCGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACg                                      >  1:7022088/1‑68 (MQ=255)

                               aCGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACg                                      >  1:7697690/1‑68 (MQ=255)

                               aCGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTc                                   <  1:3442668/71‑1 (MQ=255)

                                cgcgCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCTTAACGGTGACGATGACGCTca                                  <  1:2823407/71‑1 (MQ=255)

                                cgcgCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTca                                  <  1:5355813/71‑1 (MQ=255)

                                  cgcAGACCTGCATCTTCAATCTGTTTCAACAGCATgttgtt                                                              <  1:7097214/41‑1 (MQ=255)

                                  cgcAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTcaa                                 >  1:2424399/1‑70 (MQ=255)

                                   gcAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTcaaca                               <  1:2738138/71‑1 (MQ=255)

                                    cAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTcaacaa                              <  1:4895008/71‑1 (MQ=255)

                                    cAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTcaacaa                              <  1:7534325/71‑1 (MQ=255)

                                    cAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTcaacaa                              >  1:2726579/1‑71 (MQ=255)

                                     aGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAAt                             >  1:3576466/1‑71 (MQ=255)

                                     aGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAAt                             >  1:2507358/1‑71 (MQ=255)

                                       aCCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTg                           >  1:2950971/1‑71 (MQ=255)

                                       aCCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTg                           <  1:5332639/71‑1 (MQ=255)

                                         cTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTc                         >  1:3050151/1‑71 (MQ=255)

                                           gCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAAtt                            >  1:7348667/1‑66 (MQ=255)

                                            cATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCgg                       <  1:6435342/70‑1 (MQ=255)

                                            cATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCg                        <  1:183456/69‑1 (MQ=255)

                                              tCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGt                                                <  1:437329/43‑1 (MQ=255)

                                               cTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAAt                             <  1:3941320/61‑1 (MQ=255)

                                               cTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAAt                             <  1:4167129/61‑1 (MQ=255)

                                                ttCAATCTGTTTCAACAGTATGTTGTTGACTTCGTa                                                     <  1:102238/36‑1 (MQ=255)

                                                ttCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCAt                  >  1:3910395/1‑71 (MQ=255)

                                                 tCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCg                        >  1:6213926/1‑64 (MQ=255)

                                                  cAATCTGTTTCAACAGCATGTTGTTGACTTCGTaa                                                    >  1:3046522/1‑35 (MQ=255)

                                                  cAATCTGTTTCAACAGCATGTTGTTGACTTCGTaa                                                    >  1:588689/1‑35 (MQ=255)

                                                  cAATCTGTTTCAACAGCATGTTGTTGACTTCGTaa                                                    >  1:3513558/1‑35 (MQ=255)

                                                  cAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGgcgc                    <  1:6518760/67‑1 (MQ=255)

                                                  cAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCa                   <  1:7540822/68‑1 (MQ=255)

                                                  cAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTg                >  1:2961718/1‑71 (MQ=255)

                                                   aaTCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCAttc                <  1:3285729/70‑2 (MQ=255)

                                                   aaTCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTg                >  1:1308341/1‑70 (MQ=255)

                                                    aTCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTa              <  1:1437165/71‑1 (MQ=255)

                                                    aTCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTa              <  1:7565780/71‑1 (MQ=255)

                                                     tCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTc                         <  1:3841762/59‑1 (MQ=255)

                                                     tCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTc                         <  1:4864269/59‑1 (MQ=255)

                                                      cTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTcaa                                 >  1:5927059/1‑50 (MQ=255)

                                                      cTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTcaa                                 >  1:5001340/1‑50 (MQ=255)

                                                      cTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTcaa                                 >  1:6181188/1‑50 (MQ=255)

                                                      cTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTa              >  1:5274504/1‑69 (MQ=255)

                                                       tGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATg                                         <  1:166255/41‑1 (MQ=255)

                                                       tGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATg                                         <  1:2887387/41‑1 (MQ=255)

                                                         tttCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGTCGCATTGTACAGc          <  1:100590/70‑1 (MQ=255)

                                                         tttCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCt         <  1:6062504/71‑1 (MQ=255)

                                                         tttCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCt         <  1:219382/71‑1 (MQ=255)

                                                          ttCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTg        >  1:6431195/1‑71 (MQ=255)

                                                          ttCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTg        >  1:4878247/1‑71 (MQ=255)

                                                           tCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGt                          >  1:1628610/1‑52 (MQ=255)

                                                           tCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCg                        >  1:7759790/1‑54 (MQ=255)

                                                           tCAACAACATGCTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTgg       >  1:138587/1‑71 (MQ=255)

                                                           tCAACAACATGCTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTgg       >  1:6267084/1‑71 (MQ=255)

                                                           tCAACAACATGCTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTgg       >  1:7376325/1‑71 (MQ=255)

                                                           tCAACAACATGCTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTgg       >  1:5456301/1‑71 (MQ=255)

                                                           tCAACAACATGCTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTg        >  1:4780415/1‑70 (MQ=255)

                                                            cAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTgg       <  1:7153450/70‑1 (MQ=255)

                                                            cAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTgg       <  1:5043070/70‑1 (MQ=255)

                                                            cAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTgg       <  1:6261595/70‑1 (MQ=255)

                                                            cAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTgg       >  1:6212924/1‑70 (MQ=255)

                                                            cAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTgg       >  1:5672405/1‑70 (MQ=255)

                                                             aacaacATGCTGTTGACTTCGTAACGGTGACGATGACGCt                                    >  1:2471326/1‑40 (MQ=37)

                                                             aacaacATGCTGTTGACTTCGTAACGGTGACGATGACGCt                                    >  1:5200085/1‑40 (MQ=37)

                                                             aaCAGCATGTTGTTGACTTCGTAACGGTGACGATGa                                        >  1:6796493/1‑36 (MQ=255)

                                                             aaCAGCATGTTGTTGACTTCGTAACGGTGACGATGa                                        >  1:2289187/1‑36 (MQ=255)

                                                             aaCAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTgg       <  1:206018/69‑1 (MQ=255)

                                                             aaCAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGc      >  1:6413560/1‑70 (MQ=255)

                                                              aCAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGa    <  1:2160110/71‑1 (MQ=255)

                                                              aCAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGa    <  1:3673381/71‑1 (MQ=255)

                                                               cAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGaa   >  1:6805119/1‑71 (MQ=255)

GACCAACTGATCATCCCCGGAACGGCCCGCAACGCGCAGACCTGCATCTTCAATCTGTTTCAACAGCATGTTGTTGACTTCGTAACGGTGACGATGACGCTCAACAATTGTCGGCGCATTGTACAGCTGGCGAAC  >  W3110S.gb/2906829‑2906963

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