BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,744,958	G→A	G82E (GGA→GAA)	lyxK →	L‑xylulose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,744,958	0	G	A	99.7%	1150.4 / NA	310	G82E (GGA→GAA)	lyxK	L‑xylulose kinase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (140/169); total (140/170)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 333 total aligned reads displayed.

GGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGGAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAATCGTTCGTC  >  NC_000913/3744868‑3745047

 gCAATGCTGCATGTCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGaaaa                                                                                        <  2:6102909/93‑1 (MQ=255)

 gCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGaaaa                                                                                        <  2:2350961/93‑1 (MQ=255)

 gCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGaaaa                                                                                        <  2:2760979/93‑1 (MQ=255)

  cAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAg                                                                                       <  2:4157542/93‑1 (MQ=255)

   aaTGCTGCATGGCTGTGATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAgg                                                                                      <  1:543359/93‑1 (MQ=255)

   aaTGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAgg                                                                                      >  2:3010029/1‑93 (MQ=255)

   aaTGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAgg                                                                                      >  2:2530664/1‑93 (MQ=255)

    aTGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAggg                                                                                     >  1:2091002/1‑93 (MQ=255)

     tgctgcATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGc                                                                                    <  1:4521238/93‑1 (MQ=255)

      gctgcATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCt                                                                                   <  2:1555592/93‑1 (MQ=255)

      gctgcATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCt                                                                                   <  1:1125786/93‑1 (MQ=255)

      gctgcATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCt                                                                                   <  2:5487304/93‑1 (MQ=255)

       ttgcATGGCTGTCATTCGCGCCCTTCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCtt                                                                                  <  1:1299582/92‑1 (MQ=255)

       ctgcATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCtt                                                                                  <  1:1403647/93‑1 (MQ=255)

       ctgcATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCtt                                                                                  <  2:6565625/93‑1 (MQ=255)

       ctgcATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCtt                                                                                  <  1:2565166/93‑1 (MQ=255)

        tgcATGTCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTg                                                                                 <  1:3251451/93‑1 (MQ=255)

        tgcATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTg                                                                                 <  2:1966614/93‑1 (MQ=255)

         gcATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGt                                                                                >  1:5957919/1‑93 (MQ=255)

         gcATCGCTGTCATTCGCGCCCTGCTTACCCATTCTGGTGTTAGCGGGGAACAAATTGTCGCTATCGGCATCTCCGCACAGGAAAAGGGCTTGt                                                                                <  1:1930454/93‑1 (MQ=255)

          cATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGtt                                                                               <  1:731928/93‑1 (MQ=255)

          cATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGtt                                                                               <  2:757751/93‑1 (MQ=255)

           aTGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGttt                                                                              <  2:6421812/93‑1 (MQ=255)

           aTGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGttt                                                                              >  2:5298964/1‑93 (MQ=255)

            tgtcTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGtttt                                                                             <  2:464969/90‑1 (MQ=255)

            tGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACTAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGtttt                                                                             <  1:3349984/93‑1 (MQ=255)

             ggCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGttttt                                                                            >  2:1158712/1‑93 (MQ=255)

             ggCTGTCATTCGCGCCATGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGttttt                                                                            <  1:2192235/93‑1 (MQ=255)

              gCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTg                                                                           >  2:4970623/1‑93 (MQ=255)

               cTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGc                                                                          >  2:5978378/1‑93 (MQ=255)

               cTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGc                                                                          >  1:4446075/1‑93 (MQ=255)

                 gTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTg                                                                        >  2:5651570/1‑93 (MQ=255)

                 gTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCATGTTTTTGCTg                                                                        >  1:788203/1‑93 (MQ=255)

                 atCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTg                                                                        >  2:5385920/2‑93 (MQ=255)

                  tCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTgg                                                                       <  1:3081433/93‑1 (MQ=255)

                  tCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTgg                                                                       <  2:4791202/93‑1 (MQ=255)

                    aTTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGAt                                                                     >  2:3407976/1‑93 (MQ=255)

                    aTTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGAt                                                                     <  1:3726993/93‑1 (MQ=255)

                    aTTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGAt                                                                     >  1:1212776/1‑93 (MQ=255)

                     ttCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATa                                                                    >  2:599065/1‑93 (MQ=255)

                        gcgcCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATaaaa                                                                 >  2:5014464/1‑93 (MQ=255)

                         cgcCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATaaaaa                                                                <  2:2847040/93‑1 (MQ=255)

                         cgcCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATaaaaa                                                                <  1:1067796/93‑1 (MQ=255)

                          gcCCTGCTTACTCATTCTGGTGTTAGCTGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAAc                                                               <  1:406192/93‑1 (MQ=255)

                            ccTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAACAACGa                                                             >  1:3616878/1‑93 (MQ=255)

                            ccTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGa                                                             <  2:1118109/93‑1 (MQ=255)

                            ccTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGa                                                             >  1:6670622/1‑93 (MQ=255)

                             cTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGAc                                                            <  2:1291131/93‑1 (MQ=255)

                              tGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACCGGAAAAGGGCTTGTTTTTGCTGGATAAAAAAGACa                                                           >  1:5938924/1‑93 (MQ=255)

                              tGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACa                                                           >  2:1697275/1‑93 (MQ=255)

                               gCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACaa                                                          >  1:6358667/1‑93 (MQ=255)

                               gCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACaa                                                          >  2:6282888/1‑93 (MQ=255)

                               gCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACaa                                                          >  1:3173265/1‑93 (MQ=255)

                                 ttACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAAc                                                        <  1:2358628/93‑1 (MQ=255)

                                    cTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGc                                                     >  1:1469193/1‑93 (MQ=255)

                                     tCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCt                                                    <  1:5681673/93‑1 (MQ=255)

                                     tCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCt                                                    >  2:3263431/1‑93 (MQ=255)

                                     tCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCt                                                    >  2:245600/1‑93 (MQ=255)

                                     tCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCt                                                    >  1:5594791/1‑93 (MQ=255)

                                      cATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTc                                                   <  2:4405705/93‑1 (MQ=255)

                                       atacTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCg                                                  <  2:5734446/90‑1 (MQ=255)

                                       aTTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCg                                                  >  1:2894525/1‑93 (MQ=255)

                                        ttCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCgg                                                 >  1:3573567/1‑93 (MQ=255)

                                         tCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCggg                                                <  2:5165913/93‑1 (MQ=255)

                                           tGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGaa                                              <  1:5186885/93‑1 (MQ=255)

                                            ggTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAACACGACAAACCGCTggggtgt                                             >  2:1368584/1‑90 (MQ=255)

                                            ggTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAAt                                             <  1:6271821/93‑1 (MQ=255)

                                            ggTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAAt                                             <  1:6467554/93‑1 (MQ=255)

                                              tgtTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGc                                           >  1:2634692/1‑93 (MQ=255)

                                               gtTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCt                                          <  1:2307129/93‑1 (MQ=255)

                                                ttAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTa                                         <  2:3959599/93‑1 (MQ=255)

                                                ttAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTa                                         <  1:2863695/93‑1 (MQ=255)

                                                  aGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTAtt                                       <  1:1794015/93‑1 (MQ=255)

                                                   gCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTAttt                                      >  1:4129962/1‑93 (MQ=255)

                                                    cGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTAtttt                                     >  1:2443265/1‑93 (MQ=255)

                                                       ggAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTc                                  >  2:1023527/1‑93 (MQ=255)

                                                        gAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTcc                                 <  1:4290156/93‑1 (MQ=255)

                                                        gAACAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTcc                                 <  1:1005282/93‑1 (MQ=255)

                                                         aaCAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCt                                >  1:3893941/1‑93 (MQ=255)

                                                           cAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCg                              >  1:4771215/1‑93 (MQ=255)

                                                           cAAATTGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCg                              <  1:2479874/93‑1 (MQ=255)

                                                               ttGTCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGAcc                          >  2:2802709/1‑93 (MQ=255)

                                                                  tCGGTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACAGCTCGGGAATGCTATTTTGTCCTCGGAccgcc                       <  1:4888845/93‑1 (MQ=255)

                                                                     gTATCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCggg                    <  2:4212553/93‑1 (MQ=255)

                                                                        tCGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGa                 >  1:6028406/1‑93 (MQ=255)

                                                                         cGGCATCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGAt                <  2:3811310/93‑1 (MQ=255)

                                                                             aTCTCCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGaa            >  2:868786/1‑93 (MQ=255)

                                                                              tctcCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGaaa           >  2:2219732/1‑93 (MQ=255)

                                                                               ctcCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAt          <  2:1812883/93‑1 (MQ=255)

                                                                               ctcCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAt          <  2:6693072/93‑1 (MQ=255)

                                                                                tcCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTTGGGAATGCTATTCTGTCCTCGGACCGCCGGGCGATGGAAAtc         >  1:5755993/1‑93 (MQ=255)

                                                                                tcCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtc         <  2:5587594/93‑1 (MQ=255)

                                                                                tcCGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtc         <  1:2391466/93‑1 (MQ=255)

                                                                                 ccGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtcg        >  2:2651708/1‑93 (MQ=255)

                                                                                 ccGCACAGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtcg        <  1:5258467/93‑1 (MQ=255)

                                                                                   gCACAGGAAAAGGTCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtcgtt      <  1:4983452/93‑1 (MQ=255)

                                                                                     acaGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtcgttcg    <  2:5824711/93‑1 (MQ=255)

                                                                                      caGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtcgttcgt   >  2:37095/1‑93 (MQ=255)

                                                                                      caGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAAtcgttcgt   >  2:3175890/1‑93 (MQ=255)

                                                                                       aGGAAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAATCGTtcgtc  >  1:2750529/1‑93 (MQ=255)

GGCAATGCTGCATGGCTGTCATTCGCGCCCTGCTTACTCATTCTGGTGTTAGCGGGGAACAAATTGTCGGTATCGGCATCTCCGCACAGGGAAAGGGCTTGTTTTTGCTGGATAAAAACGACAAACCGCTCGGGAATGCTATTTTGTCCTCGGACCGCCGGGCGATGGAAATCGTTCGTC  >  NC_000913/3744868‑3745047

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