BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	608,214:1	+A	100%	intergenic (‑139/+32)	ydgF ← / ← dinB	putative amino acid permease/nuclease inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	608,214	1	.	A	90.8%	1042.3 / 79.4	305	intergenic (‑139/+32)	ydgF/dinB	putative amino acid permease/nuclease inhibitor
Reads supporting (aligned to +/- strand): ref base . (5/23); new base A (173/104); total (178/127)
Fisher's exact test for biased strand distribution p-value = 8.10e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 334 total aligned reads displayed.

GTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑G‑AA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCTGTGATATTGCCGCGGTGGTATGTCCCGT  >  BS168_reference_genome/608067‑608358

  cTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑G‑aa‑aa                                                                                                                                                >  1:1631309/1‑150 (MQ=255)

   tGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGCTTATTATAACACTAAAATTTCTATCTTTTTTACTAAAATAGCATTGATGAAAAGAATTTTAATGTGATTTGAACTATTTCCC‑G‑ta‑ata                                                                                                                                               >  2:1528066/1‑148 (MQ=255)

   tGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑G‑aa‑aaa                                                                                                                                               >  1:4614139/1‑150 (MQ=255)

    gTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTAATATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGCCTATTTCCC‑G‑aa‑aaaa                                                                                                                                              <  1:1779570/150‑2 (MQ=255)

     tCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑G‑aa‑aaaag                                                                                                                                             <  2:2589674/150‑3 (MQ=255)

      cACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑G‑aa‑aaaagc                                                                                                                                            >  1:2582473/1‑147 (MQ=255)

        cTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑G‑aa‑aaaagcgg                                                                                                                                          >  2:4276392/1‑145 (MQ=255)

        cTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑G‑aa‑aaaaaagg                                                                                                                                          >  1:1277738/1‑145 (MQ=255)

            aTGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTACTTCCC‑GAAA‑AAA‑GC‑GGCTcc                                                                                                                                     <  2:3360026/150‑1 (MQ=255)

             tgtgTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCa                                                                                                                                    >  2:4880839/1‑150 (MQ=255)

              gtgtTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCAt                                                                                                                                   <  2:5686693/150‑1 (MQ=255)

                gtTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATTTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCAttt                                                                                                                                 <  1:361575/150‑1 (MQ=255)

                   cctccATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑AC‑GGCTCCATTTAgg                                                                                                                              >  2:3143492/1‑150 (MQ=255)

                     tccATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGa                                                                                                                            >  2:2364975/1‑150 (MQ=255)

                        aTTGAAAATTGTTTATTGGGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑CC‑GGCACCATATAGGGAGcc                                                                                                                         >  2:5270238/1‑150 (MQ=255)

                        aTTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGcc                                                                                                                         >  2:3500234/1‑150 (MQ=255)

                        aTTGAAAATTGTTTATTGAGACAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGcc                                                                                                                         <  2:2340705/150‑1 (MQ=255)

                         ttGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCg                                                                                                                        <  1:4032918/150‑1 (MQ=255)

                          tGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAATATTTCTATTTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGc                                                                                                                       <  1:2861420/150‑1 (MQ=255)

                            aaaaTTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTTTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCtt                                                                                                                     <  1:5193682/150‑1 (MQ=255)

                            aaaaTTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑CGATCCATTTAGGGAGccacct                                                                                                                     >  2:4118914/1‑146 (MQ=255)

                                ttGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGGTCCATTTAGGGAGCCGCtttttt                                                                                                                 >  1:933871/1‑150 (MQ=255)

                                ttGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCtttttt                                                                                                                 >  1:5287545/1‑150 (MQ=255)

                                 tGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTa                                                                                                                >  1:5791780/1‑150 (MQ=255)

                                 tGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑G‑AA‑AAACGC‑GGCTCCATTTAGGGAGCCGCTTTTTTa                                                                                                                <  2:3745644/150‑1 (MQ=255)

                                   tttATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑G‑AA‑AAAGGC‑GGCTCCATTTAGGGAGCCGCTTTTTTATg                                                                                                              <  2:73257/150‑1 (MQ=255)

                                    ttATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGc                                                                                                             <  2:5405221/150‑1 (MQ=255)

                                     tATTGAGTCAGCAATCTATTTTTAATGATTATTATAACATTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCt                                                                                                            <  2:3076717/150‑1 (MQ=255)

                                      aTTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCtt                                                                                                           >  2:200508/1‑150 (MQ=255)

                                       ttGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCttt                                                                                                          >  2:3666172/1‑150 (MQ=255)

                                       ttGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GA‑GGCTCCATTTAGGGAGCCGCTTTTTTTTGCttt                                                                                                          >  1:2643748/1‑150 (MQ=255)

                                           gTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑G‑AA‑AAA‑ACAGGCTCCATTTAGGGAGCCGCTTTTTTATGGTTTTTCa                                                                                                      >  1:1241149/1‑150 (MQ=255)

                                              aGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGtt                                                                                                   <  2:3867054/150‑1 (MQ=255)

                                               gCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGttt                                                                                                  <  1:289178/150‑1 (MQ=255)

                                                cAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGtttt                                                                                                 <  2:5832070/150‑1 (MQ=255)

                                                    cTATTTTTAATGATTATTATAACACTAAAATTTCTATTTTTTTTACTAAAATAGGATTAATAAAAATAATTTTAATGTGATTTGGACTAATTCCC‑C‑AA‑AAAGGG‑GGCCCCTTTTAGGAAGCCGCTTTTATATGCTTTTTCAGTTTTTGTc                                                                                             <  1:757695/150‑1 (MQ=255)

                                                          ttAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAAc                                                                                       >  2:552669/1‑150 (MQ=255)

                                                          ttAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCACTTTTTTATGCTTTTTCAGTTTTTGTCATAAAc                                                                                       <  2:4839968/150‑1 (MQ=255)

                                                            aaTGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTCCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAg                                                                                     <  2:4396343/150‑1 (MQ=255)

                                                             aTGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAgg                                                                                    >  2:5617536/1‑150 (MQ=255)

                                                              tGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GA‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGt                                                                                   >  1:1961523/1‑150 (MQ=255)

                                                                attattATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTag                                                                                 >  2:2875618/1‑150 (MQ=255)

                                                                   attatATCACTAAAATTTCTATCTTTTTTATTAAAAAAGGATTAATGAAAAGAATTTTAATGTGATTTGGACTATTTCCCAA‑AA‑AAA‑TC‑GGCTTCATTTAGGGAGCCCTTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGc                                                                              <  1:3103973/150‑1 (MQ=255)

                                                                   attatAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGc                                                                              >  1:983206/1‑150 (MQ=255)

                                                                    ttatAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAAGTAGTGcc                                                                             >  2:2522392/1‑150 (MQ=255)

                                                                     tataACACTAAAATTTCTATCTTTTTTACTAAAAAAGGTTTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCCGTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGccc                                                                            <  2:4249271/150‑1 (MQ=255)

                                                                       taAAACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGt                                                                          <  1:5151036/150‑1 (MQ=255)

                                                                         acacTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTaa                                                                        >  2:5078147/1‑150 (MQ=255)

                                                                         acacTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTaa                                                                        >  1:3268508/1‑150 (MQ=255)

                                                                           acTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GA‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTTATAAACAGGTAGAGCCCGTAATc                                                                      >  1:535760/1‑150 (MQ=255)

                                                                             tAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGGTCCATTTAGGGAGCCGCATTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGAAATCTg                                                                    >  1:2213012/1‑150 (MQ=255)

                                                                             tAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTg                                                                    >  1:2476698/1‑150 (MQ=255)

                                                                               aaaTTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTc                                                                  <  1:2915571/150‑1 (MQ=255)

                                                                                 aTTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTGCCC‑GAAA‑AAA‑GC‑GGCTCCATTTACGGAGCCACTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCgg                                                                <  1:4118914/150‑1 (MQ=255)

                                                                                      tATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑G‑AA‑AAAGGC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCgg                                                           <  2:3309724/150‑1 (MQ=255)

                                                                                       aTCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGa                                                          >  1:5081630/1‑150 (MQ=255)

                                                                                       aTCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GA‑TGATCCATTTCGGGAGGCGCGTTTTTAAGCTTTTTCAGTTTTTGTCCAAATCCTGGAGAGGCCCCAATCAGTCGGTgtcgt                                                          >  1:1514523/1‑146 (MQ=255)

                                                                                          ttttttACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GG‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGGCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGc                                                       >  1:230026/1‑150 (MQ=255)

                                                                                          ttttttACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTCCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCTCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGc                                                       <  2:4940824/150‑1 (MQ=255)

                                                                                               tACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAgc                                                  >  2:2936930/1‑150 (MQ=255)

                                                                                                aCTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAgcc                                                 >  1:4224467/1‑150 (MQ=255)

                                                                                                 taaaaaaaGGTTTAATAAAAAGAATTTTAATGTGATTTGAATTATTTCCCAA‑AA‑AAA‑GC‑GGTTTCTTTTAGGGAGCCTTTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAgcct                                                <  1:4432088/148‑1 (MQ=255)

                                                                                                 cTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCCTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGGTGCGTAgcct                                                >  2:2695609/1‑150 (MQ=255)

                                                                                                  tAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCCTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAgcctg                                               >  1:1052547/1‑150 (MQ=255)

                                                                                                  tAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑G‑AA‑AAAGGC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAgcctg                                               <  2:1612219/150‑1 (MQ=255)

                                                                                                   aaaaaaGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAgcctgc                                              >  2:1258392/1‑150 (MQ=255)

                                                                                                    aaaaaGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAgcctgcc                                             >  2:4454198/1‑150 (MQ=255)

                                                                                                      aaaGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCATTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGcctgcctg                                           >  1:3175478/1‑150 (MQ=255)

                                                                                                       aaGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCctgcctgc                                          >  1:3756172/1‑150 (MQ=255)

                                                                                                        aGGTTTGATAAAAAAATTTTTAATGTTATTTGCATTATTCCAC‑G‑AA‑AAAGGC‑GGCTCAATTTAGGGAGCCTCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCtgcctgcc                                         <  1:5008935/150‑1 (MQ=255)

                                                                                                        aGGATTGATGAAAAGATTTTTAATTTAATTTGAACTATTCACC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGTTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCtgcctgcc                                         <  1:1010541/150‑1 (MQ=255)

                                                                                                         ggATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCTCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCg                                        <  1:2364975/150‑1 (MQ=255)

                                                                                                          gATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGATCCATTTAGGGAGCCGCCTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCgc                                       >  1:4333525/1‑150 (MQ=255)

                                                                                                              gatgaAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCgcagca                                   >  1:3876171/1‑150 (MQ=255)

                                                                                                                 gaAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTg                                >  1:598907/1‑150 (MQ=255)

                                                                                                                 gaAAAGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTg                                >  2:2120215/1‑150 (MQ=255)

                                                                                                                  aaaaGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTTGTGCCCGGAAGCTGTCGGGGCGGAGGCGTAGCCTGCCTGCCGCAGCAATGc                               >  2:136102/1‑150 (MQ=255)

                                                                                                                   aaaGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCt                              >  2:4016335/1‑150 (MQ=255)

                                                                                                                   aaaGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCt                              >  2:5779415/1‑150 (MQ=255)

                                                                                                                   aaaGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCATTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCt                              >  1:1594192/1‑150 (MQ=255)

                                                                                                                    aaGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GG‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTTGGTGCCGATGCGTAGCCTGCCTGCCGCAGCATTGCtg                             >  1:4784738/1‑150 (MQ=255)

                                                                                                                     aGAATTTTAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCtgt                            <  2:2476698/150‑1 (MQ=255)

                                                                                                                         ttttAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCTGTGata                        >  2:2753255/1‑150 (MQ=255)

                                                                                                                         ttttAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCTGTGata                        <  2:538123/150‑1 (MQ=255)

                                                                                                                          tttAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCTGTGatat                       <  2:1204471/150‑1 (MQ=255)

                                                                                                                           ttAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCTGTGATAtt                      >  2:2165986/1‑150 (MQ=255)

                                                                                                                           ttAATGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCTGTGATAtt                      <  2:1664777/150‑1 (MQ=255)

                                                                                                                             atttttaTTTGGACTTTTTCCC‑G‑AA‑AAAGGC‑GTCCCCATTTAGGGAGCCACTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCTGTGATATTGc                    <  1:2803010/144‑1 (MQ=255)

                                                                                                                             aaTGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCTGTGATATTGc                    >  2:974700/1‑150 (MQ=255)

                                                                                                                              aTGTGATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCTGTGATATTGcc                   >  1:2088575/1‑150 (MQ=255)

                                                                                                                              aTGTGATTTGGACTATTTCCC‑G‑AAGAAA‑GG‑GGCTCCAATTAGGGAGCCGCATTTTTACGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGTCCGTGCTCTGGCGGTGCGGGTGCGTAGGCTGCCTGCCGCAGCATGGCTGTGATACTGcc                   >  2:4414187/1‑150 (MQ=255)

                                                                                                                               tgtgATTTGGACTATTTGCA‑G‑AA‑AAATGT‑GGTTCCTTTTAGGGAGCCTTTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGTTGCGGATGCGTAACCTGCCTGCCGCAGCATTGCTGTGATATTGCcg                  <  1:1528070/150‑1 (MQ=255)

                                                                                                                                 tgATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCTGTGATATTGCcgcg                >  1:3409252/1‑150 (MQ=255)

                                                                                                                                 tgATTTGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCTGTGATATTGCcgcg                >  2:4761404/1‑150 (MQ=255)

                                                                                                                                      tGGACTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTCTTTTATGCTTGTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGGTGCGTAGCCTGCCTGCCGCAGCATTGCTGTGATATTGCCGCggtggt           >  2:3414647/1‑150 (MQ=255)

                                                                                                                                         aCTATTTCCG‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCACCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCTGTGATATTGCCGCGGTGGTATg        <  1:4880839/150‑1 (MQ=255)

                                                                                                                                         aCTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCTGTGATATTGCCGCGGTGGTATg        >  1:5597875/1‑150 (MQ=255)

                                                                                                                                          cTATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCTGTGATATTGCCGCGGTGGTATGt       <  2:4160753/150‑1 (MQ=255)

                                                                                                                                           tATTTCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCTGTGATATTGCCGCGGTGGTATGTc      >  2:82473/1‑150 (MQ=255)

                                                                                                                                             tttCCC‑GAAA‑AAA‑GC‑GGCTCCATTTAGGGAGCCTCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCTGTGATATTGCCGCGGTGGTATGTccc    <  2:710950/150‑1 (MQ=255)

GTCTGTCACTTTATGTGTTCCTCCATTGAAAATTGTTTATTGAGTCAGCAATCTATTTTTAATGATTATTATAACACTAAAATTTCTATCTTTTTTACTAAAAAAGGATTGATGAAAAGAATTTTAATGTGATTTGGACTATTTCCC‑G‑AA‑AAA‑GC‑GGCTCCATTTAGGGAGCCGCTTTTTTATGCTTTTTCAGTTTTTGTCATAAACAGGTAGAGCCCGTAATCTGTCGGTGCGGATGCGTAGCCTGCCTGCCGCAGCATTGCTGTGATATTGCCGCGGTGGTATGTCCCGT  >  BS168_reference_genome/608067‑608358

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