BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,977,219	T→C	H59H (CAT→CAC)	lysR →	LysR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,977,219	0	T	C	99.5%	1277.2 / ‑3.4	394	H59H (CAT→CAC)	lysR	LysR
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (157/235); total (157/237)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.14e-01

Only 100 of 409 total aligned reads displayed.

CACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTCTGC  >  NC_000913/2977129‑2977307

cacCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACAccc                                                                                        <  2:10037386/93‑1 (MQ=255)

cacCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACAccc                                                                                        <  2:7147382/93‑1 (MQ=255)

   cTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTAc                                                                                     <  1:9443517/93‑1 (MQ=255)

   cTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTAc                                                                                     <  1:428096/93‑1 (MQ=255)

   cTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTAc                                                                                     <  2:2840857/93‑1 (MQ=255)

    tCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTAcc                                                                                    <  2:9371375/93‑1 (MQ=255)

    tCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTAcc                                                                                    <  2:2075920/93‑1 (MQ=255)

    tCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTAcc                                                                                    <  2:2603074/93‑1 (MQ=255)

     cacaGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCg                                                                                   <  2:8684728/93‑1 (MQ=255)

        aGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGc                                                                                >  2:1909057/1‑93 (MQ=255)

         gCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCa                                                                               <  1:6816702/93‑1 (MQ=255)

           cAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTTCCGTGCAAg                                                                             <  2:8557045/93‑1 (MQ=255)

            aaCCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAgg                                                                            <  1:6507915/93‑1 (MQ=255)

            aaCCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAgg                                                                            <  1:876949/93‑1 (MQ=255)

            aaCCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAgg                                                                            >  1:5039005/1‑93 (MQ=255)

            aaCCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAgg                                                                            <  1:2022529/93‑1 (MQ=255)

             aCCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGa                                                                           >  2:3482154/1‑93 (MQ=255)

             aCCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGa                                                                           <  2:8299785/93‑1 (MQ=255)

             aCCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGa                                                                           >  2:1280844/1‑93 (MQ=255)

             aCCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGa                                                                           <  2:9519450/93‑1 (MQ=255)

              ccGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGTCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGAc                                                                          <  1:2900402/93‑1 (MQ=255)

               cGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACt                                                                         <  1:778323/93‑1 (MQ=255)

                gTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTg                                                                        >  2:3149460/1‑93 (MQ=255)

                gTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTg                                                                        <  2:1576417/93‑1 (MQ=255)

                gTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTg                                                                        >  2:987160/1‑93 (MQ=255)

                 tCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGc                                                                       <  1:9194772/93‑1 (MQ=255)

                  cAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCg                                                                      >  1:2108143/1‑93 (MQ=255)

                   aGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGt                                                                     <  2:3229233/93‑1 (MQ=255)

                   aGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGt                                                                     >  1:1828508/1‑93 (MQ=255)

                    gCCGCGATCGTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTc                                                                    <  2:1104208/93‑1 (MQ=255)

                    gCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTc                                                                    <  1:7292963/93‑1 (MQ=255)

                    gCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTc                                                                    <  1:2639785/93‑1 (MQ=255)

                    gCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTc                                                                    >  2:1690392/1‑93 (MQ=255)

                       gcgAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGt                                                                 >  2:7848930/1‑93 (MQ=255)

                        cgAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGtt                                                                <  1:4494273/93‑1 (MQ=255)

                         gAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACGCCCTACCGTGCAAGGACTGCGTCTGttt                                                               <  2:8805367/93‑1 (MQ=255)

                         gAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGttt                                                               <  2:5735823/93‑1 (MQ=255)

                            cTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTgaa                                                            >  2:137033/1‑93 (MQ=255)

                             ttGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTgaag                                                           >  2:6103565/1‑93 (MQ=255)

                             ttGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTgaag                                                           >  1:368078/1‑93 (MQ=255)

                              tGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTgaaga                                                          <  2:7356309/93‑1 (MQ=255)

                               gcgcgcTTTGAGAAGGTTATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTgaagaa                                                         <  2:4365966/93‑1 (MQ=255)

                                 gcgcTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGt                                                       <  2:3380197/93‑1 (MQ=255)

                                    cTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCa                                                    >  2:639877/1‑93 (MQ=255)

                                     tttGAGAAGGTGATCGGGCTGAAATTGTTTGATCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCaa                                                   >  1:9392409/1‑93 (MQ=255)

                                     tttGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCaa                                                   <  1:1408041/93‑1 (MQ=255)

                                      ttGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAAc                                                  <  2:5584819/93‑1 (MQ=255)

                                      ttGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAAc                                                  <  1:6271736/93‑1 (MQ=255)

                                      ttGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAAc                                                  <  1:6976744/93‑1 (MQ=255)

                                        gagaAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGa                                                <  2:5452014/93‑1 (MQ=255)

                                        gagaAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGa                                                <  2:6017837/93‑1 (MQ=255)

                                         agaAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCACCGAt                                               >  1:8459713/1‑93 (MQ=255)

                                         agaAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGAt                                               >  1:8535905/1‑93 (MQ=255)

                                         agaAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGAt                                               >  2:309318/1‑93 (MQ=255)

                                          gagggTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGAGCAAGGACTGCGGCTGTTTGAAGAAGTGCAACGATc                                              >  2:4542295/4‑93 (MQ=255)

                                          gaAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATc                                              <  2:6288918/93‑1 (MQ=255)

                                          gaAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATc                                              <  1:4326335/93‑1 (MQ=255)

                                          gaAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATc                                              <  1:6655959/93‑1 (MQ=255)

                                          gaAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATc                                              >  2:9136016/1‑93 (MQ=255)

                                           aaGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATcc                                             <  2:1471709/93‑1 (MQ=255)

                                           aaGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATcc                                             <  1:4441476/93‑1 (MQ=255)

                                           aaGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATcc                                             <  1:3430994/93‑1 (MQ=255)

                                            aGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCt                                            >  1:9767960/1‑93 (MQ=255)

                                             ggTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTg                                           <  1:8912265/93‑1 (MQ=255)

                                               tGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGCCTGACTCTGTTTGAAGAAGTGCAACGATCCTgga                                         >  1:769598/1‑92 (MQ=255)

                                                    gggCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGa                                    >  2:8187282/1‑93 (MQ=255)

                                                     ggCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGAc                                   <  1:4733490/93‑1 (MQ=255)

                                                        tGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTgg                                <  1:2336146/93‑1 (MQ=255)

                                                        tGAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTgg                                >  1:10043851/1‑93 (MQ=255)

                                                         gAAATTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGa                               >  2:5150223/1‑93 (MQ=255)

                                                          aaaTTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGAt                              >  2:4654612/1‑93 (MQ=255)

                                                           aaTTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATc                             <  2:3720415/93‑1 (MQ=255)

                                                           aaTTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATc                             <  1:1747469/93‑1 (MQ=255)

                                                           aaTTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATc                             <  2:5264605/93‑1 (MQ=255)

                                                           aaNTGTTTGNGCGCGTACGTGNNCGNTTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATc                             <  1:605527/93‑1 (MQ=255)

                                                            aTTGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCg                            <  2:5001018/93‑1 (MQ=255)

                                                              tGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCa                          >  2:4210796/1‑93 (MQ=255)

                                                              tGTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCa                          >  1:1323042/1‑93 (MQ=255)

                                                               gTTTGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCAt                         <  1:8159364/93‑1 (MQ=255)

                                                                 ttGAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTg                       <  2:2934542/93‑1 (MQ=255)

                                                                   gAGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTc                     <  2:6716174/93‑1 (MQ=255)

                                                                    aGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCa                    <  2:8427677/93‑1 (MQ=255)

                                                                    aGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCa                    >  1:5453614/1‑93 (MQ=255)

                                                                    aGCGCGTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCa                    >  2:7698305/1‑93 (MQ=255)

                                                                       gcgTACGTGTGCGATTACACCCTGCCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGTACTGGATCGCATTGTCAgcg                 <  2:6837273/93‑1 (MQ=255)

                                                                       gcgTACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAgcg                 >  2:2322712/1‑93 (MQ=255)

                                                                          tACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGcgccg              >  2:4115920/1‑93 (MQ=255)

                                                                          tACGTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGcgccg              >  1:5208640/1‑93 (MQ=255)

                                                                             gTGGGCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAg           <  1:1137345/93‑1 (MQ=255)

                                                                               gggCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGaa         >  1:302912/1‑93 (MQ=255)

                                                                                ggCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGaaa        >  1:95611/1‑93 (MQ=255)

                                                                                ggCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGaaa        <  1:4823685/93‑1 (MQ=255)

                                                                                 gCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAg       <  1:3335446/93‑1 (MQ=255)

                                                                                 gCGATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAg       >  2:8898540/1‑93 (MQ=255)

                                                                                   gATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTc     <  1:5885215/93‑1 (MQ=255)

                                                                                   gATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTc     <  1:3874074/93‑1 (MQ=255)

                                                                                   gATTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTc     <  1:3984533/93‑1 (MQ=255)

                                                                                    aTTACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTCt    >  1:8229326/1‑93 (MQ=255)

                                                                                     ttACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTCTg   >  1:7872019/1‑93 (MQ=255)

                                                                                     ttACACCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTCTg   <  1:8352801/93‑1 (MQ=255)

CACCTCACAGCCAACCGTCAGCCGCGAACTTGCGCGCTTTGAGAAGGTGATCGGGCTGAAATTGTTTGAGCGCGTACGTGGGCGATTACATCCTACCGTGCAAGGACTGCGTCTGTTTGAAGAAGTGCAACGATCCTGGTACGGACTGGATCGCATTGTCAGCGCCGCAGAAAGTCTGC  >  NC_000913/2977129‑2977307

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