BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,021,983	(A)_9→8	coding (702/1140 nt)	ECB_02821 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,021,975	0	A	.	90.0%	155.1 / 16.3	41	coding (710/1140 nt)	ECB_02821	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/3); new base . (21/15); total (23/18)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTCCAAAAAACAAAATAATAATATTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTCAAAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGCATCTCTATGTCCTAAATTA  >  REL606/3021901‑3022037

cTTTCCAAAAAACAAAATAATAATATTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTCaa                                                               >  1:1930172/1‑76 (MQ=255)

 tttCCAAAAAACAAAATAATAATATTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTCaaa                                                              <  1:1841064/76‑1 (MQ=255)

    ccAAAAAACAAAATAATAATATTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTCaaaaaa                                                           >  1:1881017/1‑76 (MQ=255)

      aaaaaaCAAAATAATAATATTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTCaaaaaaaa                                                         >  2:660746/1‑76 (MQ=255)

       aaaaaCAAAATAATAATATTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTCaaaaaaaat                                                        >  2:448956/1‑75 (MQ=255)

         aaaCAAAATAATAATATTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTCaaaaaaaatca                                                      <  1:967181/76‑4 (MQ=255)

         aaaCAAAATAATAATATTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTCaaaaaaaatca                                                      <  2:569885/76‑4 (MQ=255)

           aCAAAATAATAATATTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAAt                                                   <  2:355921/76‑1 (MQ=255)

              aaaTAATAATATTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGca                                                >  2:506127/1‑76 (MQ=255)

                  aataataTTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTCAAAAAAAAATCAATGcattc                                             <  2:43102/76‑1 (MQ=255)

                  aataataTTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCGATGcattca                                            <  1:1996230/76‑1 (MQ=255)

                  aataataTTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGcattca                                            <  1:787382/76‑1 (MQ=255)

                    taataTTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGcattcatt                                          >  1:275306/1‑76 (MQ=255)

                    taataTTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGcattcatt                                          >  2:266468/1‑76 (MQ=255)

                    taataTTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATC‑CCAAAAAAAATCAATGcattcatt                                          >  2:159598/1‑76 (MQ=255)

                     aataTTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCattcattc                                         <  1:544298/76‑1 (MQ=255)

                      atatTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTCAAAAAAAAATCAATGCattcattc                                         <  1:133492/76‑1 (MQ=255)

                      atatTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCAttcattca                                        >  2:210665/1‑76 (MQ=255)

                      atatTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCAttcattca                                        >  2:1498894/1‑76 (MQ=255)

                      atatTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCAttcattca                                        >  2:1358861/1‑76 (MQ=255)

                        atTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCaaa                                      >  1:687002/1‑76 (MQ=255)

                         ttGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAaat                                     <  2:1976473/76‑1 (MQ=255)

                         ttGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAaat                                     <  2:1522436/76‑1 (MQ=255)

                          tGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTCAAAAAAAAATCAATGCATTCATTCAaat                                     >  2:1624854/1‑76 (MQ=255)

                          tGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAaata                                    >  2:654442/1‑76 (MQ=255)

                           gTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAaataa                                   >  1:74593/1‑76 (MQ=255)

                            tCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAaataat                                  <  1:1579264/76‑1 (MQ=255)

                               tACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATc                               >  2:427904/1‑76 (MQ=255)

                                aCAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCg                              >  2:1064536/1‑76 (MQ=255)

                                 cAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGa                             <  1:1717199/76‑1 (MQ=255)

                                     tccatccTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGATACt                         <  2:714736/76‑1 (MQ=255)

                                        atccTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTg                      <  2:1301062/76‑1 (MQ=255)

                                         tccTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGc                     >  1:1466075/1‑76 (MQ=255)

                                          ccTCTACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGCa                    >  2:475901/1‑76 (MQ=255)

                                             ctACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGCAtct                 >  2:1838111/1‑76 (MQ=255)

                                              tACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGCAtctc                >  1:290967/1‑76 (MQ=255)

                                              tACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGCAtctc                >  1:16331/1‑76 (MQ=255)

                                              tACCAAGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGCAtctc                >  1:1305244/1‑76 (MQ=255)

                                                  aaGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGCATCTCTATg            >  1:1246902/1‑76 (MQ=255)

                                                  aaGAAAACAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGCATCTCTATg            >  2:1341775/1‑76 (MQ=255)

                                                       aaCAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGCATCTCTATGTCCTa       <  2:1191686/76‑1 (MQ=255)

                                                       aaCAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGCATCTCTATGTCCTa       >  1:474617/1‑76 (MQ=255)

                                                       aaCAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGCATCTCTATGTCCTa       >  1:1051061/1‑76 (MQ=255)

                                                        aCAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGCATCTCTATGTCCTaa      <  2:98063/76‑1 (MQ=255)

                                                         cAGTTAAATTTGATCTCAAAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGCATCTCTATGTCCTaa      <  2:214703/76‑1 (MQ=255)

                                                         cAGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGCATCTCTATGTCCTaaa     <  1:1301763/76‑1 (MQ=255)

                                                          aGTTAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGCATCTCTATGTCCTAAAt    <  1:1766592/76‑1 (MQ=255)

                                                            ttAAATTTGATCTC‑AAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGCATCTCTATGTCCTAAATTa  <  2:815497/76‑1 (MQ=255)

CTTTCCAAAAAACAAAATAATAATATTGTCGTACAAATCCATCCTCTACCAAGAAAACAGTTAAATTTGATCTCAAAAAAAAATCAATGCATTCATTCAAATAATATCGATACTTTGCATCTCTATGTCCTAAATTA  >  REL606/3021901‑3022037

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