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	97.6%	635.7 / 3.3	206	T408T (ACT→ACC)	yrfF	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (2/3); new base C (94/107); total (96/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

Only 100 of 213 total aligned reads displayed.

AGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATT  >  NC_000913/3525240‑3525419

aGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTc                                                                                         >  1:862339/1‑93 (MQ=255)

aGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTc                                                                                         <  2:1265495/93‑1 (MQ=255)

 gCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTcc                                                                                        <  2:1071108/93‑1 (MQ=255)

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

     ccaGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGaa                                                                                    >  2:4260947/1‑93 (MQ=255)

     ccaGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGaa                                                                                    <  1:4166383/93‑1 (MQ=255)

      caGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAAc                                                                                   <  1:3723988/93‑1 (MQ=255)

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

       aGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAAcc                                                                                  <  1:3786968/93‑1 (MQ=255)

        gCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCt                                                                                 <  2:1539310/93‑1 (MQ=255)

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

          gTAAAACAACTGGCTGACGCTGTCGTGCGGGTAGGCTATGCATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTgg                                                                               <  2:1663770/93‑1 (MQ=255)

            aaaaCAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTATTATTCGAAACTCCGGAACCTGGAg                                                                             >  1:3142426/1‑93 (MQ=255)

            aaaaCAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAg                                                                             >  1:3234704/1‑93 (MQ=255)

              aacaacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAgcg                                                                           >  2:4480232/1‑93 (MQ=255)

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

               acaacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAgcgc                                                                          <  1:3317673/93‑1 (MQ=255)

                caacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGcgcg                                                                         <  2:657802/93‑1 (MQ=255)

                   cTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGaaa                                                                      >  2:467702/1‑93 (MQ=255)

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

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

                      gCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAAcc                                                                   >  1:1871161/1‑93 (MQ=255)

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

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

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

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

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

                         gACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAAt                                                                <  1:4490764/93‑1 (MQ=255)

                         gACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAAt                                                                <  2:4003255/93‑1 (MQ=255)

                          aCGCTGGCGTGCGGGTATGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGTGCGCGAAAACCAAtt                                                               >  2:4211158/1‑93 (MQ=255)

                          aCGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAAtt                                                               <  1:2664368/93‑1 (MQ=255)

                              tGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAATTCAcc                                                           >  2:1561708/1‑93 (MQ=255)

                              tGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCAcc                                                           >  2:2959750/1‑93 (MQ=255)

                               ggCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCt                                                          >  2:311822/1‑93 (MQ=255)

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

                                 ggTGGGGGTAGGCGATACATGTCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCttt                                                        <  1:1379888/92‑1 (MQ=255)

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

                                    gCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCtttttt                                                     <  1:2841110/93‑1 (MQ=255)

                                    gCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCtttttt                                                     <  1:3448961/93‑1 (MQ=255)

                                     cGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAGCCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTa                                                    <  2:1761086/93‑1 (MQ=255)

                                     cGGGTAGGCAATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCttttttt                                                    <  2:2836918/93‑2 (MQ=255)

                                      gggTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTAc                                                   >  1:1565141/1‑93 (MQ=255)

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

                                       ggTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTAcc                                                  <  1:2703526/93‑1 (MQ=255)

                                       ggTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTAcc                                                  <  1:4661599/93‑1 (MQ=255)

                                         tAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCCATTAACCTTTTTTATCGt                                                >  2:3818078/1‑93 (MQ=255)

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

                                          aGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACTTTTTTTACCGtt                                               >  1:1521330/1‑93 (MQ=255)

                                          aGGCGAGACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGtt                                               <  2:409297/93‑1 (MQ=255)

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

                                              gATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGAc                                           <  2:1178798/93‑1 (MQ=255)

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

                                              gATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGAc                                           <  1:1173283/93‑1 (MQ=255)

                                              gATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGAc                                           <  1:1061766/93‑1 (MQ=255)

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

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

                                                tACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCACGAAAACCAATTCACCTTTTTTACCGTTTGACTg                                         >  2:3939917/1‑93 (MQ=255)

                                                 aCATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGc                                        <  1:592266/93‑1 (MQ=255)

                                                  cATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCt                                       >  2:3454904/1‑93 (MQ=255)

                                                   aTTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGACAACCAATTCACCTTTTTTACCGTTTGACTGCTc                                      >  1:3052940/1‑93 (MQ=255)

                                                   aTTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTc                                      <  2:4497442/93‑1 (MQ=255)

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

                                                    ttGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCg                                     >  2:1538570/1‑93 (MQ=255)

                                                    ttGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCg                                     >  1:4419507/1‑93 (MQ=255)

                                                    ttGCGCATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCg                                     >  2:4724289/1‑93 (MQ=255)

                                                       ggTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAg                                  <  1:118292/92‑1 (MQ=255)

                                                       ggTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAg                                  <  1:3881164/92‑1 (MQ=255)

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

                                                          aTTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGatc                               >  1:4291051/1‑93 (MQ=255)

                                                             tgCGGTACGGGAATTTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGatcatc                            <  2:239362/92‑1 (MQ=255)

                                                             agaggTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTTACTGCTCGCAGatcatc                            <  2:66569/90‑1 (MQ=255)

                                                              gCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCt                           <  2:1967531/93‑1 (MQ=255)

                                                                 gTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGa                        <  2:2019828/93‑1 (MQ=255)

                                                                 gTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGa                        >  1:2799754/1‑93 (MQ=255)

                                                                     gggAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAGCCAATGCACCGTTTTTGCCGTTTGACTGCTCGCAGATCATCTGGAacga                    <  1:3919686/93‑1 (MQ=255)

                                                                     gggAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAacga                    >  2:4748515/1‑93 (MQ=255)

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

                                                                        aaTGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGc                 >  1:3992445/1‑93 (MQ=255)

                                                                        aaTGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGc                 >  2:3693391/1‑93 (MQ=255)

                                                                        aaTGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGc                 >  1:721684/1‑93 (MQ=255)

                                                                            tgtAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGc             >  2:2482620/1‑93 (MQ=255)

                                                                            tgtAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGc             >  2:2177295/1‑93 (MQ=255)

                                                                            tgtAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGc             <  1:2291058/93‑1 (MQ=255)

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

                                                                             gtAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCt            <  1:774766/93‑1 (MQ=255)

                                                                              tAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGCGCCCCGCTc           >  1:2635717/1‑93 (MQ=255)

                                                                              tAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCAACTGGAACGACGCCCGCTc           >  1:4541366/1‑93 (MQ=255)

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

                                                                                atatTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCAt         <  1:2553350/93‑1 (MQ=255)

                                                                                atatTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCAt         >  1:2237243/1‑93 (MQ=255)

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

                                                                                  atTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTg       >  2:2082321/1‑93 (MQ=255)

                                                                                   ttCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGc      >  1:3732457/1‑93 (MQ=255)

                                                                                     cGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCa    <  1:1128090/93‑1 (MQ=255)

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

                                                                                       aaCCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCAtt  >  1:1756796/1‑93 (MQ=255)

                                                                                       aaCCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCAtt  >  1:1483682/1‑93 (MQ=255)

                                                                                       aaCCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCAtt  >  2:3327224/1‑93 (MQ=255)

                                                                                       aaCCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCAtt  <  1:3854313/93‑1 (MQ=255)

                                                                                       aaCCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCAtt  >  2:2647561/1‑93 (MQ=255)

AGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATT  >  NC_000913/3525240‑3525419

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