BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	518,231	T→C	K48R (AAA→AGA)	ybbO ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	518,231	0	T	C	87.1%	713.3 / 92.8	295	K48R (AAA→AGA)	ybbO	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (18/20); new base C (115/142); total (133/162)
Fisher's exact test for biased strand distribution p-value = 8.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.27e-02

Only 100 of 296 total aligned reads displayed.

TCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGTAATTAA  >  NC_000913/518141‑518319

 cgGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTcc                                                                                       >  1:7322450/1‑93 (MQ=255)

        gcgATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAAcc                                                                                >  1:5426006/1‑93 (MQ=255)

        gcgATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAAcc                                                                                <  2:5280190/93‑1 (MQ=255)

         cgATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCt                                                                               >  2:5716771/1‑93 (MQ=255)

         cgATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCt                                                                               <  1:6431385/93‑1 (MQ=255)

         cgATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCt                                                                               >  1:2083568/1‑93 (MQ=255)

         cgATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCt                                                                               >  2:7290907/1‑93 (MQ=255)

          gATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTg                                                                              >  2:4026575/1‑93 (MQ=255)

            tCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCCACATCATCCGGTCTCCGGCAACCTGcc                                                                            <  2:5173391/93‑1 (MQ=255)

            tCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGcc                                                                            >  1:1109695/1‑93 (MQ=255)

                cacTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGca                                                                        <  1:3450464/93‑1 (MQ=255)

                 acTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGcac                                                                       >  1:5295959/1‑93 (MQ=255)

                  cTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGcaca                                                                      <  2:6297293/93‑1 (MQ=255)

                  cTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGcaca                                                                      <  2:3608936/93‑1 (MQ=255)

                  cTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGcaca                                                                      >  2:3309102/1‑93 (MQ=255)

                   tttCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACAt                                                                     >  1:2054962/1‑93 (MQ=255)

                   tttCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACAt                                                                     >  2:3811472/1‑93 (MQ=255)

                   tttCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACAt                                                                     <  2:1720379/93‑1 (MQ=255)

                    ttCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATg                                                                    >  2:3211433/1‑93 (MQ=255)

                    ttCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATg                                                                    >  2:2638649/1‑93 (MQ=255)

                     tCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGa                                                                   >  1:1500846/1‑93 (MQ=255)

                     tCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGa                                                                   <  2:7194073/93‑1 (MQ=255)

                      cTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGaa                                                                  <  2:2425495/93‑1 (MQ=255)

                        ggTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGaaaa                                                                <  1:5462139/93‑1 (MQ=255)

                           gAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAAccc                                                             <  1:813359/93‑1 (MQ=255)

                           gAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAAccc                                                             >  2:5072587/1‑93 (MQ=255)

                            aaTCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCt                                                            <  2:501556/93‑1 (MQ=255)

                            aaTCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCt                                                            >  2:7216069/1‑93 (MQ=255)

                            aaTCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCt                                                            <  1:1903469/93‑1 (MQ=255)

                             aTCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTg                                                           <  2:3990545/93‑1 (MQ=255)

                             aTCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTg                                                           <  2:1349162/93‑1 (MQ=255)

                              tCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCATCATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTgg                                                          >  2:1731770/1‑93 (MQ=255)

                               ccAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGc                                                         <  2:6461973/93‑1 (MQ=255)

                                cAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCg                                                        <  2:1012143/93‑1 (MQ=255)

                                   atcgatcAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGttt                                                     <  2:1083791/93‑1 (MQ=255)

                                   atcgatcAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGttt                                                     >  1:5345891/1‑93 (MQ=255)

                                   atcgatcAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGttt                                                     >  2:2221961/1‑93 (MQ=255)

                                      gatcAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTaa                                                  >  2:6036997/1‑93 (MQ=255)

                                       atcAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAAt                                                 >  1:4798290/1‑93 (MQ=255)

                                       atcAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAAt                                                 <  1:6144148/93‑1 (MQ=255)

                                       atcAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAAt                                                 <  2:5929127/93‑1 (MQ=255)

                                       atcAACACGCCGGTAAATCCAATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAAt                                                 <  2:3037494/93‑1 (MQ=255)

                                        tcAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAAtt                                                <  1:6373416/93‑1 (MQ=255)

                                        tcAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAAtt                                                >  1:4536488/1‑93 (MQ=255)

                                        tcAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAAtt                                                <  2:1604179/93‑1 (MQ=255)

                                        tcAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAAtt                                                >  1:2942764/1‑93 (MQ=255)

                                         cAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTc                                               <  2:2001169/93‑1 (MQ=255)

                                           acacGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGa                                             <  1:7503217/93‑1 (MQ=255)

                                           acacGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGa                                             >  2:4437837/1‑93 (MQ=255)

                                            cacGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAg                                            <  2:5501335/93‑1 (MQ=255)

                                            cacGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAg                                            <  1:2593052/93‑1 (MQ=255)

                                             acGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAgc                                           <  1:1158284/93‑1 (MQ=255)

                                               gCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAgcgc                                         <  2:1817549/93‑1 (MQ=255)

                                                 cGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAgcgcgc                                       <  1:2128130/93‑1 (MQ=255)

                                                  ggTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCt                                      >  1:1720995/1‑93 (MQ=255)

                                                   gTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCtt                                     >  2:4266061/1‑93 (MQ=255)

                                                     aaaTCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTc                                   <  1:6646395/93‑1 (MQ=255)

                                                       aTCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCa                                 <  2:2485612/93‑1 (MQ=255)

                                                        tCCCATGCTGTTCATGCGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAg                                <  2:4755140/93‑1 (MQ=255)

                                                        tCCCATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAg                                <  1:2706790/93‑1 (MQ=255)

                                                           cATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGcc                             >  2:6568359/1‑93 (MQ=255)

                                                           cATGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGcc                             <  1:3755000/93‑1 (MQ=255)

                                                             tGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCaa                           <  2:6994167/93‑1 (MQ=255)

                                                             tGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCaa                           <  2:2144714/93‑1 (MQ=255)

                                                             tGCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCaa                           >  2:3494008/1‑93 (MQ=255)

                                                              gCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCATCCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAAt                          <  2:4578137/93‑1 (MQ=255)

                                                              gCTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAAt                          <  2:3176011/93‑1 (MQ=255)

                                                               cTGTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAAtt                         >  2:63114/1‑93 (MQ=255)

                                                                tGTTCATGCGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTc                        >  1:512491/1‑93 (MQ=255)

                                                                 gTTCATGCGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTcc                       >  1:1650709/1‑93 (MQ=255)

                                                                 gTTCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTcc                       >  1:1399955/1‑93 (MQ=255)

                                                                  ttCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCa                      <  1:4204635/93‑1 (MQ=255)

                                                                  ttCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCa                      <  1:3099404/93‑1 (MQ=255)

                                                                  ttCATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCa                      >  2:7044148/1‑93 (MQ=255)

                                                                    cATGCGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACt                    <  2:5544578/93‑1 (MQ=255)

                                                                    cATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACt                    <  1:4323120/93‑1 (MQ=255)

                                                                    cATGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACt                    <  1:6950963/93‑1 (MQ=255)

                                                                     aTGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTg                   <  1:2345591/93‑1 (MQ=255)

                                                                     aTGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTg                   >  2:1883353/1‑93 (MQ=255)

                                                                     aTGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTg                   <  1:5611168/93‑1 (MQ=255)

                                                                     aTGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTg                   <  2:3391797/93‑1 (MQ=255)

                                                                      tGCGCTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTgg                  <  2:7509505/93‑1 (MQ=255)

                                                                         gcTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAAc               <  1:3274104/93‑1 (MQ=255)

                                                                         gcTCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAAc               <  1:2857702/93‑1 (MQ=255)

                                                                           tCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACAt             <  2:2291041/93‑1 (MQ=255)

                                                                           tCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACAt             >  2:4352965/1‑93 (MQ=255)

                                                                           tCAACATCATCCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACAt             >  1:6865277/1‑93 (MQ=255)

                                                                                atcatcCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGt        <  1:824760/93‑1 (MQ=255)

                                                                                atcatcCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGt        >  2:1276957/1‑93 (MQ=255)

                                                                                atcatcCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGt        <  2:6691033/93‑1 (MQ=255)

                                                                                 tcatcCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGTa       >  2:789736/1‑93 (MQ=255)

                                                                                 tcatcCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGTa       >  1:7163095/1‑93 (MQ=255)

                                                                                  catcCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGTaa      <  1:4673347/93‑1 (MQ=255)

                                                                                   atcCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGTAAt     <  1:5836411/93‑1 (MQ=255)

                                                                                    tcCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGTAAtt    <  1:505826/93‑1 (MQ=255)

                                                                                    tcCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGTAAtt    >  2:2856303/1‑93 (MQ=255)

                                                                                    tcCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGTAAtt    >  2:1406776/1‑93 (MQ=255)

                                                                                    tcCGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGTAAtt    <  2:4861264/93‑1 (MQ=255)

                                                                                     ccGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGTAATTa   <  2:5387667/93‑1 (MQ=255)

                                                                                     ccGGTCTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGTAATTa   <  1:4859526/93‑1 (MQ=255)

TCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGTAATTAA  >  NC_000913/518141‑518319

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