BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,654,241	A→G	R46R (CGT→CGC)	b2523 ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,654,241	0	A	G	85.7%	634.0 / 85.2	266	R46R (CGT→CGC)	b2523	putative peptidase
Reads supporting (aligned to +/- strand): ref base A (19/19); new base G (118/110); total (137/129)
Fisher's exact test for biased strand distribution p-value = 8.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 273 total aligned reads displayed.

TTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCG  >  NC_000913/2654151‑2654330

ttGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCcgcgcg                                                                                         >  2:5891219/1‑93 (MQ=255)

   ccAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTg                                                                                      <  2:4951829/93‑1 (MQ=255)

    cAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTgg                                                                                     >  1:5050612/1‑93 (MQ=255)

      gAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGAt                                                                                   <  1:5969138/93‑1 (MQ=255)

       aaTGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATc                                                                                  >  2:3719259/1‑93 (MQ=255)

       aaTGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATc                                                                                  >  2:6422270/1‑93 (MQ=255)

        aTGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCa                                                                                 <  2:6022589/93‑1 (MQ=255)

        aTGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCa                                                                                 <  2:6533704/93‑1 (MQ=255)

           cccAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGcc                                                                              >  1:3630983/1‑93 (MQ=255)

              agcagcGATCCGCGTCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCa                                                                           <  2:1861734/93‑1 (MQ=255)

               gcagcGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAg                                                                          >  2:386392/1‑93 (MQ=255)

                cagcGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGa                                                                         >  1:235226/1‑93 (MQ=255)

                 agcGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGAt                                                                        >  2:6369170/1‑93 (MQ=255)

                 agagATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGAt                                                                        >  1:836401/4‑93 (MQ=255)

                  gcGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATc                                                                       >  2:3479687/1‑93 (MQ=255)

                  gcGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATc                                                                       >  1:655170/1‑93 (MQ=255)

                  gcGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATc                                                                       >  2:1756388/1‑93 (MQ=255)

                   cGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCg                                                                      <  2:291600/93‑1 (MQ=255)

                     aTCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTc                                                                    <  1:3812532/93‑1 (MQ=255)

                     aTCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTc                                                                    >  2:6810353/1‑93 (MQ=255)

                     aTCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTc                                                                    <  1:4594660/93‑1 (MQ=255)

                      tCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCt                                                                   >  1:5125271/1‑93 (MQ=255)

                      tCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCt                                                                   >  2:4843538/1‑93 (MQ=255)

                      tCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCt                                                                   >  2:900512/1‑93 (MQ=255)

                       ccGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTg                                                                  >  2:5645393/1‑93 (MQ=255)

                        cGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGc                                                                 >  2:2746124/1‑93 (MQ=255)

                        cGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGc                                                                 >  2:3159571/1‑93 (MQ=255)

                         gCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGcc                                                                <  1:3345565/93‑1 (MQ=255)

                          cATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGccc                                                               <  2:7337147/93‑1 (MQ=255)

                            tCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCg                                                             >  2:6104087/1‑93 (MQ=255)

                               cAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTc                                                          >  1:6687883/1‑93 (MQ=255)

                               cAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTc                                                          <  1:1580134/93‑1 (MQ=255)

                                 gCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCaa                                                        >  1:4769101/1‑93 (MQ=255)

                                  ccTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCaaa                                                       <  2:6696966/93‑1 (MQ=255)

                                  ccTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCaaa                                                       <  2:2235204/93‑1 (MQ=255)

                                   cTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAAt                                                      <  1:6950505/93‑1 (MQ=255)

                                     tCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGc                                                    >  1:181130/1‑93 (MQ=255)

                                      cACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCa                                                   <  1:5788237/93‑1 (MQ=255)

                                      cACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCa                                                   <  2:5233378/93‑1 (MQ=255)

                                       aCCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAg                                                  <  2:4511210/93‑1 (MQ=255)

                                       aCCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAg                                                  >  2:2593624/1‑93 (MQ=255)

                                       aCCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAg                                                  <  1:4863209/93‑1 (MQ=255)

                                       aCCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAg                                                  >  2:4863430/1‑93 (MQ=255)

                                        ccGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAgg                                                 >  1:5740125/1‑93 (MQ=255)

                                         cGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAggg                                                <  1:4263380/93‑1 (MQ=255)

                                         cGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAggg                                                <  1:437929/93‑1 (MQ=255)

                                          gCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGt                                               >  1:893810/1‑93 (MQ=255)

                                          gCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGt                                               <  1:1105200/93‑1 (MQ=255)

                                          gCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGt                                               <  1:1381879/93‑1 (MQ=255)

                                           cTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTa                                              <  2:2150563/93‑1 (MQ=255)

                                            ttAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTaa                                             >  2:7447962/1‑93 (MQ=255)

                                            ttAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTaa                                             >  2:33589/1‑93 (MQ=255)

                                             tAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAAt                                            <  1:718760/93‑1 (MQ=255)

                                                cTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGcc                                         >  1:394825/1‑93 (MQ=255)

                                                 tGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCg                                        >  1:6640984/1‑93 (MQ=255)

                                                   aaCATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTc                                      <  2:4642121/93‑1 (MQ=255)

                                                    aCATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCa                                     >  1:1793735/1‑93 (MQ=255)

                                                     cATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCat                                    >  2:4352728/1‑93 (MQ=255)

                                                     cATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCat                                    >  1:6042470/1‑93 (MQ=255)

                                                      aTGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCatt                                   >  1:547061/1‑93 (MQ=255)

                                                       tGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCatta                                  <  2:249937/93‑1 (MQ=255)

                                                         cTTGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCattatt                                <  1:3472872/93‑1 (MQ=255)

                                                         cTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCattatt                                <  1:2612254/93‑1 (MQ=255)

                                                          ttGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCAttatta                               <  2:1676832/93‑1 (MQ=255)

                                                           tGATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTaa                              >  1:2764434/1‑93 (MQ=255)

                                                            gATGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAAt                             >  2:3580577/1‑93 (MQ=255)

                                                             aTGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATg                            >  1:6489993/1‑93 (MQ=255)

                                                              tGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGc                           >  1:6385152/1‑93 (MQ=255)

                                                              tGCCCAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGc                           <  2:103709/93‑1 (MQ=255)

                                                              tGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGc                           <  1:3120077/93‑1 (MQ=255)

                                                                cccAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTg                         <  2:1393134/93‑1 (MQ=255)

                                                                cccAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTg                         <  1:2293/93‑1 (MQ=255)

                                                                 ccAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGt                        <  1:5956882/93‑1 (MQ=255)

                                                                 ccAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGt                        >  1:7332272/1‑93 (MQ=255)

                                                                  cAGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTa                       <  2:2007371/93‑1 (MQ=255)

                                                                   aGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTaa                      <  2:7307704/93‑1 (MQ=255)

                                                                   aGACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTaa                      >  2:552217/1‑93 (MQ=255)

                                                                    gACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAg                     >  2:3918432/1‑93 (MQ=255)

                                                                    gACCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAg                     >  1:6788327/1‑93 (MQ=255)

                                                                     aCCGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGt                    >  1:4621489/1‑93 (MQ=255)

                                                                       cGTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTg                  >  2:1325842/1‑93 (MQ=255)

                                                                        gTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGc                 <  1:3910357/93‑1 (MQ=255)

                                                                        gTCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGc                 >  1:2110741/1‑93 (MQ=255)

                                                                         tCAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCt                >  2:1490309/1‑93 (MQ=255)

                                                                          cAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCtt               <  1:32953/93‑1 (MQ=255)

                                                                          cAATCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCtt               <  1:761492/93‑1 (MQ=255)

                                                                           aaTCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCttt              <  2:4251342/93‑1 (MQ=255)

                                                                            aTCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCTCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCtttt             >  1:1030975/1‑93 (MQ=255)

                                                                            aTCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCtttt             >  2:5909485/1‑93 (MQ=255)

                                                                            aTCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCtttt             <  2:5371052/93‑1 (MQ=255)

                                                                             tCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCttttt            >  2:7172625/1‑93 (MQ=255)

                                                                             tCTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCttttt            >  1:6202586/1‑93 (MQ=255)

                                                                              cTTGCGCGCCGCGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTtc           <  1:1493900/93‑1 (MQ=255)

                                                                              cTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTtc           <  2:6722540/93‑1 (MQ=255)

                                                                                    cgccgcGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCa     <  2:6856571/93‑1 (MQ=255)

                                                                                    cgccgcGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCa     <  1:5495150/93‑1 (MQ=255)

                                                                                     gccgcGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAg    >  2:6295762/1‑93 (MQ=255)

                                                                                     gccgcGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAg    <  2:5194219/93‑1 (MQ=255)

                                                                                      ccgcGCGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAgc   <  1:4877777/93‑1 (MQ=255)

                                                                                       cgcgcgCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAgcg  >  2:5715204/1‑93 (MQ=255)

TTGCCAGAATGCCCAGCAGCGATCCGCATCCCAGCCTTCACCGCTTAACTGAACATGCTTGATGCCCAGACCGTCAATCTTGCGCGCCGCACGCTGGATCAGCCCCAGATCGTCTGCCCCGTTCAAATGCAGGGTAATGCCGTCATTATTAATGCTGTAAGTTGCTTTTTCTCCCCAGCG  >  NC_000913/2654151‑2654330

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