BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,414,329	C→T	I521I (ATC→ATT)	pta →	phosphotransacetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,414,329	0	C	T	99.5%	778.7 / ‑5.1	207	I521I (ATC→ATT)	pta	phosphotransacetylase
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (113/93); total (114/93)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.19e-01

Only 100 of 239 total aligned reads displayed.

AGACAACGTGGTGCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATCCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGT  >  NC_000913/2414239‑2414418

aGACAACGTGGTGCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCAttcg                                                                                         <  1:4704059/93‑1 (MQ=255)

aGACAACGTGGTGCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCAttcg                                                                                         >  1:6110875/1‑93 (MQ=255)

aGACAACGTGGTGCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCAttcg                                                                                         >  2:5518000/1‑93 (MQ=255)

 gACAACGTGGTGCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCAttcgt                                                                                        <  2:1958259/93‑1 (MQ=255)

  aCAACGTGGTGCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCAttcgtc                                                                                       <  2:2891736/93‑1 (MQ=255)

  aCAACGTGGTGCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCAttcgtc                                                                                       <  1:3610037/93‑1 (MQ=255)

   cAACGTGGTGCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCAttcgtcc                                                                                      >  1:436519/1‑93 (MQ=255)

   cAACGTGGTGCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCAttcgtcc                                                                                      >  2:3498162/1‑93 (MQ=255)

   cAACGTGGTGCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCAttcgtcc                                                                                      <  2:5696625/93‑1 (MQ=255)

     aCGTGGTGCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTccgc                                                                                    <  2:148462/93‑1 (MQ=255)

     aCGTGGTGCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTccgc                                                                                    >  1:4983267/1‑93 (MQ=255)

     aCGTGGTGCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTccgc                                                                                    <  2:4257667/93‑1 (MQ=255)

      cGTGGTGCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTccgcc                                                                                   >  2:4152901/1‑93 (MQ=255)

        tggtgCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCcgccgc                                                                                 <  1:6384085/93‑1 (MQ=255)

        tggtgCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCcgccgc                                                                                 <  1:3966763/93‑1 (MQ=255)

            gCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCa                                                                             >  2:2477108/1‑93 (MQ=255)

              tCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGc                                                                           >  1:2738111/1‑93 (MQ=255)

                  tACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGAt                                                                       <  2:4047293/93‑1 (MQ=255)

                    cGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCa                                                                     <  1:6104244/93‑1 (MQ=255)

                    cGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCa                                                                     <  2:6130364/93‑1 (MQ=255)

                     gCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCaa                                                                    >  1:697033/1‑93 (MQ=255)

                      cTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCaaa                                                                   >  1:5912967/1‑93 (MQ=255)

                      cTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCaaa                                                                   >  1:1594727/1‑93 (MQ=255)

                      cTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCaaa                                                                   >  1:15129/1‑93 (MQ=255)

                       tGATGCTGGAACAGGATGAAGTTGATGGTCTGTTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCaaaa                                                                  <  2:2701837/93‑1 (MQ=255)

                       tGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCaaaa                                                                  >  1:6010218/1‑93 (MQ=255)

                        gATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAAc                                                                 >  2:2611183/1‑93 (MQ=255)

                        gATGCTGGAACAGGATGAAGTTGAGGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAAc                                                                 <  1:1894859/93‑1 (MQ=255)

                         aTGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACt                                                                <  1:5370720/93‑1 (MQ=255)

                            cTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCa                                                             <  2:6458731/93‑1 (MQ=255)

                             tGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAGAACTGCAc                                                            <  1:6325755/93‑1 (MQ=255)

                              ggAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCAcc                                                           >  1:5297607/1‑93 (MQ=255)

                               gAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCg                                                          >  2:113526/1‑93 (MQ=255)

                               gAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCg                                                          >  2:5192474/1‑93 (MQ=255)

                                aaCAGGATGAAGTTGATGGTCTGGTTACCGGTGCTGTTCACACTACCGCAAACACTATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCgg                                                         <  1:886449/93‑1 (MQ=255)

                                 aCAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCggg                                                        >  1:2584801/1‑93 (MQ=255)

                                 aCAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCggg                                                        >  2:233360/1‑93 (MQ=255)

                                   aGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCa                                                      >  1:5151955/1‑93 (MQ=255)

                                    ggATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAg                                                     >  2:5281166/1‑93 (MQ=255)

                                    ggATGAAGTTGATGGTCTGGTATCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAg                                                     >  1:4093448/1‑93 (MQ=255)

                                     gATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGc                                                    <  1:203180/93‑1 (MQ=255)

                                      aTGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCt                                                   <  1:2797088/93‑1 (MQ=255)

                                        gAAGTTGATGGTCTGGTTTCCGGTCCTGTTCACGCTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTcc                                                 <  1:2274656/93‑1 (MQ=255)

                                         aaGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTccc                                                <  2:4380791/93‑1 (MQ=255)

                                         aaGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTccc                                                >  2:5078432/1‑93 (MQ=255)

                                         aaGTGGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTccc                                                <  1:2889561/93‑1 (MQ=255)

                                          aGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCt                                               <  2:1525785/93‑1 (MQ=255)

                                          aGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCt                                               <  1:849841/93‑1 (MQ=255)

                                          aGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCt                                               <  1:3662417/93‑1 (MQ=255)

                                           gTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTg                                              <  2:3288805/93‑1 (MQ=255)

                                              gATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTggca                                           >  1:1916260/1‑91 (MQ=255)

                                                tGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATc                                         <  2:5358449/93‑1 (MQ=255)

                                                     tGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCg                                    >  2:2677319/1‑93 (MQ=255)

                                                       gTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCgtg                                  <  2:1816587/93‑1 (MQ=255)

                                                       gTTTCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCgtg                                  >  1:4590489/1‑93 (MQ=255)

                                                        tttCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCgtgt                                 >  2:5584692/1‑93 (MQ=255)

                                                          tCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGttc                               >  1:5580939/1‑93 (MQ=255)

                                                          tCCGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGttc                               <  2:3037762/93‑1 (MQ=255)

                                                           ccGGTGCTTTTCACACTACCGCAAACGCCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGttct                              <  2:2655428/93‑1 (MQ=255)

                                                            cGGTGCTTTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGttctt                             <  2:4357597/93‑1 (MQ=255)

                                                            cGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGttctt                             <  2:488318/93‑1 (MQ=255)

                                                            cGGTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCCGGCAGCTCCCTGGTATCTTCCGTGttctt                             >  1:1302017/1‑93 (MQ=255)

                                                             ggTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGttcttc                            >  2:1400154/1‑93 (MQ=255)

                                                              gTGCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCa                           >  2:4444654/1‑93 (MQ=255)

                                                               tGCTGTTCACACTACCGCAAACCCCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAt                          <  1:4905348/93‑1 (MQ=255)

                                                               ggCTGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAt                          <  2:4862691/92‑1 (MQ=255)

                                                                 gtGTTCACACCACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgc                        <  2:2756751/92‑1 (MQ=255)

                                                                 cTGTTCCCACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgc                        <  2:5747116/93‑1 (MQ=255)

                                                                 cTGTTCACACTACCGCAAACACCATGCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgc                        <  1:3205500/93‑1 (MQ=255)

                                                                 cTGTCCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgc                        <  1:3051162/93‑1 (MQ=255)

                                                                  tGTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgct                       >  1:4748986/1‑93 (MQ=255)

                                                                   gTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgctg                      >  2:2806538/1‑93 (MQ=255)

                                                                   gTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgctg                      >  2:3644164/1‑93 (MQ=255)

                                                                   gTTCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgctg                      >  1:1493970/1‑93 (MQ=255)

                                                                    ttCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgctgc                     <  1:4803555/93‑1 (MQ=255)

                                                                    ttCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgctgc                     >  2:5416695/1‑93 (MQ=255)

                                                                    ttCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgctgc                     >  2:3940235/1‑93 (MQ=255)

                                                                    ttCACACTACCGCAAACACAATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgctgc                     <  1:2464126/93‑1 (MQ=255)

                                                                     tCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgctgct                    <  1:5641438/93‑1 (MQ=255)

                                                                     tCACACTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgctgct                    >  2:3897010/1‑93 (MQ=255)

                                                                      cacaCTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgctgctg                   >  1:1133522/1‑93 (MQ=255)

                                                                       acacTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgctgctgc                  <  2:6242841/93‑1 (MQ=255)

                                                                        cacTACCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGcc                 <  1:6237977/93‑1 (MQ=255)

                                                                            aCCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGaa             >  2:3008283/1‑93 (MQ=255)

                                                                            aCCGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTACTGCCGGaa             >  2:6414181/1‑93 (MQ=255)

                                                                             ccGCAAACACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAAc            <  1:2666722/93‑1 (MQ=255)

                                                                              cccAAACCCCATTCGGCCGCCCGTGCAGCCGATCAAAACCGCCCCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACa           <  2:1022953/91‑1 (MQ=255)

                                                                               ccAACCACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAg          <  1:6051002/92‑1 (MQ=255)

                                                                                 aaaCCCCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGgt        <  1:21522/93‑1 (MQ=255)

                                                                                 aaaCACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGgt        >  1:1035984/1‑93 (MQ=255)

                                                                                 aaaCACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGgt        >  1:3499456/1‑93 (MQ=255)

                                                                                  aaCACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGgtt       <  2:5549645/93‑1 (MQ=255)

                                                                                  aaCACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGgtt       >  2:5986432/1‑93 (MQ=255)

                                                                                  aaCACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGgtt       <  1:3317659/93‑1 (MQ=255)

                                                                                  aaCACCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGgtt       >  2:6440174/1‑93 (MQ=255)

                                                                                   acacCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGgttt      >  2:1597492/1‑93 (MQ=255)

                                                                                   acacCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGgttt      <  1:421492/93‑1 (MQ=255)

                                                                                   acacCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGgttt      <  1:2407091/93‑1 (MQ=255)

                                                                                   acacCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGgttt      <  2:873378/93‑1 (MQ=255)

                                                                                    cacCATTCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGgttta     <  2:4568675/93‑1 (MQ=255)

AGACAACGTGGTGCTCGGTACGCTGATGCTGGAACAGGATGAAGTTGATGGTCTGGTTTCCGGTGCTGTTCACACTACCGCAAACACCATCCGTCCGCCGCTGCAGCTGATCAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGT  >  NC_000913/2414239‑2414418

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