BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,187,979	C→T	A341T (GCC→ACC)	phoQ ←	PhoQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,187,979	0	C	T	93.5%	369.6 / 6.2	125	A341T (GCC→ACC)	phoQ	PhoQ
Reads supporting (aligned to +/- strand): ref base C (4/4); new base T (54/62); total (59/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

Only 100 of 132 total aligned reads displayed.

TGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGA  >  NC_000913/1187889‑1188068

tGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGa                                                                                         >  1:2663210/1‑93 (MQ=255)

 gATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGAc                                                                                        <  2:1654930/93‑1 (MQ=255)

   tctcTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCg                                                                                      >  1:671323/1‑93 (MQ=255)

     tctGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGa                                                                                    >  2:1801255/1‑93 (MQ=255)

       tGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATg                                                                                  >  2:646176/1‑93 (MQ=255)

         gCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCa                                                                                <  2:1692065/93‑1 (MQ=255)

           gAAATATCGAGAGCGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGc                                                                              >  2:3638374/1‑93 (MQ=255)

           gAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGc                                                                              <  2:1283140/93‑1 (MQ=255)

            aaaTATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCt                                                                             <  2:4344562/93‑1 (MQ=255)

             aaTATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTc                                                                            >  2:2695222/1‑93 (MQ=255)

              atatCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTcg                                                                           >  1:401520/1‑93 (MQ=255)

                 tCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCgg                                                                        >  1:954177/1‑93 (MQ=255)

                 tCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCgg                                                                        >  2:3219646/1‑93 (MQ=255)

                 tCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCgg                                                                        >  2:666241/1‑93 (MQ=255)

                  cGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGc                                                                       <  1:4050886/93‑1 (MQ=255)

                   gagagagaTATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCt                                                                      >  1:2432795/1‑93 (MQ=255)

                    agagagaTATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTg                                                                     >  1:1301471/1‑93 (MQ=255)

                    agagagaTATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTg                                                                     >  2:1159219/1‑93 (MQ=255)

                     gagagaTATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTgg                                                                    >  1:1545497/1‑92 (MQ=255)

                     gagagaTATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGa                                                                    >  1:2744343/1‑93 (MQ=255)

                     gagagaTATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGa                                                                    <  1:3962886/93‑1 (MQ=255)

                       gagaTATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGc                                                                  >  1:2542282/1‑93 (MQ=255)

                       gagaTATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGc                                                                  <  2:502410/93‑1 (MQ=255)

                        agaTATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCa                                                                 <  2:1947367/93‑1 (MQ=255)

                            atTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACAGGATGCAGCTCGCGGCTGAGCAATGt                                                             <  2:673756/93‑1 (MQ=255)

                               gACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTccc                                                          <  1:1976576/93‑1 (MQ=255)

                                aCCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGGAATGTCccg                                                         >  1:3956082/1‑93 (MQ=255)

                                aCCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGAGGCAGCTCGCGGCTGAGCAATGTCccg                                                         >  2:3795983/1‑93 (MQ=255)

                                 ccccTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCccgc                                                        >  2:1524234/1‑93 (MQ=255)

                                 ccccTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCccgc                                                        >  1:3728691/1‑93 (MQ=255)

                                 ccccTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCccgc                                                        >  1:3511660/1‑93 (MQ=255)

                                  cccTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCccgcc                                                       <  2:1077055/93‑1 (MQ=255)

                                  cccTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGGGACCGGATGCGGCTCGCGGCTGAGCAGTGGccccgc                                                       >  2:963892/1‑90 (MQ=255)

                                   ccTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCccgccg                                                      >  2:4660038/1‑93 (MQ=255)

                                    cTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCcgccgc                                                     >  1:1053654/1‑93 (MQ=255)

                                    cTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCcgccgc                                                     <  1:2841351/93‑1 (MQ=255)

                                     tttGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCcgcg                                                    >  2:2547982/1‑93 (MQ=255)

                                       tGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCa                                                  <  2:2852198/93‑1 (MQ=255)

                                          gTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATAc                                               >  1:2955255/1‑93 (MQ=255)

                                             gATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTgg                                            >  1:3088492/1‑93 (MQ=255)

                                             gATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTgg                                            <  1:676754/93‑1 (MQ=255)

                                             gATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCACGCCGCGGCTACTgg                                            >  2:2961966/1‑93 (MQ=255)

                                               tACACTTTGTTCAGCGCTGTGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCa                                          >  2:4562495/1‑93 (MQ=255)

                                               tACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCa                                          <  1:2054425/93‑1 (MQ=255)

                                                acacTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCAc                                         <  2:1295044/93‑1 (MQ=255)

                                                  acTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGTAATGTCCCGCCGCGCATACTGGCACgt                                       >  2:1436601/1‑92 (MQ=255)

                                                  acTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGa                                       >  1:3658079/1‑93 (MQ=255)

                                                  acTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGa                                       <  2:3325203/93‑1 (MQ=255)

                                                  acTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGa                                       <  2:3725426/93‑1 (MQ=255)

                                                   cTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGAt                                      <  2:2586805/93‑1 (MQ=255)

                                                     ttGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTTACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGc                                    <  2:3543816/93‑1 (MQ=255)

                                                           aGCGCTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGgtag                              <  2:3959385/93‑1 (MQ=255)

                                                            gcgcTGAGGTGAGAATGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGgtagt                             <  2:4342063/93‑1 (MQ=255)

                                                             cgcTGAGGTGAGATTGTCCAGCAGTGGTGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGgtagta                            <  1:3302946/93‑1 (MQ=255)

                                                             cgcTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGgtagta                            <  1:3545230/93‑1 (MQ=255)

                                                             cgcTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGgtagta                            <  1:2392190/93‑1 (MQ=255)

                                                             cgcTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGgtagta                            <  1:1328619/93‑1 (MQ=255)

                                                              gcTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCTCCGCGCATACTGGCACGATGCAGGtagtag                           <  2:425771/93‑1 (MQ=255)

                                                              gcTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGtagtag                           >  2:3698995/1‑93 (MQ=255)

                                                              gcTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGtagtag                           >  1:4204267/1‑93 (MQ=255)

                                                               cTGAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGc                          <  2:3937033/93‑1 (MQ=255)

                                                                 gAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCa                        <  1:1827871/93‑1 (MQ=255)

                                                                 gAGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCa                        >  1:3947193/1‑93 (MQ=255)

                                                                  aGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCaa                       <  1:3916352/93‑1 (MQ=255)

                                                                  aGGTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCaa                       >  1:309556/1‑93 (MQ=255)

                                                                   ggTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAAt                      <  2:2406605/93‑1 (MQ=255)

                                                                   ggTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAAt                      <  2:2463645/93‑1 (MQ=255)

                                                                    gttggaTTGTCCAGCAGTTGGGTGACCGGTTGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAAtt                     <  1:460172/89‑1 (MQ=255)

                                                                    gTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAAtt                     >  2:1261808/1‑93 (MQ=255)

                                                                    gTGAGATTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAAtt                     <  2:1936901/93‑1 (MQ=255)

                                                                    gTGAGAGTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAAtt                     <  2:2165034/93‑1 (MQ=255)

                                                                       agaTTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGc                  <  1:2558171/93‑1 (MQ=255)

                                                                       agaTTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGc                  >  1:1961466/1‑93 (MQ=255)

                                                                       agaTTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGc                  <  2:3584106/93‑1 (MQ=255)

                                                                       agaTTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGGAGTAGCCAAATTGc                  >  1:4241518/1‑93 (MQ=255)

                                                                       agaTTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGc                  <  2:3706890/93‑1 (MQ=255)

                                                                        gaTTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCt                 >  2:2315419/1‑93 (MQ=255)

                                                                        gaTTGTCCAGCAGTGCGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCt                 <  2:3474117/93‑1 (MQ=255)

                                                                         aTTGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCtg                <  2:3421026/93‑1 (MQ=255)

                                                                          ttGTCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCtgt               <  2:1786705/93‑1 (MQ=255)

                                                                          ttGGCCAGCAGGCGGGTGACAGGATGCAGCTCGCGGCTGAGCAAGGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCtgt               <  2:1716020/93‑1 (MQ=255)

                                                                             tCCAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGaa            <  1:672076/93‑1 (MQ=255)

                                                                              ccAGCAGTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGaaa           <  1:513987/93‑1 (MQ=255)

                                                                               cagcagTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAAt          <  1:1933241/93‑1 (MQ=255)

                                                                               cagcagTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAAt          <  2:1913001/93‑1 (MQ=255)

                                                                               cagcagTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAAt          <  2:2051768/93‑1 (MQ=255)

                                                                                 gcagTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATgg        >  2:3460825/1‑92 (MQ=255)

                                                                                 gcagTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGc        <  1:2360162/93‑1 (MQ=255)

                                                                                 gcagTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGc        <  1:3314264/93‑1 (MQ=255)

                                                                                  cagTGGGGTGGCCGGATGCACCTCGCGGCTGAGCAATGTCCCGCCGCGCATGCTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCg       <  1:4228009/93‑1 (MQ=255)

                                                                                  cagTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCg       <  1:3360116/93‑1 (MQ=255)

                                                                                  cagTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCg       <  2:2160425/93‑1 (MQ=255)

                                                                                   agTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGGATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCgg      >  2:1509758/1‑93 (MQ=255)

                                                                                   agTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCgg      <  1:2568016/93‑1 (MQ=255)

                                                                                    gTGGGGTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGc     <  2:753826/93‑1 (MQ=255)

                                                                                    gTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGc     >  1:2268916/1‑93 (MQ=255)

                                                                                    gTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGGAGTAGCCAATTTGCTGTGAAATGCGGc     >  2:1462938/1‑93 (MQ=255)

                                                                                     gggggTGACCGGATGCACCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCt    <  2:775687/92‑1 (MQ=255)

                                                                                       gggTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGa  <  1:4384068/93‑1 (MQ=255)

                                                                                       gggTGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGa  <  2:180342/93‑1 (MQ=255)

TGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGA  >  NC_000913/1187889‑1188068

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