BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,096,973	A→C	F200L (TTT→TTG)	glpF ←	glycerol facilitator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,096,973	0	A	C	98.2%	206.6 / ‑3.7	54	F200L (TTT→TTG)	glpF	glycerol facilitator
Reads supporting (aligned to +/- strand): ref base A (0/1); new base C (24/29); total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

TAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAGCAAGGGAGCCAAAGG  >  REL606/4096900‑4097044

tAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCcaa                                                                       <  1:1673646/76‑1 (MQ=255)

 aaaGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCcaaa                                                                      >  1:1307159/1‑76 (MQ=255)

 aaaGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCcaaa                                                                      >  1:1572527/1‑76 (MQ=255)

   aGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAAcc                                                                    <  1:1318256/76‑1 (MQ=255)

      cGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGt                                                                 <  2:1638152/76‑1 (MQ=255)

       gACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTc                                                                >  2:1030783/1‑76 (MQ=255)

         cATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCaa                                                              <  1:949177/76‑1 (MQ=255)

          aTTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAAt                                                             <  2:1460868/76‑1 (MQ=255)

             gCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATggg                                                          <  2:1081916/76‑1 (MQ=255)

             gCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATggg                                                          <  2:218164/76‑1 (MQ=255)

               cccAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCCAgggggc                                                        >  1:243460/1‑74 (MQ=255)

               cccAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGCGTTCATGGCCAAACCCTTCACTGTGcc                                                        >  1:243456/1‑76 (MQ=255)

                ccAGCCCGCCAGCCAGGCAATGACTTGCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGccc                                                       <  2:1225307/76‑1 (MQ=255)

                 cAGCCCGCCAGCCAGGCAGAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCa                                                      <  2:345312/76‑1 (MQ=255)

                 cAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCa                                                      >  2:1871962/1‑76 (MQ=255)

                  aGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCTCGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCAt                                                     <  1:1835943/76‑1 (MQ=255)

                    cccGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAg                                                   <  1:400481/76‑1 (MQ=255)

                              agGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAAt                                         >  2:50115/1‑76 (MQ=255)

                                gCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGccn                                                       <  1:1046485/60‑2 (MQ=255)

                                gCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGccc                                                       >  2:1046485/1‑60 (MQ=255)

                                 cAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGAc                                      >  2:1635354/1‑76 (MQ=255)

                                  aaaGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGAcc                                     <  2:1864707/76‑1 (MQ=255)

                                  aaaGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGAcc                                     >  2:397986/1‑76 (MQ=255)

                                  aaaGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGAcc                                     >  2:762036/1‑76 (MQ=255)

                                  aaaGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGAcc                                     >  1:1217917/1‑76 (MQ=255)

                                   aaGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCg                                    >  1:828630/1‑76 (MQ=255)

                                      aCTTTCGGACCGAAGTCACGAGCTGGGTTCAAGGCCAAAACCGTCAGTGGATTAATAGATGCGCCAATGATAGCaa                                 >  2:359967/1‑76 (MQ=255)

                                       cTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAAt                                <  2:789263/76‑1 (MQ=255)

                                        tttCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATc                               <  1:981927/76‑1 (MQ=255)

                                        tttCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATc                               <  1:1241391/76‑1 (MQ=255)

                                        tttCGGACCGAAATCGCGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATc                               <  2:1900369/76‑1 (MQ=255)

                                        tttCGGACCCAAGACTCGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATc                               <  1:476876/76‑1 (MQ=255)

                                        tttCGGAACGAAGTTACGCGCTGGGTTCAAGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATc                               <  1:1010604/76‑1 (MQ=255)

                                        tttCGAACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATc                               <  2:1665662/76‑1 (MQ=255)

                                         ttCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCa                              <  1:1702764/76‑1 (MQ=255)

                                          tCGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAg                             >  1:1693659/1‑76 (MQ=255)

                                           cGGACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGt                            >  1:249774/1‑76 (MQ=255)

                                            ggACCGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTa                           >  2:1471699/1‑76 (MQ=255)

                                               ccGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGAc                        <  1:504597/76‑1 (MQ=255)

                                               ccGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGAc                        <  2:1126938/76‑1 (MQ=255)

                                                cGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCTAATGACCGCAATCAGTAGAcc                       <  2:1850442/76‑1 (MQ=255)

                                                cGAAGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGAcc                       >  1:541818/1‑76 (MQ=255)

                                                  aaGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCaa                     <  2:1041565/76‑1 (MQ=255)

                                                  aaGTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGAGGCGCCAATGACCGCAATCAGTAGACCaa                     <  2:1225568/76‑1 (MQ=255)

                                                    gTCACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATc                   <  2:1282722/76‑1 (MQ=255)

                                                      cACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAg                 <  2:248600/76‑1 (MQ=255)

                                                      cACGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAg                 <  2:541818/76‑1 (MQ=255)

                                                       aCGCGCTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAGc                <  1:1741300/76‑1 (MQ=255)

                                                           gcTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAGCAAgg            >  2:957509/1‑76 (MQ=255)

                                                           gcTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAGCAAgg            >  1:385017/1‑76 (MQ=255)

                                                           gcTGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAGCAAgg            <  1:1165399/76‑1 (MQ=255)

                                                             tGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAGCAAGGGa          >  2:1433649/1‑76 (MQ=255)

                                                             tGGGTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAGCAAGGGa          >  2:1204577/1‑76 (MQ=255)

                                                              gggTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAGCAAGGGAg         >  2:1859823/1‑76 (MQ=255)

                                                               ggTTCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAGCAAGGGAGc        >  2:854599/1‑76 (MQ=255)

                                                                 ttCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAGCAAGGGAGCCa      <  2:635507/76‑1 (MQ=255)

                                                                  tCATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAGCAAGGGAGCCaa     <  1:624318/76‑1 (MQ=255)

                                                                   cacggCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAGCAAGGGAGCCaaa    >  2:1385295/4‑76 (MQ=255)

                                                                   cATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAGCAAGGGAGCCaaa    >  2:1858059/1‑76 (MQ=255)

                                                                   cATGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAGCAAGGGAGCCaaa    >  2:645071/1‑76 (MQ=255)

                                                                     tGGCCAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAGCAAGGGAGCCAAAgg  >  2:891422/1‑76 (MQ=255)

TAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAACCCGTCAATGGGCCCATAGATGCGCCAATGACCGCAATCAGTAGACCAATCAGCAAGGGAGCCAAAGG  >  REL606/4096900‑4097044

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 20 ≤` `ATCG/ATCG` `< 28 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`