BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,595,328	(G)_8→9	noncoding (72/87 nt)	leuV ←	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,595,320	1	.	G	91.7%	323.3 / 24.6	96	noncoding (80/87 nt)	leuV	tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base . (2/6); new base G (54/34); total (56/40)
Fisher's exact test for biased strand distribution p-value = 6.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 128 total aligned reads displayed.

CTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAT  >  REL606/4595245‑4595385

cTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑Ga                                                                      <  1:496081/76‑1 (MQ=255)

cTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑Ga                                                                      <  2:3628180/76‑1 (MQ=255)

cTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑Ga                                                                      <  2:1758885/76‑1 (MQ=255)

cTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑Ga                                                                      >  1:937380/1‑76 (MQ=255)

 tgatgaTAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑g                                                                   <  1:4075478/76‑1 (MQ=255)

 tgatgaTAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑g                                                                   <  2:1506363/76‑1 (MQ=255)

  gatgatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑gg                                                                  <  1:2769268/76‑1 (MQ=255)

   atgatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑ggg                                                                 >  1:6061351/1‑76 (MQ=255)

      atAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑gggggg                                                              >  2:748133/1‑76 (MQ=255)

       tAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑ggggggg                                                             >  1:519600/1‑76 (MQ=255)

       tAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑ggggggg                                                             >  2:1476500/1‑76 (MQ=255)

        aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑gggggggg                                                            <  1:212378/76‑1 (MQ=255)

        aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑gggggggg                                                            <  2:244288/76‑1 (MQ=255)

         aTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑ggggggggg                                                           <  2:398579/76‑2 (MQ=255)

         aTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑ggggggggg                                                           <  2:3033696/76‑2 (MQ=255)

         aTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑ggggggggg                                                           <  2:587466/76‑2 (MQ=255)

          tGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑ggggggggga                                                          <  2:1863830/76‑3 (MQ=255)

          tGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑ggggggggga                                                          <  2:5645735/76‑3 (MQ=255)

          tGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑ggggggggga                                                          >  2:2676152/1‑74 (MQ=255)

          tGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑ggggggggga                                                          >  1:4762992/1‑74 (MQ=255)

           gAATTTGAATCAGGGATGAAAAATCAAGGAAGAAAAAAGAAAGAAAGTAAAGATAATTGGT‑GC‑GA‑‑gggggggggac                                                         <  1:6277857/76‑4 (MQ=255)

             aTTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑GGGGGGGG‑ACTTg                                                      <  1:130964/76‑1 (MQ=255)

              tttGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTg                                                      >  1:1130084/1‑76 (MQ=255)

               ttGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGa                                                     <  1:4543946/76‑1 (MQ=255)

               ttGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGa                                                     >  1:5414864/1‑76 (MQ=255)

                tGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGaa                                                    >  1:1473986/1‑76 (MQ=255)

                tGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGaa                                                    <  2:3986136/76‑1 (MQ=255)

                tGAATAAGGAATGAAAATTAAAGAAAAAAACAAGAAAGGAAGAAAAGAAAATTGGG‑GCGGG‑‑GGGGGGGG‑ACTTGaa                                                    <  2:812345/76‑1 (MQ=255)

                 gAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAAc                                                   <  1:2054468/76‑1 (MQ=255)

                 gAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAAc                                                   >  1:5664100/1‑76 (MQ=255)

                  aaTCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAAcc                                                  >  1:1325354/1‑76 (MQ=255)

                  aaTCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAAcc                                                  <  2:3372456/76‑1 (MQ=255)

                  aaTCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAAcc                                                  >  2:433283/1‑76 (MQ=255)

                   aTCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAAccc                                                 >  2:1109361/1‑76 (MQ=255)

                     cAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAAccccc                                               <  2:81588/76‑1 (MQ=255)

                     cAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAAccccc                                               >  1:2173769/1‑76 (MQ=255)

                     cAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAAccccc                                               <  2:5738589/76‑1 (MQ=255)

                      aGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAATTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCa                                              >  1:3157125/1‑76 (MQ=38)

                      aGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGGACACCCa                                              >  2:535916/1‑76 (MQ=255)

                      aGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCa                                              >  1:513494/1‑76 (MQ=255)

                      aGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCa                                              <  1:4702280/76‑1 (MQ=255)

                      aGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCa                                              >  1:5494734/1‑76 (MQ=255)

                      aGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCa                                              <  2:345905/76‑1 (MQ=255)

                      aGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCa                                              <  1:4675269/76‑1 (MQ=255)

                      aGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCa                                              >  1:1568367/1‑76 (MQ=255)

                       gggATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCAc                                             <  1:4726236/76‑1 (MQ=255)

                       gggATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCAc                                             <  1:3212252/76‑1 (MQ=255)

                       gggATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCAc                                             <  1:1535540/76‑1 (MQ=255)

                       gggATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCAc                                             <  2:1110929/76‑1 (MQ=255)

                        ggATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAGGGGGGGGGG‑ACTTGAACCCCCAc                                             >  1:1067541/1‑76 (MQ=38)

                        ggATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACg                                            <  2:4728681/76‑1 (MQ=255)

                         gATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGt                                           >  1:3256241/1‑76 (MQ=255)

                          aTGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTc                                          <  1:1266173/76‑1 (MQ=255)

                           tGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTcc                                         <  2:1400311/76‑1 (MQ=255)

                           tGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTcc                                         <  2:1890379/76‑1 (MQ=255)

                            gAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCg                                        >  1:416971/1‑76 (MQ=255)

                            gAAAAAAAAAGAAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCg                                        <  2:3167534/76‑1 (MQ=25)

                             aaaaaTAAAGGAAAAAAAAAAAAAGGAAGAAAATATAATTGGTGGG‑GG‑‑GGGGGGGG‑TCTTGACCCCCCACGTCCGt                                       <  2:65304/76‑1 (MQ=16)

                              aaaaTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑CCTTGAACCCCCACGCCCGTa                                      >  1:543968/1‑76 (MQ=37)

                              aaaaTCAAGAAAGAAACAAGAAAGAAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGACCCCCCACGTCCGTa                                      <  2:2462461/76‑1 (MQ=25)

                               aaaTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTaa                                     <  2:2587996/76‑1 (MQ=255)

                               aaaTCAAGAAAGAAACAAGAAAGGAAGTAAAGATAATGGGT‑GC‑GAG‑GGGGGGGG‑ACTTGACCCCCCACGTCCGTaa                                     <  1:6005387/76‑1 (MQ=25)

                               aaaGCAAGAAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGACCCCCCACGTCCGTaa                                     <  1:1138293/76‑1 (MQ=25)

                                agccaagaaagaaacaagaaagaaaTAAAAGATAATGGGG‑CC‑GA‑‑GGGGGGGGTATTTGACCCCCCACGTCCGTAAg                                    <  2:3263190/73‑1 (MQ=16)

                                aaTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GG‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAg                                    >  1:5307296/1‑76 (MQ=37)

                                aaTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAg                                    >  1:1029081/1‑76 (MQ=255)

                                aaTCAAGAAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGACCCCCCACGTCCGTAAg                                    <  2:4089296/76‑1 (MQ=25)

                                  tCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGa                                  >  1:2048040/1‑76 (MQ=255)

                                  ccAAGGAAGAAACAAGAAAGGAAGTAAAGATAATCGGG‑GC‑GAG‑GGGGGGGG‑TCTTGACCCCCCACGTCCGTAAGGa                                  <  2:3041202/75‑1 (MQ=25)

                                   cAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAgggc                                 >  1:4692634/1‑74 (MQ=25)

                                      ggAAGAAACAAAAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑CCTTGACCCCCCACGTCCGTAAGGACACt                              <  2:1327753/76‑1 (MQ=21)

                                        aagaaACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTaa                            >  2:1616557/1‑76 (MQ=255)

                                        aagaaACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTaa                            <  1:6264265/76‑1 (MQ=255)

                                        aagaaACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTaa                            >  1:3561910/1‑76 (MQ=255)

                                          gaaACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGAAACTAaca                          >  2:2776716/1‑76 (MQ=21)

                                               aaGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑GCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGa                     >  2:3589381/1‑76 (MQ=18)

                                               aaGAAAGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGa                     >  2:2288770/1‑76 (MQ=255)

                                                  aaaGGAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGc                  >  1:6003438/1‑76 (MQ=255)

                                                  aaaGGAAGTAAAGATAATTGGT‑GC‑GA‑‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCt                 >  2:4229909/1‑76 (MQ=255)

                                                    aGGAAGTAAAGATAATTGGT‑GCGGG‑‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTa                <  2:2503360/76‑1 (MQ=17)

                                                     ggAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg               >  1:692150/1‑76 (MQ=255)

                                                      gAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgc              <  2:1090793/76‑1 (MQ=17)

                                                      gAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgc              >  2:295963/1‑76 (MQ=255)

                                                      gAAGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgc              >  2:411797/1‑76 (MQ=255)

                                                       aaGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcg             >  2:350318/1‑76 (MQ=255)

                                                       aaGTAAAGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcg             >  1:3409329/1‑76 (MQ=255)

                                                           aaaGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTc         >  2:2196071/1‑76 (MQ=255)

                                                            aaGATAATTGGT‑GC‑GA‑‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTa       >  2:5871042/1‑76 (MQ=40)

                                                             aGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTa       >  1:4844166/1‑76 (MQ=255)

                                                             aGATAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTa       >  1:2388475/1‑76 (MQ=255)

                                                              gagaaTGGGT‑GC‑GAG‑GGGGGGGG‑ACGTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAc      <  2:4723595/73‑1 (MQ=11)

                                                              aaTAATTGGG‑GC‑GAG‑GGGGGGGG‑ACTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAc      <  2:5327534/75‑1 (MQ=11)

                                                               aTAATTGGT‑GC‑GAG‑GGGGGGGG‑CCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGCCTAca     >  2:2947101/1‑75 (MQ=11)

                                                               aTAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAcc     >  2:2767841/1‑76 (MQ=255)

                                                                tAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCa    <  2:4246683/76‑1 (MQ=11)

                                                                tAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCa    >  1:1735927/1‑76 (MQ=255)

                                                                tAATTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCa    >  2:4526521/1‑76 (MQ=255)

                                                                 aaTTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa   <  2:4678061/76‑1 (MQ=11)

                                                                 aaTTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa   >  1:3869065/1‑76 (MQ=255)

                                                                  aTTGGT‑GC‑GAG‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAt  >  2:1572751/1‑76 (MQ=9)

CTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGT‑GC‑GA‑‑GGGGGGGG‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAT  >  REL606/4595245‑4595385

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