BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	866,644	A→G	N285S (AAT→AGT)	ybiK →	putative asparaginase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	866,644	0	A	G	98.8%	756.5 / ‑1.4	244	N285S (AAT→AGT)	ybiK	putative asparaginase
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (130/110); total (131/113)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

Only 100 of 259 total aligned reads displayed.

AGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAG  >  NC_000913/866553‑866733

 gTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTg                                                                                         <  2:9719679/93‑1 (MQ=255)

    tcGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTcg                                                                                      <  1:428920/93‑1 (MQ=255)

     cgcgGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTcgc                                                                                     <  2:1066326/93‑1 (MQ=255)

     cgcgGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTcgc                                                                                     <  1:7431176/93‑1 (MQ=255)

       cgGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCgcgc                                                                                   >  2:3829174/1‑93 (MQ=255)

        ggAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCt                                                                                  >  1:2723835/1‑93 (MQ=255)

         gAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTa                                                                                 <  1:1403912/93‑1 (MQ=255)

           aGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTAcc                                                                               <  1:1134337/93‑1 (MQ=255)

           aGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTAcc                                                                               <  1:3364531/93‑1 (MQ=255)

           aGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTAcc                                                                               <  1:2140381/93‑1 (MQ=255)

           aGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGAGTCGCGCTAcc                                                                               >  2:1729972/1‑93 (MQ=255)

            gCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCg                                                                              <  2:9388988/93‑1 (MQ=255)

               tGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGttt                                                                           <  2:2899221/93‑1 (MQ=255)

               tGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGttt                                                                           <  2:7581135/93‑1 (MQ=255)

                gCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTa                                                                          >  1:7152807/1‑93 (MQ=255)

                gCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTa                                                                          >  2:935424/1‑93 (MQ=255)

                gCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTa                                                                          >  2:5314118/1‑93 (MQ=255)

                 cGAGCGGGTAGTAATGGAAAACCTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTaa                                                                         <  1:4955127/93‑1 (MQ=255)

                   aGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAaca                                                                       <  1:10010513/93‑1 (MQ=255)

                   aGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAaca                                                                       >  1:6446814/1‑93 (MQ=255)

                   aGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGACGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAaca                                                                       <  1:9669934/93‑1 (MQ=255)

                    gCGGGTAGTAATAGAGAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAacac                                                                      <  2:8953674/93‑1 (MQ=255)

                      gggTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCg                                                                    >  1:2032716/1‑93 (MQ=255)

                      gggTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCg                                                                    >  2:4239352/1‑93 (MQ=255)

                       ggTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGa                                                                   <  2:10033262/93‑1 (MQ=255)

                        gtagtaATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGaa                                                                  >  1:6431241/1‑93 (MQ=255)

                        gtagtaATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGaa                                                                  >  2:9696732/1‑93 (MQ=255)

                            taATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGaa                                                              >  1:5295550/1‑93 (MQ=255)

                               tGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGt                                                           <  2:8353222/93‑1 (MQ=255)

                               tGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGt                                                           >  1:3842752/1‑93 (MQ=255)

                                ggAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTa                                                          <  1:5164593/93‑1 (MQ=255)

                                ggAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTa                                                          <  2:1261176/93‑1 (MQ=255)

                                 gAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGGTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTAt                                                         <  1:1817849/93‑1 (MQ=255)

                                 gAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTAt                                                         >  1:6825210/1‑93 (MQ=255)

                                 gAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAAGGTAt                                                         >  2:1293492/1‑93 (MQ=255)

                                  aaaaaCTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATc                                                        <  2:3239255/93‑1 (MQ=255)

                                  aaaaaCTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATc                                                        <  1:957054/93‑1 (MQ=255)

                                    aaaCTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATcgc                                                      >  1:3251908/1‑93 (MQ=255)

                                    aaaCTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATcgc                                                      >  1:8753164/1‑93 (MQ=255)

                                        tCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCt                                                  <  2:7766218/93‑1 (MQ=255)

                                        tCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCt                                                  <  2:3106676/93‑1 (MQ=255)

                                        tCCCTGCGCTTGGCGGTAGCGGTGGCTTAATAGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATGGCGccg                                                  >  2:5373508/1‑92 (MQ=255)

                                         cccTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGGGCCTg                                                 >  2:4050540/1‑93 (MQ=255)

                                         cccTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTg                                                 <  2:3724739/93‑1 (MQ=255)

                                         cccTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTAGCGCTACCGTTTAACACCGAAGGAATGTATCACGCCTg                                                 >  2:5691828/1‑93 (MQ=255)

                                          ccTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTgg                                                >  1:7920629/1‑93 (MQ=255)

                                          ccTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTgg                                                >  1:8588881/1‑93 (MQ=255)

                                          ccTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTgg                                                <  1:3558109/93‑1 (MQ=255)

                                             gcgcTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGc                                             <  2:5299718/93‑1 (MQ=255)

                                              ggcTGGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGGCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCt                                            <  1:3959018/92‑1 (MQ=255)

                                               gcTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTa                                           <  1:606785/93‑1 (MQ=255)

                                                cTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTAc                                          <  1:9330251/93‑1 (MQ=255)

                                                 ttGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACg                                         <  2:7893143/93‑1 (MQ=255)

                                                  tGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGc                                        >  2:4393494/1‑93 (MQ=255)

                                                  tGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGc                                        >  1:5755582/1‑93 (MQ=255)

                                                  tGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGc                                        >  2:6312445/1‑93 (MQ=255)

                                                  tGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGc                                        >  1:6582067/1‑93 (MQ=255)

                                                     cGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAgg                                     >  2:6518734/1‑93 (MQ=255)

                                                      ggTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGc                                    <  2:6780314/93‑1 (MQ=255)

                                                      ggTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGc                                    >  2:2509075/1‑93 (MQ=255)

                                                       gTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAggag                                   >  2:7065051/1‑91 (MQ=255)

                                                       gTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCg                                   >  2:5501684/1‑93 (MQ=255)

                                                       gTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGCGCTACGCAGGCg                                   >  2:3463697/1‑93 (MQ=255)

                                                        gaGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGa                                  <  1:8092736/92‑1 (MQ=255)

                                                          gCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATa                                <  2:9013039/93‑1 (MQ=255)

                                                          gCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATa                                <  1:5027693/93‑1 (MQ=255)

                                                          gCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGC‑GGCGATAc                               >  2:6900814/1‑93 (MQ=255)

                                                             gTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGc                             >  1:9043034/1‑93 (MQ=255)

                                                             gTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGc                             >  1:7841390/1‑93 (MQ=255)

                                                               ggCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCa                           <  2:6184660/93‑1 (MQ=255)

                                                               ggCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCa                           >  1:7344836/1‑93 (MQ=255)

                                                               ggCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCa                           >  1:6143293/1‑93 (MQ=255)

                                                                gCTTAATCTCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCaa                          <  1:4544840/93‑1 (MQ=255)

                                                                gCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCaa                          >  2:1974169/1‑93 (MQ=255)

                                                                 cTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAac                         >  2:4639163/1‑93 (MQ=255)

                                                                 cTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAac                         >  1:5519199/1‑93 (MQ=255)

                                                                 cTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAac                         <  1:4652748/93‑1 (MQ=255)

                                                                  ttAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAacc                        >  1:7495397/1‑93 (MQ=255)

                                                                   tAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAaccc                       >  2:2022051/1‑92 (MQ=255)

                                                                   tAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAacca                       >  1:4161643/1‑93 (MQ=255)

                                                                    aaTCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAaccac                      >  1:4840926/1‑93 (MQ=255)

                                                                    aaTCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAaccac                      <  2:5634447/93‑1 (MQ=255)

                                                                     aTCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAaccacc                     >  2:8090238/1‑93 (MQ=255)

                                                                      tCGCTATCGACGATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCg                    <  2:1898847/93‑1 (MQ=255)

                                                                         cTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTa                 >  2:1360411/1‑93 (MQ=255)

                                                                          tATCGACCATGAAGGGAGTGTCGCGCTACCTTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTAt                <  2:7334879/93‑1 (MQ=255)

                                                                          tATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTAt                >  2:4880594/1‑93 (MQ=255)

                                                                                ccATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCg          >  2:8494003/1‑93 (MQ=255)

                                                                                  aTGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTg        >  1:6971323/1‑93 (MQ=255)

                                                                                  aTGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTg        >  1:1964199/1‑93 (MQ=255)

                                                                                   tGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGa       >  1:8481332/1‑93 (MQ=255)

                                                                                   tGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGa       <  2:8825945/93‑1 (MQ=255)

                                                                                    gAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGaa      >  1:1995093/1‑93 (MQ=255)

                                                                                    gAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGaa      >  1:2247707/1‑93 (MQ=255)

                                                                                      aGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGaaaa    <  1:2658657/93‑1 (MQ=255)

                                                                                      aGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGaaaa    <  2:9358843/93‑1 (MQ=255)

                                                                                      aGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGaaaa    <  1:7160797/93‑1 (MQ=255)

                                                                                       gggAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGaaaaa   >  1:8299109/1‑93 (MQ=255)

                                                                                       gggAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGaaaaa   <  2:4294611/93‑1 (MQ=255)

                                                                                        ggAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGACAAAg  >  1:9929403/1‑93 (MQ=255)

AGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAG  >  NC_000913/866553‑866733

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