BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	612,751	T→C	A238A (GCT→GCC)	fes →	enterochelin esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	612,751	0	T	C	87.3%	701.7 / 88.6	283	A238A (GCT→GCC)	fes	enterochelin esterase
Reads supporting (aligned to +/- strand): ref base T (24/12); new base C (121/126); total (145/138)
Fisher's exact test for biased strand distribution p-value = 5.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

Only 100 of 283 total aligned reads displayed.

TATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACT  >  NC_000913/612661‑612839

tATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCAt                                                                                        <  1:4981407/93‑1 (MQ=255)

 aTCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCAtt                                                                                       >  1:1543626/1‑93 (MQ=255)

  tCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTg                                                                                      >  2:3413683/1‑93 (MQ=255)

  tCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTg                                                                                      >  2:7151027/1‑93 (MQ=255)

  tCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTg                                                                                      <  2:5465153/93‑1 (MQ=255)

    gACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGcc                                                                                    >  2:1557084/1‑93 (MQ=255)

    gACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGcc                                                                                    <  2:4701739/93‑1 (MQ=255)

     acacCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGccc                                                                                   >  1:6563446/1‑93 (MQ=255)

     acacCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGccc                                                                                   <  2:1458620/93‑1 (MQ=255)

     acacCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGccc                                                                                   <  2:5153115/93‑1 (MQ=255)

      caccacGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGcccc                                                                                  <  2:131945/93‑1 (MQ=255)

      caccacGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGcccc                                                                                  >  2:3075699/1‑93 (MQ=255)

       accacGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCt                                                                                 >  2:2720457/1‑93 (MQ=255)

        ccacGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCtt                                                                                >  2:2291940/1‑93 (MQ=255)

          acGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCtttt                                                                              >  2:3298057/1‑93 (MQ=255)

             cACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGc                                                                           <  1:7306274/93‑1 (MQ=255)

              aCCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCg                                                                          >  1:1119034/1‑93 (MQ=255)

              aCCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCg                                                                          >  1:7170335/1‑93 (MQ=255)

              aCCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCg                                                                          >  1:2150343/1‑93 (MQ=255)

               ccGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGa                                                                         <  1:2929256/93‑1 (MQ=255)

               ccGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGa                                                                         <  2:4307847/93‑1 (MQ=255)

               ccGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGa                                                                         >  1:5235624/1‑93 (MQ=255)

                cgcgCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGAt                                                                        >  1:533600/1‑93 (MQ=255)

                  cgcCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCg                                                                      <  1:6252968/93‑1 (MQ=255)

                  cgcCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCg                                                                      >  2:5282570/1‑93 (MQ=255)

                   gcCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGt                                                                     <  2:2158075/93‑1 (MQ=255)

                   gcCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGt                                                                     <  2:7297035/93‑1 (MQ=255)

                   gcCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGt                                                                     <  2:1115483/93‑1 (MQ=255)

                    cccACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTg                                                                    >  1:3375211/1‑93 (MQ=255)

                         gAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGAt                                                               >  1:2712136/1‑93 (MQ=255)

                         gAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGAt                                                               >  2:1408494/1‑93 (MQ=255)

                         gAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGAt                                                               >  1:7436858/1‑93 (MQ=255)

                         gAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGAt                                                               <  2:4805136/93‑1 (MQ=255)

                          aaCTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATc                                                              <  1:4403811/93‑1 (MQ=255)

                          aaCTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTTCCGATc                                                              >  1:6877604/1‑93 (MQ=255)

                          aaCTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATc                                                              <  2:6754175/93‑1 (MQ=255)

                              gCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCAc                                                          >  2:2530728/1‑93 (MQ=255)

                                cGTGTAATGCGGATTTCTGGCTCACAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCg                                                        >  2:4612001/1‑93 (MQ=255)

                                 gtgtAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGt                                                       <  1:1182866/93‑1 (MQ=255)

                                  tgtAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTg                                                      >  2:7024592/1‑93 (MQ=255)

                                  tgtAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTg                                                      >  1:2437735/1‑93 (MQ=255)

                                  tgtAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTg                                                      <  2:402591/93‑1 (MQ=255)

                                  tgtAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTg                                                      >  1:2758592/1‑93 (MQ=255)

                                   gtAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTgg                                                     >  1:785355/1‑93 (MQ=255)

                                    tAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGt                                                    <  1:4173348/93‑1 (MQ=255)

                                    tAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGt                                                    <  1:851553/93‑1 (MQ=255)

                                       tGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGc                                                 <  1:1924494/93‑1 (MQ=255)

                                       tGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGc                                                 <  1:3710704/93‑1 (MQ=255)

                                        gCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGcc                                                <  2:5717292/93‑1 (MQ=255)

                                        gCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGcc                                                <  2:2412111/93‑1 (MQ=255)

                                        gCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGcc                                                <  1:2082742/93‑1 (MQ=255)

                                          ggATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCgg                                              <  1:239335/93‑1 (MQ=255)

                                          ggATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCgg                                              >  2:3602869/1‑93 (MQ=255)

                                             tttCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCa                                           >  1:1297206/1‑93 (MQ=255)

                                             tttCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCa                                           >  1:3206427/1‑93 (MQ=255)

                                             tttCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCa                                           <  2:7457683/93‑1 (MQ=255)

                                             tttCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCa                                           <  2:5865980/93‑1 (MQ=255)

                                             nntCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCa                                           <  2:1299217/91‑1 (MQ=255)

                                                 tGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGt                                       <  1:5418434/93‑1 (MQ=255)

                                                   gCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGttt                                     <  1:571032/93‑1 (MQ=255)

                                                    cTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGtttt                                    <  1:400261/93‑1 (MQ=255)

                                                    cTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGtttt                                    <  2:6592641/93‑1 (MQ=255)

                                                     tCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtg                                   >  2:3498724/1‑93 (MQ=255)

                                                      cGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtgg                                  <  1:941265/93‑1 (MQ=255)

                                                      cGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtgg                                  <  1:1759926/93‑1 (MQ=255)

                                                       gCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggt                                 >  1:6457026/1‑93 (MQ=255)

                                                       gCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggt                                 <  1:3773678/93‑1 (MQ=255)

                                                       gCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggt                                 >  1:3491400/1‑93 (MQ=255)

                                                       gCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggt                                 <  1:5857498/93‑1 (MQ=255)

                                                        cAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggtg                                >  2:5400365/1‑93 (MQ=255)

                                                         aGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggtgg                               >  2:6976293/1‑93 (MQ=255)

                                                         aGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggtgg                               <  2:4973965/93‑1 (MQ=255)

                                                         aGTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTtggtgg                               <  1:6075279/93‑1 (MQ=255)

                                                          gTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTggg                              >  2:5592021/1‑93 (MQ=255)

                                                          gTACAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTggg                              >  2:5332380/1‑93 (MQ=255)

                                                            aCAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCt                            >  1:5595034/1‑93 (MQ=255)

                                                            aCAGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCt                            >  1:7049775/1‑93 (MQ=255)

                                                              aGCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCttt                          >  2:2989244/1‑93 (MQ=255)

                                                               gCAAGAGTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTc                         <  1:6669543/93‑1 (MQ=255)

                                                                  agagTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCcgc                      >  2:1782849/1‑93 (MQ=255)

                                                                    agTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCgcgc                    <  2:184248/93‑1 (MQ=255)

                                                                    agTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCgcgc                    <  1:2601894/93‑1 (MQ=255)

                                                                     gTTATTACCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCt                   >  2:3518087/1‑93 (MQ=255)

                                                                      ttattaCCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTg                  >  2:1623036/1‑93 (MQ=255)

                                                                      ttattaCCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTg                  <  2:6340788/93‑1 (MQ=255)

                                                                      ttattaCCCCTGGGGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTg                  >  1:990309/1‑93 (MQ=255)

                                                                        attaCCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTa                <  2:5937002/93‑1 (MQ=255)

                                                                        attaCCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTa                <  2:6212150/93‑1 (MQ=255)

                                                                           aCCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGc             <  1:2373626/93‑1 (MQ=255)

                                                                           aCCCCTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGc             <  1:4668049/93‑1 (MQ=255)

                                                                             cccTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCg           <  2:1031657/93‑1 (MQ=255)

                                                                             cccTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCg           <  1:3123668/93‑1 (MQ=255)

                                                                               cTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGa         >  1:6208012/1‑93 (MQ=255)

                                                                               cTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGa         >  2:6156884/1‑93 (MQ=255)

                                                                               cTGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGa         <  1:2019718/93‑1 (MQ=255)

                                                                                tGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGAc        <  2:4564269/93‑1 (MQ=255)

                                                                                tGGTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGAc        >  2:2637871/1‑93 (MQ=255)

                                                                                 ggTGAAAGCCATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACt       >  1:1614333/1‑93 (MQ=255)

                                                                                  gTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTg      >  2:1341695/1‑93 (MQ=255)

                                                                                      aaGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACt  >  2:3847503/1‑93 (MQ=255)

TATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGCACCGTGGTTGCCGGGCAGAGTTTTGGTGGGCTTTCCGCGCTGTATGCCGGACTGCACT  >  NC_000913/612661‑612839

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