BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,455,159	C→T	9.0%	V113I (GTT→ATT)	yfcO ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,455,159	0	C	T	9.0%	199.4 / 8.5	67	V113I (GTT→ATT)	yfcO	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (35/25); new base T (2/4); total (37/30)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

AGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTTGCTTTGTTGAT  >  W3110S.gb/2455096‑2455222

agTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCa                                                             <  1:1386496/68‑1 (MQ=255)

agTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCa                                                             <  1:1995097/68‑1 (MQ=255)

      tcatAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGta                                                       <  1:942932/68‑1 (MQ=255)

      tcatAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGta                                                       <  1:2218997/68‑1 (MQ=255)

      tcatAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGta                                                       <  1:3673653/68‑1 (MQ=255)

       catAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATg                                                         <  1:2809457/65‑1 (MQ=255)

       catAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATg                                                         <  1:3126424/65‑1 (MQ=255)

       catAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATg                                                         <  1:525843/65‑1 (MQ=255)

       catAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGtat                                                      >  1:105353/1‑68 (MQ=255)

       catAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGtat                                                      >  1:313181/1‑68 (MQ=255)

       catAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGtat                                                      >  1:2399219/1‑68 (MQ=255)

       catAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGtat                                                      >  1:1874644/1‑68 (MQ=255)

       catAAAAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTc                                                              >  1:2202754/1‑60 (MQ=255)

        ttaataaCCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGtat                                                      <  1:192020/66‑1 (MQ=255)

           ataaCCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATAc                                                 >  1:2877301/1‑69 (MQ=255)

            taaCCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGtat                                                      <  1:2109383/63‑1 (MQ=255)

              accaaccaTGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATAcc                                                >  1:391994/1‑67 (MQ=255)

                caaccaTGGACAAGAAAAGTAAGATTCAATCTCTTCAGAGACAATAATCTCAAATGTATAATACCCGc                                             <  1:3314044/68‑1 (MQ=255)

                caaccaTGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAATCTCAAATGTATAATACCCGc                                             <  1:59858/68‑1 (MQ=255)

                caaccaTGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAATCTCAAATGTATAATACCCGc                                             >  1:2986114/1‑68 (MQ=255)

                caaccaTGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAATCTCAAATGTATAATACCCGc                                             <  1:1208442/68‑1 (MQ=255)

                caaccaTGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAATCTCAAATGTATAATACCCGc                                             <  1:2988592/68‑1 (MQ=255)

                caaccaTGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAATCTCAAATGTATAATACCCGc                                             >  1:3728548/1‑68 (MQ=255)

                caaccaTGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCa                                            >  1:579103/1‑69 (MQ=255)

                      tcgACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGtat                                                      <  1:1501027/51‑1 (MQ=255)

                      tGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCa                                      >  1:1337670/1‑69 (MQ=255)

                        gACAAGAAAAGTAAGATTGAATCTCTGCAGAGACAATAACCTCAAATGtat                                                      <  1:3319915/51‑1 (MQ=255)

                        gACAAGAAAAGTAAGATTGAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCa                                      >  1:3297002/1‑67 (MQ=255)

                        gACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCa                                      >  1:1833390/1‑67 (MQ=255)

                        gACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCa                                      >  1:1418467/1‑67 (MQ=255)

                          cAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCa                                  >  1:169963/1‑69 (MQ=255)

                            aGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCa                                  >  1:1663243/1‑67 (MQ=255)

                              aaaaGTAAGATTCAATCTCTTCAGAGACAATAACCTCAAATGt                                                        >  1:2702021/1‑43 (MQ=255)

                              aaaaGTAAGATTCAATCTCTTCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGAc                               >  1:1112874/1‑68 (MQ=255)

                               aaaGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACg                              >  1:2832691/1‑68 (MQ=255)

                                aaGTAAGATTCAATCTCTGCAGAGACAATAAGCTCAAATGta                                                       <  1:966920/42‑1 (MQ=255)

                                aaGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGta                                                       <  1:2739043/42‑1 (MQ=255)

                                aaGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGta                                                       <  1:1594738/42‑1 (MQ=255)

                                aaGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTcc                                       >  1:3195871/1‑58 (MQ=255)

                                aaGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTcc                                       >  1:628915/1‑58 (MQ=255)

                                aaGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTcc                                       >  1:64062/1‑58 (MQ=255)

                                aaGTAAGATTCAATCTCTCCAGAGACAATAACCTCAAATGta                                                       <  1:1783021/42‑1 (MQ=255)

                                 aGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGt                            <  1:29922/68‑1 (MQ=255)

                                  gTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTaa                          <  1:3202469/69‑1 (MQ=255)

                                  gTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTaa                          <  1:1653020/69‑1 (MQ=255)

                                  gTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTaa                          <  1:811345/69‑1 (MQ=255)

                                  gTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTaa                          <  1:3254356/69‑1 (MQ=255)

                                  gTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTaa                          >  1:886612/1‑69 (MQ=255)

                                       aTTCAATCTCTGCAGAGACAATAACCTCAAATGTataata                                                  <  1:2922196/40‑1 (MQ=255)

                                         tCAATCTCTGCAGAGACAATAACCTCAAATGTATAATAc                                                 >  1:685946/1‑39 (MQ=255)

                                         tCAATCTCTGCAGAGACAATAACCTCAAATGTATAATAc                                                 >  1:1138597/1‑39 (MQ=255)

                                         tCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGt                                          >  1:2162971/1‑46 (MQ=255)

                                         tCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAgg                    >  1:296018/1‑68 (MQ=255)

                                         tCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAgg                    >  1:1953440/1‑68 (MQ=255)

                                             tctcTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAAc               >  1:2699794/1‑69 (MQ=255)

                                               tctGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACtt             >  1:2656771/1‑69 (MQ=255)

                                                  gCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATc                        <  1:2105619/55‑1 (MQ=255)

                                                   cAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCGAGACGGTAATCCAGGACAACTTGCtt         >  1:1475833/1‑69 (MQ=255)

                                                   cAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTa                           >  1:3313609/1‑51 (MQ=255)

                                                   cAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTa                           >  1:3181856/1‑51 (MQ=255)

                                                   cAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTTGc           <  1:2725201/67‑1 (MQ=255)

                                                   cAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTTGCtt         >  1:3398828/1‑69 (MQ=255)

                                                      agaCAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTTGCTttgt      <  1:603780/69‑1 (MQ=255)

                                                          aaTAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTTGCTTTGTTGat  >  1:3390319/1‑69 (MQ=255)

                                                          aaTAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTTGCTTTGTTGat  >  1:632900/1‑69 (MQ=255)

                                                          aaTAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTTGCTTTGTTGat  >  1:2729833/1‑69 (MQ=255)

                                                            tAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACa                 <  1:430285/52‑1 (MQ=255)

AGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTTGCTTTGTTGAT  >  W3110S.gb/2455096‑2455222

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