BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	BS168_reference_genome	2,271,424	T→C	R78R (AGA→AGG)	uvrX ←	lesion bypass phage DNA polymerase; phage SPbeta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	2,271,424	0	T	C	99.6%	833.1 / NA	252	R78R (AGA→AGG)	uvrX	lesion bypass phage DNA polymerase; phage SPbeta
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (139/112); total (139/113)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 252 total aligned reads displayed.

TTCAGGATCTTCTTTTCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATTCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAACACTACACTTCCCTGTCTATCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCC  >  BS168_reference_genome/2271277‑2271570

ttcAGGATCTTCTTTTCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCt                                                                                                                                                  >  2:2239543/1‑150 (MQ=255)

ttcAGGATCTTCTTTTCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCt                                                                                                                                                  <  1:27990/150‑1 (MQ=255)

 tcAGGATCTTCTTTTCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCtt                                                                                                                                                 <  1:68247/150‑1 (MQ=255)

  cAGGATCTTCTTTTCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTg                                                                                                                                                <  1:1278200/150‑1 (MQ=255)

  cAGGATCTTCTTTTCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTg                                                                                                                                                >  2:536979/1‑150 (MQ=255)

  cAGGACCTCCTCTCCCTGGATCTAAAAACGTTTCAGCAATTGAATACGTATGGACACATTTTTCAGGTTCAAACCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGACTTACAATGTGTATCCTTg                                                                                                                                                <  2:2452537/150‑1 (MQ=11)

   aGGATCTTCTTTTCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTgg                                                                                                                                               <  1:2532403/150‑1 (MQ=255)

    ggATCTTCTTTTCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGa                                                                                                                                              >  2:2908319/1‑150 (MQ=255)

     gATCTTCTTTTCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGAt                                                                                                                                             >  2:3720940/1‑150 (MQ=255)

      aTCTTCTTTTCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATc                                                                                                                                            >  2:928564/1‑150 (MQ=255)

            ttttCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAgg                                                                                                                                      <  1:1141784/150‑1 (MQ=255)

            ttttCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAgg                                                                                                                                      <  2:3194951/150‑1 (MQ=255)

             tttCCTGCATCTAAAATAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTTAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAGAAGCTTCATCTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGt                                                                                                                                     <  1:1438042/150‑1 (MQ=255)

             tttCCTGCATCTAAAAAAGATTCATCAAGTGAATAAGTATGGACTTATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGAGGATTTACAATGTGTATCCTTGGTTCTTCTGGt                                                                                                                                     >  1:1895552/1‑150 (MQ=14)

              ttCCTGCATCTAAAAAAGATTCATCAATTGAATCCGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTa                                                                                                                                    <  1:1694394/150‑1 (MQ=255)

               tCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTAt                                                                                                                                   <  1:2348958/150‑1 (MQ=255)

                ccTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTAtc                                                                                                                                  <  1:105069/150‑1 (MQ=255)

                ccTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTAtc                                                                                                                                  <  2:838073/150‑1 (MQ=255)

                ccTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTAtc                                                                                                                                  >  2:1165708/1‑150 (MQ=255)

                  tGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTAtctc                                                                                                                                >  1:3215143/1‑150 (MQ=255)

                   gCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCa                                                                                                                               >  1:2448777/1‑150 (MQ=255)

                     aTCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCaaa                                                                                                                             >  2:1120561/1‑150 (MQ=255)

                       cTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATa                                                                                                                           <  2:3227389/150‑1 (MQ=255)

                        tAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAg                                                                                                                          >  2:2318766/1‑150 (MQ=255)

                          aaaaaGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGtc                                                                                                                        >  2:301513/1‑150 (MQ=255)

                            aaaGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGtctc                                                                                                                      >  2:2717498/1‑150 (MQ=255)

                            aaaGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGtctc                                                                                                                      >  2:3324418/1‑150 (MQ=255)

                            aaaGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGtctc                                                                                                                      >  2:664177/1‑150 (MQ=255)

                               gATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTAACTCAAATAGTCTCGAt                                                                                                                   <  1:373810/150‑1 (MQ=255)

                                aTTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATc                                                                                                                  <  2:143199/150‑1 (MQ=255)

                                 ttCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATcc                                                                                                                 >  1:2855798/1‑150 (MQ=255)

                                 tccatcaATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATcc                                                                                                                 <  1:460687/148‑1 (MQ=255)

                                   catcaATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATcccg                                                                                                               >  1:265080/1‑148 (MQ=255)

                                   catcaATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTg                                                                                                               >  2:72159/1‑150 (MQ=255)

                                      caATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGttt                                                                                                            <  2:3527342/150‑1 (MQ=255)

                                       aaTTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAGCTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGtttt                                                                                                           >  1:2906524/1‑150 (MQ=255)

                                           gAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGAtt                                                                                                       >  1:1995445/1‑150 (MQ=255)

                                             aTACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTcc                                                                                                     <  1:2125674/150‑1 (MQ=255)

                                              tACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATGCCa                                                                                                    >  1:1915095/1‑150 (MQ=255)

                                                 gTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCaaaa                                                                                                 >  1:1714656/1‑150 (MQ=255)

                                                          caTTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTa                                                                                        <  2:608597/150‑1 (MQ=37)

                                                           aTTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTaa                                                                                       >  2:2091348/1‑150 (MQ=37)

                                                             ttCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGt                                                                                     >  2:3461804/1‑150 (MQ=37)

                                                               ctcAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGc                                                                                   <  1:323988/150‑1 (MQ=25)

                                                               ctcAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGc                                                                                   <  1:477222/150‑1 (MQ=25)

                                                                tcAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCa                                                                                  <  1:1581143/150‑1 (MQ=21)

                                                                 cAGGCACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAg                                                                                 <  1:3761214/150‑1 (MQ=21)

                                                                 cAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAg                                                                                 >  1:2969643/1‑150 (MQ=21)

                                                                  aGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAgg                                                                                >  2:3902131/1‑150 (MQ=21)

                                                                   ggAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTCCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGga                                                                               <  1:829828/150‑1 (MQ=21)

                                                                     aaCAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGgaga                                                                             <  1:899937/150‑1 (MQ=21)

                                                                     aaCAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGgaga                                                                             >  1:2892873/1‑150 (MQ=21)

                                                                      aCAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGAt                                                                            >  1:2241859/1‑150 (MQ=21)

                                                                      aCAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGAt                                                                            <  2:1606158/150‑1 (MQ=21)

                                                                        aaaTCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGGTGc                                                                          >  2:3946582/1‑150 (MQ=21)

                                                                         aaTCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCa                                                                         >  2:182549/1‑150 (MQ=21)

                                                                            cTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCt                                                                      <  1:714052/150‑1 (MQ=21)

                                                                                aaaaCAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTaata                                                                  <  2:415680/150‑3 (MQ=21)

                                                                                aaaaCAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTACCTCACATAGTCTCGATCCTGTTTAGATTCCCAAATCCTTTTTAAGTGCAGGAGATGCTGCTaaga                                                                  >  2:3418014/1‑148 (MQ=11)

                                                                                  aaCAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAatact                                                                <  2:674096/150‑1 (MQ=21)

                                                                                      gCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGGGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGACTCGATCCTGTTTTGATTCCAAAATCCTTTTTAAGTGCAGGAGATGCAGCTAatactacac                                                            >  1:1459117/1‑150 (MQ=11)

                                                                                             aaTTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTg                                                     >  1:3650619/1‑150 (MQ=21)

                                                                                             aaTTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTACCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTg                                                     <  2:1342613/150‑1 (MQ=255)

                                                                                              aTTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGt                                                    >  2:3760020/1‑150 (MQ=21)

                                                                                              aTTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGGGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGt                                                    >  1:3106421/1‑150 (MQ=255)

                                                                                               tttCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTc                                                   >  1:3929448/1‑150 (MQ=18)

                                                                                                ttCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAATTCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCt                                                  >  2:3114132/1‑150 (MQ=12)

                                                                                                ttCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCt                                                  >  1:1814360/1‑150 (MQ=17)

                                                                                                 tCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCtc                                                 >  1:217781/1‑149 (MQ=17)

                                                                                                  cAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCtct                                                >  2:551815/1‑148 (MQ=12)

                                                                                                   aGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGtctctc                                               <  2:2327855/150‑1 (MQ=12)

                                                                                                   aGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAAGAGATGCAGCTAATACTACACTTCCCTGtctctc                                               >  1:3720745/1‑150 (MQ=255)

                                                                                                        aaaCTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTAt                                          <  1:1802658/150‑1 (MQ=12)

                                                                                                          aCTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTc                                        <  1:3274286/150‑1 (MQ=12)

                                                                                                          aCTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTc                                        >  2:1945065/1‑150 (MQ=12)

                                                                                                          aCTCTGATGAAAAGCTTCATTTGTGGATATACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCCATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTAtac                                        <  2:1587077/150‑3 (MQ=11)

                                                                                                           ctctGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTcc                                       >  2:3587085/1‑150 (MQ=12)

                                                                                                           ctctGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTcc                                       >  2:590462/1‑150 (MQ=12)

                                                                                                            tctGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTccc                                      <  2:1616712/150‑1 (MQ=12)

                                                                                                              tgatgaAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGGATTCCCTa                                    >  2:3643346/1‑150 (MQ=255)

                                                                                                                 tgaAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTacaa                                 <  1:3343783/150‑1 (MQ=12)

                                                                                                                 tgaAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTacaa                                 <  1:3886706/150‑1 (MQ=12)

                                                                                                                   aaaaGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCCCTCCGTATTCCCTacaaca                               >  1:3472435/1‑150 (MQ=255)

                                                                                                                      aGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCa                            >  2:1712829/1‑150 (MQ=255)

                                                                                                                         ttCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGa                         >  2:3144107/1‑150 (MQ=255)

                                                                                                                          tCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAAGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGAt                        >  2:745159/1‑150 (MQ=255)

                                                                                                                            aTTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTACTACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAg                      >  2:1503745/1‑150 (MQ=255)

                                                                                                                            aTTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAg                      <  1:1375016/150‑1 (MQ=255)

                                                                                                                            aTTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACGACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGaatg                      >  2:1905377/1‑147 (MQ=11)

                                                                                                                               tgtgGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCAt                   <  2:2843555/150‑1 (MQ=255)

                                                                                                                                gtgGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATg                  <  1:3902131/150‑1 (MQ=255)

                                                                                                                                gtgGATTTACAATGTGTATCCTTGGACCCTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGTAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATg                  <  1:3369700/150‑1 (MQ=11)

                                                                                                                                   gATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTa               >  2:254600/1‑150 (MQ=255)

                                                                                                                                   gATTTACAATGTGTACCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTa               <  2:2741679/150‑1 (MQ=255)

                                                                                                                                     tttACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGGTCCTGTTTTGATTCCAAAAGCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTaaa             >  1:631969/1‑150 (MQ=11)

                                                                                                                                       tACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAgg           <  2:4041495/150‑1 (MQ=255)

                                                                                                                                            tgtgTACCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTaa      <  1:2318766/150‑1 (MQ=255)

                                                                                                                                             gtgtATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAg     <  2:894290/150‑1 (MQ=255)

                                                                                                                                              tgtATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCACCAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGGTAAGc    >  1:3980185/1‑150 (MQ=11)

                                                                                                                                                tATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCTTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGccc  <  1:1545851/150‑1 (MQ=255)

TTCAGGATCTTCTTTTCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATTCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAACACTACACTTCCCTGTCTATCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCC  >  BS168_reference_genome/2271277‑2271570

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