BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,051,230	T→C	V54A (GTA→GCA)	sfa →	Sfa

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,051,230	0	T	C	90.1%	629.4 / 48.7	214	V54A (GTA→GCA)	sfa	Sfa
Reads supporting (aligned to +/- strand): ref base T (7/14); new base C (97/95); total (105/109)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 217 total aligned reads displayed.

TCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGTTAAAAAAACAAGCC  >  NC_000913/1051140‑1051319

tcttcAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCAc                                                                                         <  1:1862584/93‑1 (MQ=255)

 cttcAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACTTGACCAGGGCAcc                                                                                        <  1:6381344/93‑1 (MQ=255)

 cttcAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCAcc                                                                                        <  1:2864845/93‑1 (MQ=255)

  ttcAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCg                                                                                       <  2:710275/93‑1 (MQ=255)

   tcAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCgg                                                                                      <  1:2484739/93‑1 (MQ=255)

    cAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCggg                                                                                     >  1:4476787/1‑93 (MQ=255)

     aaCGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGa                                                                                    >  1:1102736/1‑93 (MQ=255)

     aaCGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGa                                                                                    >  1:4506778/1‑93 (MQ=255)

      aCGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGaa                                                                                   >  1:4597823/1‑93 (MQ=255)

      aCGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGaa                                                                                   >  2:3528023/1‑93 (MQ=255)

        gCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAAcc                                                                                 >  2:5615966/1‑93 (MQ=255)

         cTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCa                                                                                >  2:4262408/1‑93 (MQ=255)

          tgtatgtaTAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCaa                                                                               >  2:668410/1‑93 (MQ=255)

           gtatgtatAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAac                                                                              <  1:6188034/93‑1 (MQ=255)

           gtatgtatAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAac                                                                              >  1:1121797/1‑93 (MQ=255)

           gtatgtatAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAac                                                                              <  1:6135794/93‑1 (MQ=255)

              tgtatAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACAcc                                                                           >  1:6516376/1‑93 (MQ=255)

                tataGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCtt                                                                         <  1:50798/93‑1 (MQ=255)

                 ataGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTa                                                                        <  1:2078675/93‑1 (MQ=255)

                  taGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTAc                                                                       >  1:4843929/1‑93 (MQ=255)

                    gTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTg                                                                     <  2:7417/93‑1 (MQ=255)

                    gTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTg                                                                     >  2:1785164/1‑93 (MQ=255)

                     tAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTgg                                                                    <  1:813310/93‑1 (MQ=255)

                        gCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCgt                                                                 >  2:3349290/1‑93 (MQ=255)

                        gCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCgt                                                                 >  1:5819420/1‑93 (MQ=255)

                         cGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCgtg                                                                <  1:1293170/93‑1 (MQ=255)

                         cGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCgtg                                                                <  1:3775283/93‑1 (MQ=255)

                         cGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCgtg                                                                >  1:2031196/1‑93 (MQ=255)

                          gATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTACTGGCgtgt                                                               <  1:1531776/93‑1 (MQ=255)

                           aTAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGtt                                                              >  2:4753626/1‑93 (MQ=255)

                           aTAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGtt                                                              <  2:6245271/93‑1 (MQ=255)

                            tAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGATCCAACACCTTACTGGCGTGTTg                                                             <  1:5245822/93‑1 (MQ=255)

                            tAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTg                                                             <  2:4951160/93‑1 (MQ=255)

                              aCCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTACTGGCGTGTTGCt                                                           <  2:2425265/93‑1 (MQ=255)

                              aCCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCt                                                           >  2:5092774/1‑93 (MQ=255)

                               ccTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTg                                                          <  2:3370863/93‑1 (MQ=255)

                                cTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGt                                                         >  1:3966700/1‑93 (MQ=255)

                                 tGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTACTGGCGTGTTGCTGTc                                                        >  2:3279796/1‑93 (MQ=255)

                                 tGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTc                                                        >  2:2004670/1‑93 (MQ=255)

                                   ttGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCtt                                                      >  2:2558072/1‑93 (MQ=255)

                                   ttGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCtt                                                      <  2:3470351/93‑1 (MQ=255)

                                     gATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTCCTGGCGTGTTGCTGTCtttt                                                    >  2:1712277/1‑93 (MQ=255)

                                     gATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCtttt                                                    <  2:3268398/93‑1 (MQ=255)

                                       ttattGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTaa                                                  >  2:5221700/1‑93 (MQ=255)

                                       ttattGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTaa                                                  >  2:3676741/1‑93 (MQ=255)

                                        tattGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAg                                                 <  2:6076552/93‑1 (MQ=255)

                                           tGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGAcc                                              <  1:5145058/93‑1 (MQ=255)

                                            gAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCa                                             <  1:4389353/93‑1 (MQ=255)

                                              aTCTTTCGGGGAGATGTCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCaga                                           >  1:6110635/1‑93 (MQ=255)

                                               tCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCagaa                                          <  2:6568576/93‑1 (MQ=255)

                                               tCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCagaa                                          <  1:962500/93‑1 (MQ=255)

                                                cTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCagaag                                         <  1:2572446/93‑1 (MQ=255)

                                                 tttCGGGGAGATGGCTTATAACATTTCTTACCTGACCTGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCagaaga                                        >  2:4658380/1‑93 (MQ=255)

                                                   tCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAgg                                      >  2:5849337/1‑93 (MQ=255)

                                                   tCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAgg                                      <  2:4407994/93‑1 (MQ=255)

                                                     ggggAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGtt                                    >  2:4291704/1‑93 (MQ=255)

                                                      gggAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTa                                   >  1:417740/1‑93 (MQ=255)

                                                      gggAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTa                                   >  2:2992525/1‑93 (MQ=255)

                                                      gggAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTa                                   >  2:4420674/1‑93 (MQ=255)

                                                      gggAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTa                                   <  1:1810410/93‑1 (MQ=255)

                                                       ggAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTaa                                  >  2:4472287/1‑93 (MQ=255)

                                                       ggAGATGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTaa                                  >  2:4144206/1‑93 (MQ=255)

                                                        gagaTGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAAc                                 >  1:3227337/1‑93 (MQ=255)

                                                         agaTGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACa                                <  1:3528814/93‑1 (MQ=255)

                                                          gaTGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAg                               <  1:3074169/93‑1 (MQ=255)

                                                           aTGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGt                              <  1:2261899/93‑1 (MQ=255)

                                                           aTGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGt                              <  1:57516/93‑1 (MQ=255)

                                                            tGGCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTg                             >  2:1188686/1‑93 (MQ=255)

                                                             ggCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGa                            <  1:669253/93‑1 (MQ=255)

                                                              gCTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGaa                           <  1:2011505/93‑1 (MQ=255)

                                                               cTTATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAat                          >  1:3702350/1‑93 (MQ=255)

                                                                ttATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAata                         <  1:2725548/93‑1 (MQ=255)

                                                                ttATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAata                         >  1:2891088/1‑93 (MQ=255)

                                                                ttATAACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAata                         >  2:1105017/1‑93 (MQ=255)

                                                                 tataACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAatat                        <  2:58893/93‑1 (MQ=255)

                                                                  ataACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATAtt                       <  2:3602821/93‑1 (MQ=255)

                                                                  ataACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATAtt                       <  1:2532500/93‑1 (MQ=255)

                                                                   taACATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTg                      <  2:1974439/93‑1 (MQ=255)

                                                                    aaCATTTCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGa                     >  1:4642482/1‑93 (MQ=255)

                                                                        tttCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAagag                 >  1:1611666/1‑93 (MQ=255)

                                                                        tttCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAagag                 >  2:3780673/1‑93 (MQ=255)

                                                                        tttCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAagag                 <  1:1469508/93‑1 (MQ=255)

                                                                        tttCTTACCTGACCAGGGCACCGGGAACCAAAACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGCTAACAGGGAATATTGAagag                 >  1:1685310/1‑93 (MQ=255)

                                                                          tCTTACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGtt               <  2:5996584/93‑1 (MQ=255)

                                                                            ttACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGTTaa             <  1:5453538/93‑1 (MQ=255)

                                                                             tACCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGGTAACAGTGAATATTGAAGAGTTaaa            >  1:5595194/1‑93 (MQ=255)

                                                                              aCCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGTTaaaa           >  2:5015913/1‑93 (MQ=255)

                                                                              aCCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGTTaaaa           <  1:3851843/93‑1 (MQ=255)

                                                                              aCCTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGACTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGTTaaaa           >  1:6320529/1‑93 (MQ=255)

                                                                                cTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGTTaaaaca         >  2:1600605/1‑91 (MQ=255)

                                                                                cTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGTTaaaaaa         >  2:1273678/1‑93 (MQ=255)

                                                                                cTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGTTaaaaaa         >  1:986799/1‑93 (MQ=255)

                                                                                cTGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGTTaaaaaa         >  1:4979997/1‑93 (MQ=255)

                                                                                 tGACCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGTTaaaaaaa        <  1:328825/93‑1 (MQ=255)

                                                                                   aCCAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGTTAAAAAAACa      <  1:4268009/93‑1 (MQ=255)

                                                                                    ccAGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGTTaaaaaaaaaa     >  2:3853396/1‑93 (MQ=255)

                                                                                      aGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGTTAAAAAAACAAGc   >  2:5771511/1‑93 (MQ=255)

                                                                                      aGGGCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGTTAAAAAAACAAGc   <  2:2776871/93‑1 (MQ=255)

                                                                                       gggCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGTTAAAAAAACAAGcc  <  1:2478108/93‑1 (MQ=255)

                                                                                       gggCACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGTTAAAAAAACAAGcc  >  1:1386533/1‑93 (MQ=255)

TCTTCAACGCTGTATGTATAGTAAGCGATAACCTGTTGATTATTGAATCTTTCGGGGAGATGGCTTATAACATTTCTTACCTGACCAGGGTACCGGGAACCAACACCTTACTGGCGTGTTGCTGTCTTTTAAGACCAGAAGAGGTTAACAGTGAATATTGAAGAGTTAAAAAAACAAGCC  >  NC_000913/1051140‑1051319

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