BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,999,772	T→C	H4R (CAC→CGC)	b1922 ←	sigma factor 28

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,999,772	0	T	C	99.2%	391.9 / ‑5.5	127	H4R (CAC→CGC)	b1922	sigma factor 28
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (40/86); total (40/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 136 total aligned reads displayed.

AGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCT  >  NC_000913/1999682‑1999861

aGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAgcgt                                                                                         <  2:2630153/93‑1 (MQ=255)

 gATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGtt                                                                                        >  2:1724972/1‑93 (MQ=255)

  aTCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGttt                                                                                       <  1:3976355/93‑1 (MQ=255)

    cGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTAt                                                                                     >  1:1543503/1‑93 (MQ=255)

     gTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATc                                                                                    <  1:3131407/93‑1 (MQ=255)

     gTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGATTATc                                                                                    >  2:3389079/1‑93 (MQ=255)

        aaGTTCCATGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCCGGGCATAACGCTGCCACAGCGAGCGTTTATTCAt                                                                                 >  2:4374128/1‑93 (MQ=255)

        aaGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCCCCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCAt                                                                                 <  2:934171/93‑1 (MQ=255)

         aGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCAtt                                                                                <  2:2791033/93‑1 (MQ=255)

          gTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTa                                                                               >  2:3536032/1‑93 (MQ=255)

           ttCCACGCTCTCGGGCTGTCGAACCTGCAGGCGCTATGCTTCGTGACGCACCAGGGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTac                                                                              <  2:3289395/93‑1 (MQ=255)

              cACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACAcc                                                                           <  1:67715/93‑1 (MQ=255)

                cGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGTCATAACGCTGCCACAGCGAGCGTTTATCCATTACACCtt                                                                         <  2:493824/93‑1 (MQ=255)

                cGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCtt                                                                         <  1:1681493/93‑1 (MQ=255)

                cGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCtt                                                                         <  1:2046877/93‑1 (MQ=255)

                 gCTCGCGGGCAGTCGAACCTGCAGGCACAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTc                                                                        <  2:4142500/93‑1 (MQ=255)

                  cTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCa                                                                       <  1:833108/93‑1 (MQ=255)

                    cgcgGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGc                                                                     <  1:4247389/93‑1 (MQ=255)

                      cgGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCgg                                                                   >  1:2407553/1‑93 (MQ=255)

                       gggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCgga                                                                  >  1:353432/1‑92 (MQ=255)

                       gggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGt                                                                  <  1:3414957/93‑1 (MQ=255)

                       gggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGt                                                                  >  1:995785/1‑93 (MQ=255)

                       gggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGt                                                                  <  2:3264186/93‑1 (MQ=255)

                       gggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGt                                                                  >  1:4503156/1‑93 (MQ=255)

                        ggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGta                                                                 <  2:2296801/93‑1 (MQ=255)

                        ggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGta                                                                 <  1:2145145/93‑1 (MQ=255)

                        ggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGta                                                                 <  2:3378865/93‑1 (MQ=255)

                         gCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGtat                                                                >  2:2461244/1‑93 (MQ=255)

                          cAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGtata                                                               <  1:402898/93‑1 (MQ=255)

                           aGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATag                                                              >  1:1913426/1‑93 (MQ=255)

                            gTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATaga                                                             <  2:2753812/93‑1 (MQ=255)

                            gTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATaga                                                             <  1:2770865/93‑1 (MQ=255)

                            gTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATaga                                                             <  2:2680261/93‑1 (MQ=255)

                              cGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGt                                                           >  1:973298/1‑93 (MQ=255)

                              cGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGt                                                           <  2:1668284/93‑1 (MQ=255)

                                 aCCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATCACGCTGCCACAGCGAGCGTTTATACATTACACCTTCAGCGGTATAGAGTGaa                                                        >  1:666782/1‑93 (MQ=255)

                                 aCCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGaa                                                        >  1:3603444/1‑93 (MQ=255)

                                 aCCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGaa                                                        <  2:2547868/93‑1 (MQ=255)

                                   cTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAAtt                                                      >  2:915876/1‑93 (MQ=255)

                                    tGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTc                                                     >  1:347867/1‑93 (MQ=255)

                                    tGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTc                                                     >  2:1009714/1‑93 (MQ=255)

                                     gCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCa                                                    <  1:2632528/93‑1 (MQ=255)

                                     gCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCa                                                    <  1:2609669/93‑1 (MQ=255)

                                      gaGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCAc                                                   <  2:1167187/92‑1 (MQ=255)

                                      cAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCAc                                                   <  2:1881100/93‑1 (MQ=255)

                                      cAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCAc                                                   <  1:570666/93‑1 (MQ=255)

                                      cAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCAc                                                   <  2:3576128/93‑1 (MQ=255)

                                        ggCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACgt                                                 >  2:2082185/1‑92 (MQ=255)

                                        ggCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGa                                                 >  1:3451174/1‑93 (MQ=255)

                                        ggCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGa                                                 <  2:336379/93‑1 (MQ=255)

                                          cgcAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATa                                               <  1:2554281/93‑1 (MQ=255)

                                          cgcAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATa                                               <  2:4653906/93‑1 (MQ=255)

                                           gcAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATaa                                              <  2:3766298/93‑1 (MQ=255)

                                             aaTGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAAc                                            >  1:363425/1‑93 (MQ=255)

                                              aTGCTTGGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACa                                           <  2:3501433/93‑1 (MQ=255)

                                              aTGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACa                                           <  1:1093040/93‑1 (MQ=255)

                                                  ttCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGccc                                       <  2:924807/93‑1 (MQ=255)

                                                  ttCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGccc                                       >  2:1616752/1‑93 (MQ=255)

                                                  ttCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGccc                                       >  2:4272484/1‑93 (MQ=255)

                                                   tCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCt                                      <  2:4675640/93‑1 (MQ=255)

                                                      tGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCg                                   >  2:1541824/1‑93 (MQ=255)

                                                       gACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGt                                  <  2:1123545/93‑1 (MQ=255)

                                                       gACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGt                                  <  2:1414623/93‑1 (MQ=255)

                                                       gACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGt                                  >  1:4303372/1‑93 (MQ=255)

                                                          gCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTtat                               >  2:853683/1‑93 (MQ=255)

                                                          ccACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTtat                               <  1:1432310/92‑1 (MQ=255)

                                                           cACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTtata                              <  2:1973609/93‑1 (MQ=255)

                                                           cACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTtata                              <  1:912135/93‑1 (MQ=255)

                                                           cACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTtata                              <  1:4241193/93‑1 (MQ=255)

                                                            aCCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTatat                             <  1:4548037/93‑1 (MQ=255)

                                                             ccAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATg                            <  1:926725/93‑1 (MQ=255)

                                                              cAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGa                           <  2:397175/93‑1 (MQ=255)

                                                                gCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGt                         <  2:667947/93‑1 (MQ=255)

                                                                 cGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGtt                        >  1:300450/1‑93 (MQ=255)

                                                                  gggACATAACGCTGCCACAGCGAGCGTTTATCCATTACCCCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTa                       <  1:4394493/93‑1 (MQ=255)

                                                                  gggACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTa                       >  1:704881/1‑93 (MQ=255)

                                                                   ggACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAt                      >  2:1103496/1‑93 (MQ=255)

                                                                   ggACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAt                      >  1:1286226/1‑93 (MQ=255)

                                                                   ggACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAt                      >  1:139136/1‑93 (MQ=255)

                                                                    gACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATc                     <  2:4254964/93‑1 (MQ=255)

                                                                    gACAAACCGCGGCCCGAGCGAGCGTTTATCCATTACACCTGCAGCGGGATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATc                     <  1:1238084/93‑1 (MQ=255)

                                                                     aCATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCg                    <  2:3086361/93‑1 (MQ=255)

                                                                     aCATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCg                    <  1:3838985/93‑1 (MQ=255)

                                                                        tAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCa                 <  1:337392/93‑1 (MQ=255)

                                                                        tAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCa                 <  1:1012174/93‑1 (MQ=255)

                                                                          aCGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATg               <  1:1361663/93‑1 (MQ=255)

                                                                           cGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGa              >  2:2017223/1‑93 (MQ=255)

                                                                            gCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGAt             <  2:2186400/93‑1 (MQ=255)

                                                                            gCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGAt             <  1:771474/93‑1 (MQ=255)

                                                                               gCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTAt          <  2:3887231/93‑1 (MQ=255)

                                                                                ccACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATc         <  1:1411657/93‑1 (MQ=255)

                                                                                ccACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATc         <  1:4123393/93‑1 (MQ=255)

                                                                                  acaGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATTCg       >  1:1039636/1‑93 (MQ=255)

                                                                                  acaGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCg       <  2:4120368/93‑1 (MQ=255)

                                                                                   caGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGt      <  1:2464938/93‑1 (MQ=255)

                                                                                    agcgagcgTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGtt     <  1:2016715/93‑1 (MQ=255)

                                                                                     gCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGttt    <  2:2675104/93‑1 (MQ=255)

                                                                                     gCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGttt    >  1:875925/1‑93 (MQ=255)

                                                                                       gAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCt  <  1:3291224/93‑1 (MQ=255)

                                                                                       gAGCGTTTATCCAGTACTCCTTCAGCGGTATAGATTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCTTTTCt  <  2:1403510/93‑1 (MQ=255)

AGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCT  >  NC_000913/1999682‑1999861

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