BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	796,336	(T)_8→7	intergenic (‑7/‑419)	ECB_00737 ← / → ECB_00738	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	796,329	0	T	.	96.0%	591.7 / 15.9	126	coding (1/993 nt)	ECB_00737	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/3); new base . (60/61); total (62/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

GGGGACTCTCATGTTAATGTGAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCATTTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGTTATTAAAA  >  REL606/796255‑796399

ggggACTCTCATGTTAATGTGAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCAtt                                                                       >  2:1297806/1‑76 (MQ=255)

  ggACTCTCATGTTAATGTGAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCGtttt                                                                     >  2:1326441/1‑76 (MQ=255)

  ggACTCTCATGTTAATGTGAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCAtttt                                                                     >  1:1035314/1‑76 (MQ=255)

  ggACTCTCATGTTAATGTGAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCAtttt                                                                     >  2:1080808/1‑76 (MQ=255)

  ggACTCTCATGTTAATGTGAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCAtttt                                                                     >  2:3377846/1‑76 (MQ=255)

      tctcATGTTAATGTGAAAAACAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCAtttttttt                                                                 >  1:1230066/1‑76 (MQ=255)

          aTGTTAATGTGAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTtgt                                                            <  1:2649365/76‑1 (MQ=255)

           tGTTAATGTGAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTtgtt                                                           >  1:1681030/1‑76 (MQ=255)

           tGTTAATGTGAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTtgtt                                                           >  1:1840221/1‑76 (MQ=255)

            gTTAATGTGAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTtgtta                                                          >  2:2516179/1‑76 (MQ=255)

             ttAATGTGAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTtgttat                                                         >  2:2535505/1‑76 (MQ=255)

             ttAATGTGAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTtgttat                                                         <  2:1365798/76‑1 (MQ=255)

                aTGTGAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTtgttatgtt                                                      <  2:2180681/76‑1 (MQ=255)

                aTGTGAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTtgttatgtt                                                      >  2:732962/1‑76 (MQ=255)

                 tgtgAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTtgttatgtta                                                     <  1:4642498/76‑1 (MQ=255)

                    gAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTg                                                  >  2:1907382/1‑76 (MQ=255)

                    gAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTg                                                  <  1:4123563/76‑1 (MQ=255)

                     aaaaaTAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTgg                                                 >  2:65605/1‑76 (MQ=255)

                      aaaaTAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGt                                                <  2:298771/76‑1 (MQ=255)

                      aaaaTAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGt                                                <  1:2426801/76‑1 (MQ=255)

                         aTAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCATTTTTTTTCTTGTTATGTTATTGGTAc                                              >  2:4046298/1‑76 (MQ=255)

                         aTAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACt                                             >  2:300029/1‑76 (MQ=255)

                         aTAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACt                                             >  1:1346710/1‑76 (MQ=255)

                          tAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTa                                            >  1:1085934/1‑76 (MQ=255)

                           aaaCACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAg                                           >  1:3241454/1‑76 (MQ=255)

                            aaCACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGt                                          <  2:1886677/76‑1 (MQ=255)

                             acacATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTc                                         <  1:3506346/76‑1 (MQ=255)

                              cacaTAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCa                                        <  1:2194451/76‑1 (MQ=255)

                              cacaTAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCa                                        >  2:2411867/1‑76 (MQ=255)

                                caTAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCAtt                                      >  1:602871/1‑76 (MQ=255)

                                caTAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCAtt                                      <  2:4015108/76‑1 (MQ=255)

                                  tAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTg                                    <  2:3988984/76‑1 (MQ=255)

                                  tAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTg                                    >  1:2839578/1‑76 (MQ=255)

                                  tAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTg                                    <  1:800929/76‑1 (MQ=255)

                                  aaaaaaGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTg                                    <  2:2115059/75‑1 (MQ=255)

                                    aaaaGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAg                                  >  1:4481507/1‑76 (MQ=255)

                                     aaaGCAGCTATAAATCTTAAGATTGTTACTCCATCCATTTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAg                                  >  2:3836889/1‑76 (MQ=255)

                                     aaaGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAgg                                 >  2:4359027/1‑76 (MQ=255)

                                     aaaGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAgg                                 >  2:4466980/1‑76 (MQ=255)

                                     aaaGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAgg                                 <  1:2865717/76‑1 (MQ=255)

                                     aaaGCAGCTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAgg                                 <  1:3178546/76‑1 (MQ=255)

                                       agcagcTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCt                               >  1:2078701/1‑76 (MQ=255)

                                       agcagcTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCt                               >  1:3814666/1‑76 (MQ=255)

                                       agcagcTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCt                               <  1:3060110/76‑1 (MQ=255)

                                        gcagcTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCtt                              <  2:3891649/76‑1 (MQ=255)

                                        gcagcTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCtt                              <  2:4233784/76‑1 (MQ=255)

                                         cagcTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTtc                             >  2:4424787/1‑76 (MQ=255)

                                          agcTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTtct                            >  1:816378/1‑76 (MQ=255)

                                           gcTATAAATCTTAAGATTGTTACTCCATCCATTTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTtct                            <  1:3671572/76‑1 (MQ=255)

                                            cTATAAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCa                          >  1:2295847/1‑76 (MQ=255)

                                             tataAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCAt                         <  1:1508470/76‑1 (MQ=255)

                                             tataAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCAt                         >  2:2741721/1‑76 (MQ=255)

                                             tataAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCAt                         <  2:2476690/76‑1 (MQ=255)

                                              ataAATCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATg                        >  1:298884/1‑76 (MQ=255)

                                                aaaTCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGAt                      >  2:2853448/1‑76 (MQ=255)

                                                 aaTCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGAtt                     <  2:856834/76‑1 (MQ=255)

                                                 aaTCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGAtt                     >  1:2030862/1‑76 (MQ=255)

                                                  aTCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTa                    <  1:185201/76‑1 (MQ=255)

                                                   tCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTaa                   >  2:3447269/1‑76 (MQ=255)

                                                   tCTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTaa                   <  1:79612/76‑1 (MQ=255)

                                                    cTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTaaa                  <  1:1192428/76‑1 (MQ=255)

                                                    cTTAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTaaa                  <  1:2180499/76‑1 (MQ=255)

                                                     ttAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTaaaa                 >  2:151292/1‑76 (MQ=255)

                                                     ttAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTaaaa                 <  2:669616/76‑1 (MQ=255)

                                                     ttAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTaaaa                 >  1:599172/1‑76 (MQ=255)

                                                      tAAGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAg                <  2:3985759/76‑1 (MQ=255)

                                                       aaGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAgg               >  1:3009210/1‑76 (MQ=255)

                                                       aaGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAgg               >  1:4825193/1‑76 (MQ=255)

                                                       aaGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAgg               <  2:582343/76‑1 (MQ=255)

                                                       aaGATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAgg               <  2:2034470/76‑1 (MQ=255)

                                                         gATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGtt             >  1:1443781/1‑76 (MQ=255)

                                                         gATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGtt             <  2:107773/76‑1 (MQ=255)

                                                         gATTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGtt             <  1:5061721/76‑1 (MQ=255)

                                                          aTTGTTCCTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTc            <  2:569413/76‑1 (MQ=255)

                                                          aTTGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTc            <  1:3035556/76‑1 (MQ=255)

                                                           ttGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCg           <  2:2078701/76‑1 (MQ=255)

                                                            tGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGt          >  2:4050479/1‑76 (MQ=255)

                                                            tGTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGt          >  2:1456731/1‑76 (MQ=255)

                                                             gTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGtt         <  2:2727307/76‑1 (MQ=255)

                                                             gTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGtt         <  1:2027825/76‑1 (MQ=255)

                                                             gTTACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGtt         <  2:1463971/76‑1 (MQ=255)

                                                              ttACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGtta        >  2:1434119/1‑76 (MQ=255)

                                                               tACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGttat       >  1:4967674/1‑76 (MQ=255)

                                                               tACTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGttat       >  2:4679583/1‑76 (MQ=255)

                                                                aCTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGttatt      <  1:958891/76‑1 (MQ=255)

                                                                aCTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGttatt      <  2:2976745/76‑1 (MQ=255)

                                                                 cTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGttatta     <  1:3622668/76‑1 (MQ=255)

                                                                 cTCCATCCA‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGttatta     >  1:3238989/1‑76 (MQ=255)

                                                                  tccatccaTTTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGttatta     <  2:3863734/76‑1 (MQ=255)

                                                                  tccatcca‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGTTATTaa    <  1:626398/76‑1 (MQ=255)

                                                                  tccatcca‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGTTATTaa    <  1:65605/76‑1 (MQ=255)

                                                                  tccatcca‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGTTATTaa    >  2:545496/1‑76 (MQ=255)

                                                                  tccatcca‑TTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGTTATTaa    >  2:2470679/1‑76 (MQ=255)

                                                                   ccatccatTTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGTTATTaa    <  2:3427307/76‑1 (MQ=255)

                                                                   ccatcca‑tTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGTTATTaaa   >  2:2061425/1‑76 (MQ=255)

                                                                   ccatcca‑tTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGTTATTaaa   <  1:342183/76‑1 (MQ=255)

                                                                   ccatcca‑tTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGTTATTaaa   <  2:3455334/76‑1 (MQ=255)

                                                                    catcca‑tTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGTTATTaaaa  >  1:4792584/1‑76 (MQ=255)

                                                                    catcca‑tTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGTTATTaaaa  <  2:2675422/76‑1 (MQ=255)

                                                                    catcca‑tTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGTTATTaaaa  <  1:1655206/76‑1 (MQ=255)

GGGGACTCTCATGTTAATGTGAAAAATAAACACATAAAAAGCAGCTATAAATCTTAAGATTGTTACTCCATCCATTTTTTTTCTTGTTATGTTATTGGTACTAGTCATTTGAGGCTTCTCATGATTAAAAGGTTCGTTATTAAAA  >  REL606/796255‑796399

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