BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,525,330	T→C	T408T (ACT→ACC)	yrfF →	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,525,330	0	T	C	99.6%	739.1 / ‑5.0	238	T408T (ACT→ACC)	yrfF	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (110/127); total (111/127)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

Only 100 of 247 total aligned reads displayed.

AGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATT  >  NC_000913/3525240‑3525419

    accaGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGa                                                                                     >  1:407438/1‑93 (MQ=255)

       aGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAAcc                                                                                  >  1:292579/1‑93 (MQ=255)

         cGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTg                                                                                <  2:1566914/93‑1 (MQ=255)

          gTAAAACAACTGGCTGACGCTGTCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTgg                                                                               <  1:4545835/93‑1 (MQ=255)

           tAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGa                                                                              <  2:4560270/93‑1 (MQ=255)

           tAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGa                                                                              <  2:2515536/93‑1 (MQ=255)

             aaaCAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAgc                                                                            <  1:1710319/93‑1 (MQ=255)

               acaacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAgcgc                                                                          >  1:3277833/1‑93 (MQ=255)

                 agcTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGa                                                                        <  1:1739787/91‑1 (MQ=255)

                 aacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGa                                                                        >  1:4604380/1‑93 (MQ=255)

                 aacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGa                                                                        >  1:2218199/1‑93 (MQ=255)

                 aacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAAACTCCGGAACCTGGAGCGCGa                                                                        >  1:2321387/1‑93 (MQ=255)

                  acTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGaa                                                                       <  2:2423781/93‑1 (MQ=255)

                    tGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGaaaa                                                                     <  2:466240/93‑1 (MQ=255)

                    tGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGaaaa                                                                     <  2:3813189/93‑1 (MQ=255)

                     ggCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAAc                                                                    >  2:4094987/1‑93 (MQ=255)

                     ggCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAAc                                                                    <  1:4478631/93‑1 (MQ=255)

                        tGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCaa                                                                 >  1:2087430/1‑93 (MQ=255)

                         gACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAAt                                                                >  2:3737424/1‑93 (MQ=255)

                          aCGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAAtt                                                               >  1:2958920/1‑93 (MQ=255)

                             cTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCAc                                                            >  1:1428442/1‑93 (MQ=255)

                             cTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCAc                                                            >  2:247351/1‑93 (MQ=255)

                              tGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGACCCTCCGGAACCTGGAGCGCGAAAACCAATTCAcc                                                           <  2:2968562/93‑1 (MQ=255)

                                gCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCtt                                                         <  2:3906018/93‑1 (MQ=255)

                                 cGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCttt                                                        <  2:3064518/93‑1 (MQ=255)

                                  gggCGGGTAGGCCATACATTGCGTAGTAGCGGTACGGGAATGTGTAATATTCGACCCTCGGGAACCTGGAGCGCGAAAACCAATTCCCCtttt                                                       <  2:3272608/91‑1 (MQ=255)

                                  gTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCtttt                                                       <  2:1178989/93‑1 (MQ=255)

                                   tGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGGGCGAAAACCACTTCACCttttt                                                      >  1:2604126/1‑93 (MQ=255)

                                   tGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCttttt                                                      <  1:495790/93‑1 (MQ=255)

                                   tGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCttttt                                                      <  2:2013997/93‑1 (MQ=255)

                                   tGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCttttt                                                      <  1:1049847/93‑1 (MQ=255)

                                   tGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCttttt                                                      <  1:2040794/93‑1 (MQ=255)

                                     cGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTa                                                    <  1:3483380/93‑1 (MQ=255)

                                      gggTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTAc                                                   <  1:3553362/93‑1 (MQ=255)

                                      gggTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTAc                                                   >  2:2355301/1‑93 (MQ=255)

                                       ggTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTAcc                                                  >  1:4247200/1‑93 (MQ=255)

                                         tAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGt                                                >  2:2862087/1‑93 (MQ=255)

                                          aGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGtt                                               <  2:4486754/93‑1 (MQ=255)

                                          aGGCGATACATTGCGTATTAGCGGTACGGGAATGAGCAATATTCGAACCTCCGGACCCTGGCGCGCGAAAACCAATTCACCTTTTTTAGCGtt                                               >  2:1618916/1‑93 (MQ=255)

                                           ggCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGttt                                              <  2:1483177/93‑1 (MQ=255)

                                           ggCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGttt                                              <  1:999777/93‑1 (MQ=255)

                                           ggCGATACATTGCGTAGTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGttt                                              <  1:2372324/93‑1 (MQ=255)

                                             cGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGa                                            <  2:4214560/93‑1 (MQ=255)

                                              gATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGATCCTCCGGAACCTGGAGCGCGAAAACCAATTCCCCTTTTTTACCGTTTGAc                                           <  2:3635745/93‑1 (MQ=255)

                                              gATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGAc                                           >  2:3947150/1‑93 (MQ=255)

                                               aTACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACt                                          <  2:2758298/93‑1 (MQ=255)

                                               aTACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACt                                          >  2:2324805/1‑93 (MQ=255)

                                               aTACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACt                                          >  2:1233181/1‑93 (MQ=255)

                                               aTACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACt                                          >  2:2553750/1‑93 (MQ=255)

                                                tACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTg                                         <  2:2745696/93‑1 (MQ=255)

                                                tACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTg                                         >  1:3462382/1‑93 (MQ=255)

                                                tACATTGCGTATTAGCGGTACGGGAATCTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTg                                         >  2:913027/1‑93 (MQ=255)

                                                   aTTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTc                                      >  2:4350973/1‑93 (MQ=255)

                                                   aTTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTc                                      >  2:2910645/1‑93 (MQ=255)

                                                   aTTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTc                                      >  2:1836736/1‑93 (MQ=255)

                                                   aTTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTCTTTACCGTTTGACTGCTc                                      >  2:1961111/1‑93 (MQ=255)

                                                      gCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCa                                   >  2:1288928/1‑93 (MQ=255)

                                                       cGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAg                                  >  1:1365441/1‑93 (MQ=255)

                                                        gTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGa                                 >  1:2812768/1‑93 (MQ=255)

                                                        gTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGa                                 >  1:3839024/1‑93 (MQ=255)

                                                         tATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGat                                <  2:1916075/93‑1 (MQ=255)

                                                         tATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGat                                >  2:1390589/1‑93 (MQ=255)

                                                         tATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGat                                >  1:755149/1‑93 (MQ=255)

                                                          aTTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGatc                               <  1:2143903/93‑1 (MQ=255)

                                                          aTTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGatc                               <  2:3354922/93‑1 (MQ=255)

                                                           ttcgCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGatca                              <  2:3290105/90‑1 (MQ=255)

                                                           ttAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGatca                              >  1:213342/1‑93 (MQ=255)

                                                           ttAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGatca                              <  1:1289368/93‑1 (MQ=255)

                                                           ttAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGatca                              <  2:2116739/93‑1 (MQ=255)

                                                            tAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGGACCTGGAGCGCGAAAACCAATTCAAATTTTTTACCGTTTGACTGCTCGCAGatcat                             >  2:3155925/1‑93 (MQ=255)

                                                            tAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGatcat                             >  1:4350028/1‑93 (MQ=255)

                                                            tAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGatcat                             >  1:4090249/1‑93 (MQ=255)

                                                              gCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCt                           >  1:623308/1‑93 (MQ=255)

                                                              gCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCt                           <  2:626119/93‑1 (MQ=255)

                                                              gCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCt                           <  1:2100916/93‑1 (MQ=255)

                                                              gCGGTACGGGAATGTGTAATATTCGAACCCCCGGAACCTGGAGCCAGAACACCAATTCACCTTTTTTACCGCTTGACTGCTCGCAGATCATCt                           >  2:690571/1‑93 (MQ=255)

                                                               cGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTg                          <  1:1886676/93‑1 (MQ=255)

                                                                  tACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGaa                       >  2:501590/1‑93 (MQ=255)

                                                                  tACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGaa                       <  2:4433840/93‑1 (MQ=255)

                                                                   aCGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAac                      <  1:552885/93‑1 (MQ=255)

                                                                   aCGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAac                      <  2:188897/93‑1 (MQ=255)

                                                                   aCGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAac                      >  2:2083763/1‑93 (MQ=255)

                                                                    cGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAacg                     <  1:929667/93‑1 (MQ=255)

                                                                     gggAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAacga                    <  1:1440447/93‑1 (MQ=255)

                                                                      ggAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAacgac                   <  2:541700/93‑1 (MQ=255)

                                                                       gAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAacgacg                  <  2:1077479/93‑1 (MQ=255)

                                                                        aaTGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTATCGTTTGACTGCTCGCAGATCATCTGGAACGACGc                 >  2:1684216/1‑93 (MQ=255)

                                                                         aTGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGcc                >  1:834816/1‑93 (MQ=255)

                                                                           gtgtAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCg              <  1:1610099/93‑1 (MQ=255)

                                                                             gtAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCt            >  1:646542/1‑93 (MQ=255)

                                                                               aaTATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCa          >  1:3422560/1‑93 (MQ=255)

                                                                               aaTATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGAAGCCCGCTCa          >  2:83701/1‑93 (MQ=255)

                                                                                 tatTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCAtt        <  1:2618953/93‑1 (MQ=255)

                                                                                   ttCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGc      <  1:1008064/93‑1 (MQ=255)

                                                                                   ttCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGc      <  1:1165994/93‑1 (MQ=255)

                                                                                     cGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCa    >  2:43916/1‑93 (MQ=255)

                                                                                      gAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCAt   >  2:1262236/1‑93 (MQ=255)

                                                                                      gAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCAt   >  2:2230307/1‑93 (MQ=255)

                                                                                      gAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCAt   <  1:1521525/93‑1 (MQ=255)

                                                                                       aaCCTCCGGAATCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCAtt  <  1:287858/93‑1 (MQ=255)

AGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATT  >  NC_000913/3525240‑3525419

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: `‑`