BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	BS168_reference_genome	165,826	+C	intergenic (+1/‑4)	trnI‑Asn → / → trnI‑Thr	tRNA‑Asn/tRNA‑Thr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	165,825	1	.	C	98.5%	675.5 / 3.7	202	noncoding (72/72 nt)	trnI‑Asn	tRNA‑Asn
Reads supporting (aligned to +/- strand): ref base . (2/1); new base C (97/102); total (99/103)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 213 total aligned reads displayed.

GTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATCTCATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAG‑C‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAG  >  BS168_reference_genome/165674‑165968

gTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAg                                                                                                                                                   >  1:2783808/1‑150 (MQ=255)

 tttCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAG‑c                                                                                                                                                 <  2:2910327/150‑1 (MQ=255)

  ttCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAG‑cc                                                                                                                                                >  1:697927/1‑149 (MQ=255)

  ttCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAG‑cc                                                                                                                                                <  1:338396/150‑2 (MQ=255)

     ccccTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAG‑ccatt                                                                                                                                             <  2:3333904/150‑5 (MQ=255)

      cccTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAG‑ccattg                                                                                                                                            >  1:1251693/1‑145 (MQ=21)

      cccTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAG‑ccattg                                                                                                                                            <  2:2393982/150‑6 (MQ=21)

      cccTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAG‑ccattg                                                                                                                                            <  2:2644691/150‑6 (MQ=21)

         tgtgAGAGCAGGTCGCCGCCAAGCAGATTTATTGGTTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCCTAGTGGATGAGCGATCGGCTGTTAACCAGTCGGACGTAGGTTCGACTCCTCCCTG‑TGGAGCC‑ATTGCtt                                                                                                                                         <  2:2342342/150‑1 (MQ=255)

          gtgAGAGTAGGACTCCGCCAAGCAGATTTATTAATTCTACCCTGTCTGATTATGAAGCTTTATAT‑T‑ATCCCGCAGTAGGTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTc                                                                                                                                        <  2:3662699/150‑1 (MQ=18)

          gtgAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTc                                                                                                                                        >  1:394541/1‑150 (MQ=21)

           tgAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGGAGGTTCGGATCCTACCTG‑CGGAGCC‑ATTGCTTcc                                                                                                                                       >  2:114545/1‑150 (MQ=21)

            gagaGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCa                                                                                                                                      <  1:1193069/150‑1 (MQ=21)

            gagaGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCa                                                                                                                                      >  2:4008389/1‑150 (MQ=21)

             tgagTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTTTAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGACGTAGGTTCGAAACCTACCTG‑CGGAGCC‑ATTGCTTCCAt                                                                                                                                     <  1:3293828/149‑1 (MQ=21)

              gagTAGGCCGCCGCCAAGCAGATTTATTAGTTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATa                                                                                                                                    <  2:4313303/150‑1 (MQ=21)

              gagTAGGACCCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATa                                                                                                                                    <  2:1602793/150‑1 (MQ=21)

                gTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGc                                                                                                                                  >  1:4690488/1‑150 (MQ=21)

                gTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCACAGc                                                                                                                                  >  1:1513171/1‑150 (MQ=21)

                  aGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTGCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTc                                                                                                                                <  1:2929452/150‑1 (MQ=21)

                     aCGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTcagc                                                                                                                             <  2:699614/150‑1 (MQ=21)

                      cgccgcCAAGCAGATTCATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATCCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGCCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTcagca                                                                                                                            <  1:780276/150‑1 (MQ=21)

                          gccAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTa                                                                                                                        <  2:1496279/150‑1 (MQ=21)

                            cAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTaga                                                                                                                      <  2:4505331/150‑1 (MQ=21)

                            cAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGGCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTaga                                                                                                                      >  2:3978542/1‑150 (MQ=21)

                             aaGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGGagag                                                                                                                     >  2:127991/1‑150 (MQ=21)

                                cAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCAc                                                                                                                  <  2:607150/150‑1 (MQ=21)

                                    tttATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCCTAGCTCAGCAGGTAGAGCACTTcc                                                                                                              >  2:2458135/1‑150 (MQ=21)

                                    tttATTAATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTcc                                                                                                              >  1:4188550/1‑150 (MQ=21)

                                     ttattaATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCa                                                                                                             <  1:1367983/150‑1 (MQ=21)

                                      tattaATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCAt                                                                                                            <  2:2254279/150‑1 (MQ=21)

                                        ttaATTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATgg                                                                                                          >  1:4761915/1‑150 (MQ=21)

                                          aaTTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAGTCCTACCTGCCGGGC‑C‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCAGGGTa                                                                                                        >  2:3787160/1‑150 (MQ=21)

                                           aTTCTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGCATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCGGCAGGTAGAGCACTTCCATGATaa                                                                                                       >  1:29921/1‑150 (MQ=21)

                                            ttCTGCCATGTTTGTTTATGAAGGTTT‑T‑TTTTATTCCGCAGTAGCTCAGTGGTAGAGCAATCGGCTGTTAACCGATCGGTTGTAGGTTCGAATCCTTCCTG‑CGGAG‑CTATTGCTTCCATACTTCAGCAGGTAGTGCACTTCCATGGTAAg                                                                                                      <  1:3760889/150‑1 (MQ=18)

                                            ttCAACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCCCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAg                                                                                                      <  1:585537/150‑1 (MQ=21)

                                              cTACCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCCCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGa                                                                                                    <  1:2032892/150‑1 (MQ=21)

                                                aCCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAag                                                                                                  >  2:4855209/1‑150 (MQ=21)

                                                aCCATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAag                                                                                                  <  1:2821182/150‑1 (MQ=21)

                                                  cATGTTTGTTTATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGGTCGGTCGTAGGTTCGAATCCTAGCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCTCTTCCATGGTAAGGAagag                                                                                                <  1:114545/150‑1 (MQ=21)

                                                     gtttgtttATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTc                                                                                             <  2:1780936/150‑1 (MQ=21)

                                                     gtttgtttATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTc                                                                                             >  1:1719678/1‑150 (MQ=21)

                                                      tttgtttATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCa                                                                                            <  1:406955/150‑1 (MQ=21)

                                                        tgtttATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGc                                                                                          >  1:2535437/1‑150 (MQ=21)

                                                         gtttCTGAAGCTTTTTAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCg                                                                                         <  1:2393129/150‑1 (MQ=21)

                                                         gtttATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCg                                                                                         <  1:3893409/150‑1 (MQ=21)

                                                         gtttATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCg                                                                                         <  1:3481545/150‑1 (MQ=21)

                                                         gtttATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCCCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCg                                                                                         <  1:4446549/150‑1 (MQ=21)

                                                          tttATGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGGAGAGGTCAGCgg                                                                                        >  1:2172604/1‑150 (MQ=21)

                                                            tATGAAGCTTTATAT‑T‑ATGCCGCAGTAGCTCAGTGGTAGAGCTACCGGCTGTTAACCGATCGGTCGTAGGTCCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGtt                                                                                      <  1:3229088/150‑1 (MQ=21)

                                                             aTGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTc                                                                                     >  1:1458418/1‑150 (MQ=21)

                                                             aTGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGCCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTc                                                                                     <  1:4601973/150‑1 (MQ=21)

                                                              tGAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCg                                                                                    <  2:478632/150‑1 (MQ=21)

                                                               gAAGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTATCCGATCGGTCGTAGGTCCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGa                                                                                   <  2:307290/150‑1 (MQ=21)

                                                                 aGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGc                                                                                 >  1:374077/1‑150 (MQ=21)

                                                                 aGCTTTATAT‑T‑ATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGc                                                                                 <  1:1222292/150‑1 (MQ=21)

                                                                       ttatattaTTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTTACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGc                                                                             <  2:2197053/143‑1 (MQ=25)

                                                                       ttatattaTTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGGGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGc                                                                             >  2:2488056/8‑150 (MQ=21)

                                                                          tattaTTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTACCCGATCGGTGGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTg                                                                          <  1:4119574/146‑1 (MQ=25)

                                                                          tattaTTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACTTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTg                                                                          <  2:3429892/146‑1 (MQ=255)

                                                                          tattaTTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTg                                                                          >  1:3551026/5‑150 (MQ=255)

                                                                          tattaTTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTCCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTg                                                                          <  2:1259650/146‑1 (MQ=255)

                                                                                    cAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTaaa                                                                >  1:4829251/1‑150 (MQ=255)

                                                                                    cAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTaaa                                                                >  1:1906304/1‑150 (MQ=255)

                                                                                            cAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGtattat                                                        >  2:2550888/1‑150 (MQ=255)

                                                                                             aGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTattatt                                                       >  1:5201/1‑150 (MQ=255)

                                                                                             aGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTattatt                                                       >  1:1088273/1‑150 (MQ=255)

                                                                                              gTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTAttatta                                                      >  2:1392123/1‑150 (MQ=255)

                                                                                              gTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTAttatta                                                      >  1:2283978/1‑150 (MQ=255)

                                                                                                       cTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGttt                                             <  2:827523/150‑1 (MQ=255)

                                                                                                       cTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGttt                                             <  2:173253/150‑1 (MQ=255)

                                                                                                        tATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTtc                                            >  1:2634885/1‑150 (MQ=255)

                                                                                                        tATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTtc                                            >  2:3206863/1‑150 (MQ=255)

                                                                                                         aTCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTtct                                           <  2:1652615/150‑1 (MQ=255)

                                                                                                         aTCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTtct                                           <  2:2083643/150‑1 (MQ=255)

                                                                                                          tCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTtctc                                          >  1:4619257/1‑150 (MQ=255)

                                                                                                           cGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCa                                         <  2:1071617/150‑1 (MQ=255)

                                                                                                           cGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCa                                         <  2:1339808/150‑1 (MQ=255)

                                                                                                            ggCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCAt                                        <  2:3034881/150‑1 (MQ=255)

                                                                                                            ggCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCAt                                        <  2:3353686/150‑1 (MQ=255)

                                                                                                              cTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATaa                                      <  2:392764/150‑1 (MQ=255)

                                                                                                                      cGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAg                              >  1:130623/1‑150 (MQ=255)

                                                                                                                       gATCGGCCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGc                             <  2:4172233/150‑1 (MQ=255)

                                                                                                                         ccGGCCGTAGGTCCGAGCCCTGCCCG‑GGGAGCC‑ATTGCTGCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCtt                           <  2:29921/149‑1 (MQ=255)

                                                                                                                          cGGTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGGGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCttt                          >  2:3080359/1‑150 (MQ=255)

                                                                                                                           ggTCGTAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGGGAAAGCtgtt                         >  2:3315994/1‑150 (MQ=255)

                                                                                                                               gcaGGTCCGAATCCTCCCTG‑CGGAGCC‑ATTGCTTCCATAGCCCAGCGGGTAGAGCCCCTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTa                     <  2:3403137/148‑1 (MQ=255)

                                                                                                                                tAGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTAt                    <  2:525359/150‑1 (MQ=255)

                                                                                                                                 aGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAGGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGGGAAGGCTTTTTTTAtt                   >  2:2004613/1‑150 (MQ=255)

                                                                                                                                 aGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTAtt                   >  1:2449917/1‑150 (MQ=255)

                                                                                                                                 aGGTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTAtt                   <  2:1759403/150‑1 (MQ=255)

                                                                                                                                  ggTTCGAATCCTCCCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTg                  <  2:4619037/150‑1 (MQ=255)

                                                                                                                                  ggTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTGAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTg                  >  2:612948/1‑150 (MQ=255)

                                                                                                                                  ggTTCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGCAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGACAAAGCTTTTTTTATTg                  >  1:3818099/1‑150 (MQ=255)

                                                                                                                                    ttCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCg                >  2:297112/1‑150 (MQ=255)

                                                                                                                                     tCGAATCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGa               >  1:4105730/1‑150 (MQ=255)

                                                                                                                                         aTCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAAGAGGAAGCTTTTTTTATTGCGAGATg           >  1:1985654/1‑150 (MQ=255)

                                                                                                                                          tCCTACCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGc          >  1:1614534/1‑150 (MQ=255)

                                                                                                                                           ccTCCCTG‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCg         <  2:4459321/150‑1 (MQ=255)

                                                                                                                                                  g‑CGGAGCC‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAagtag  >  2:2306970/1‑150 (MQ=255)

GTTTCCCCCTGTGAGAGTAGGACGCCGCCAAGCAGATTTATTAATTCTACCATGTTTGTTTATGAAGCTTTATATCTCATTCCGCAGTAGCTCAGTGGTAGAGCTATCGGCTGTTAACCGATCGGTCGTAGGTTCGAATCCTACCTG‑CGGAG‑C‑ATTGCTTCCATAGCTCAGCAGGTAGAGCACTTCCATGGTAAGGAAGAGGTCAGCGGTTCGAGCCCGCTTGGAAGCTTAAATGTATTATTACCAAGGTTTCTCATAAGGAGAAAGCTTTTTTTATTGCGATATGCGGAAGTAG  >  BS168_reference_genome/165674‑165968

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