BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	2,114,852	(T)_8→7	coding (187/366 nt)	ECB_02002 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	2,114,845	0	T	.	90.2%	156.1 / 14.1	41	coding (194/366 nt)	ECB_02002	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/3); new base . (17/20); total (18/23)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

AAAATATCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTATTTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACTTTGGAATT  >  REL606/2114772‑2114912

aaaaTATCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAttt                                                                   <  2:1878232/76‑1 (MQ=255)

aaaaTATCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAttt                                                                   <  1:430614/76‑1 (MQ=255)

aaaaTATCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAttt                                                                   <  2:448244/76‑1 (MQ=255)

 aaaTATCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttt                                                                  >  2:1255356/1‑76 (MQ=255)

   atatCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttttt                                                                <  2:474953/76‑1 (MQ=255)

    tatCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAttttttt                                                               <  2:1088315/76‑1 (MQ=255)

     atCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttttttc                                                              >  2:651302/1‑75 (MQ=255)

     atCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttttttc                                                              >  2:651310/1‑75 (MQ=255)

     atCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttttttc                                                              >  2:558538/1‑75 (MQ=255)

     atCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttttttc                                                              <  1:1780647/76‑2 (MQ=255)

      tCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttttttct                                                             >  2:741362/1‑74 (MQ=255)

      tCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttttttct                                                             >  2:1832014/1‑74 (MQ=255)

      tCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttttttct                                                             >  1:106692/1‑74 (MQ=255)

       ccATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttttttctc                                                            <  2:998293/76‑4 (MQ=255)

       ccATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttttttctc                                                            >  2:182725/1‑73 (MQ=255)

       ccATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTAtttttttctc                                                            >  1:1478960/1‑73 (MQ=255)

           cAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCt                                                       <  1:770005/76‑1 (MQ=255)

           cAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCt                                                       <  1:770019/76‑1 (MQ=255)

               cAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTgaga                                                   >  1:1389743/1‑76 (MQ=255)

                aaTGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGagag                                                  <  2:1548964/76‑1 (MQ=255)

                  tGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCa                                                >  2:381310/1‑76 (MQ=255)

                    aGCTAGCAAATTCCTTACCTTTTTCAGATGAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAAt                                              <  1:182163/76‑1 (MQ=255)

                     gCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATc                                             <  1:383188/76‑1 (MQ=255)

                      cTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATcc                                            >  2:674649/1‑76 (MQ=255)

                        aGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATcccc                                          >  1:1694356/1‑76 (MQ=255)

                         gCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCg                                         >  1:1110888/1‑76 (MQ=255)

                          cAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGc                                        <  2:1687943/76‑1 (MQ=255)

                           aaaTTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCa                                       <  1:1734316/76‑1 (MQ=255)

                           aaaTTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCa                                       >  1:693755/1‑76 (MQ=255)

                              ttCCTTACCTGTTTCAGATAAAAATAATTTAACGTTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATg                                    <  1:127644/76‑1 (MQ=255)

                              ttCCTTACCTGTTTCAGATAAAAATAATTTAACATTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATg                                    <  1:127620/76‑1 (MQ=255)

                               tcctttccttTTTCAGATAAAAATAATTTAACACTCTCATTATTTTTTTTCTCAGCTGAGAGCAATCCCCGCAATg                                    <  2:982026/76‑1 (MQ=255)

                                 cTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTATTTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGAc                                  <  1:1456483/76‑1 (MQ=255)

                                   tACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATa                               >  2:1034144/1‑76 (MQ=255)

                                   tACCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATa                               <  1:1407164/76‑1 (MQ=255)

                                    aCCTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATaa                              <  2:1736311/76‑1 (MQ=255)

                                     ccTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAAt                             <  2:529458/76‑1 (MQ=255)

                                      cTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAAtt                            >  1:896447/1‑76 (MQ=255)

                                      cTTTTTCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAAtt                            <  2:551099/76‑1 (MQ=255)

                                         tttCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAAttttt                         >  2:1444787/1‑76 (MQ=255)

                                           tCAGATAAAAATAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTtat                       <  1:388867/76‑1 (MQ=255)

                                             aGATAAAAATAATTTAACACTCTCATTATTTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTtata                      <  2:45416/76‑1 (MQ=255)

                                                  aaaaTAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGatat                <  1:795857/76‑1 (MQ=255)

                                                   aaaTAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATAtt               >  2:744447/1‑76 (MQ=255)

                                                   aaaTAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATAtt               >  2:81580/1‑76 (MQ=255)

                                                   aaaTAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATAtt               >  2:999238/1‑76 (MQ=255)

                                                   aaaTAATTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATAtt               <  2:120344/76‑1 (MQ=255)

                                                    aataatTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATAttt              <  1:720331/76‑1 (MQ=255)

                                                       aatTTAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCAc           <  1:806616/76‑1 (MQ=255)

                                                         tttAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACtt         >  2:1543670/1‑76 (MQ=255)

                                                         tttAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACtt         >  2:1422389/1‑76 (MQ=255)

                                                          ttAACACTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACttt        >  2:1637078/1‑76 (MQ=255)

                                                            aaCACTCTCATTATTTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACTTTg       >  2:1069770/1‑76 (MQ=255)

                                                              cacTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTACACTTTGGaa    >  2:436423/1‑76 (MQ=255)

                                                               acTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACTTTGGAAt   >  2:1621919/1‑76 (MQ=255)

                                                               acTCTCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACTTTGGAAt   <  1:364796/76‑1 (MQ=255)

                                                                ctctCATTA‑TTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACTTTGGAAtt  <  2:1871908/76‑1 (MQ=255)

AAAATATCCATCAGTCAATGAGCTAGCAAATTCCTTACCTTTTTCAGATAAAAATAATTTAACACTCTCATTATTTTTTTTCTCAGCTGAGAGCAATCCCCGCAATGACATAATTTTTATAAGATATTTCACTTTGGAATT  >  REL606/2114772‑2114912

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