BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,484,846	A→G	A1254A (GCA→GCG)	hrpA →	helicase, ATP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,484,846	0	A	G	99.2%	775.3 / ‑3.6	246	A1254A (GCA→GCG)	hrpA	helicase, ATP‑dependent
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (141/103); total (141/105)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 268 total aligned reads displayed.

GGTTGATCCACATCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCACGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCAACAACTTGGTACGCCTTA  >  NC_000913/1484756‑1484935

ggTTGATCCACATCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCcgcgcg                                                                                         <  1:3062412/93‑1 (MQ=255)

 gTTGATCCACATCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGcgt                                                                                        >  1:1650269/1‑93 (MQ=255)

 gTTGATCCACATCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGcgt                                                                                        >  1:3558700/1‑93 (MQ=255)

  ttGATCCACATCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGcgtc                                                                                       >  1:4577918/1‑93 (MQ=255)

   tGATCCACATCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGcgtcg                                                                                      <  2:5951307/93‑1 (MQ=255)

    ggtCCACATCGCGACCGGGCGCAGATGCTGAAAGTCGAAAACGTCCCGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGcgtcgt                                                                                     <  2:8288150/91‑1 (MQ=255)

    gATCCACATCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGcgtcgt                                                                                     <  2:491809/93‑1 (MQ=255)

    gATCCACATCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGcgtcgt                                                                                     >  2:3119602/1‑93 (MQ=255)

     aTCCACATCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTg                                                                                    >  1:3313630/1‑93 (MQ=255)

      tCCACATCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGa                                                                                   >  1:7028720/1‑93 (MQ=255)

       ccACATCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAg                                                                                  <  1:6201254/93‑1 (MQ=255)

       ccACATCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAg                                                                                  >  2:5900289/1‑93 (MQ=255)

       ccACATCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAg                                                                                  <  2:4551829/93‑1 (MQ=255)

       ccACATCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAg                                                                                  >  2:4362973/1‑93 (MQ=255)

         ccaTCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGa                                                                                <  1:954785/92‑1 (MQ=255)

          caTCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGAAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGAt                                                                               >  1:3424857/1‑93 (MQ=255)

          caTCACGACCGTACGCAGATGCTGAAAGTCGAAAACGTCCAGCGGGCGTGGCAGCAATGGATCAACAAACTGCCGCACGCGCGTCGTGAGGAt                                                                               <  2:8966305/93‑1 (MQ=255)

            tCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGa                                                                             >  2:4231370/1‑93 (MQ=255)

            tCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGa                                                                             <  1:4794474/93‑1 (MQ=255)

              gcgACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAg                                                                           >  2:6428659/1‑93 (MQ=255)

               cgACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGa                                                                          <  1:5647568/93‑1 (MQ=255)

                gACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCCATGGATCAACACACTGCCGCCCGCGCGCCGTGAGGATGAacac                                                                         >  2:6175960/1‑93 (MQ=255)

                gACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGAc                                                                         <  2:5292623/93‑1 (MQ=255)

                gACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGAGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGAc                                                                         >  1:3882737/1‑93 (MQ=255)

                 aCCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACg                                                                        >  1:9222543/1‑93 (MQ=255)

                    gTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGa                                                                     <  2:216792/93‑1 (MQ=255)

                    gTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGa                                                                     >  1:3667914/1‑93 (MQ=255)

                     tGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGaa                                                                    >  1:8631769/1‑93 (MQ=255)

                       cgcAGATGCTGAAAGTCGAAAACGTCCAGCAGGCTTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAag                                                                  <  1:9861053/93‑1 (MQ=255)

                       cgcAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAag                                                                  >  1:8074780/1‑93 (MQ=255)

                       cgcAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAag                                                                  <  1:8180547/93‑1 (MQ=255)

                        gcAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAaga                                                                 <  2:7972906/93‑1 (MQ=255)

                         cAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAagag                                                                >  2:811326/1‑93 (MQ=255)

                          aGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAgaga                                                               >  1:6120825/1‑93 (MQ=255)

                          aGATGCTCAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAgaga                                                               <  1:7836391/93‑1 (MQ=255)

                           gATGCTGAAAGTCGAAAACGTCCAGCAGGAGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGAt                                                              >  1:7372284/1‑93 (MQ=255)

                            aTGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATc                                                             <  1:2061240/93‑1 (MQ=255)

                             tGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATcc                                                            <  2:7367442/93‑1 (MQ=255)

                              gCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCg                                                           >  1:251262/1‑93 (MQ=255)

                               cTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGt                                                          <  1:9070710/93‑1 (MQ=255)

                                tGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGTGTCGTGAGGATGAAGACGTGAAAGAGATCCGtt                                                         >  1:9671858/1‑93 (MQ=255)

                                tGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGtt                                                         >  1:8349940/1‑93 (MQ=255)

                                tGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGtt                                                         >  2:9229671/1‑93 (MQ=255)

                                 gAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTg                                                        >  1:283463/1‑93 (MQ=255)

                                 gAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTg                                                        >  2:6886268/1‑93 (MQ=255)

                                 gAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTg                                                        <  1:8499493/93‑1 (MQ=255)

                                 gAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTg                                                        >  1:3869713/1‑93 (MQ=255)

                                  aaaGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTgg                                                       <  1:6461339/93‑1 (MQ=255)

                                   aaGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCACGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGga                                                      <  2:8720509/93‑1 (MQ=255)

                                    aGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGgat                                                     >  1:1110374/1‑93 (MQ=255)

                                      tCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGgatga                                                   <  2:402680/93‑1 (MQ=255)

                                        gAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATa                                                 <  1:8134007/93‑1 (MQ=255)

                                        gAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATa                                                 >  2:540113/1‑93 (MQ=255)

                                         aaaaCGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATaa                                                >  2:5619748/1‑92 (MQ=255)

                                          aaaCGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATaga                                               >  2:4274239/1‑93 (MQ=255)

                                          aaaCGTCCAGAAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGGTGGATGATaga                                               >  2:3805236/1‑93 (MQ=255)

                                            aCGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGACGACGTGAAAGAGATCCGTTGGATGATagaag                                             >  2:9830583/1‑93 (MQ=255)

                                            aCGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATagaag                                             >  2:482521/1‑93 (MQ=255)

                                            aCGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGCTGGATGATagaag                                             >  1:8953565/1‑93 (MQ=255)

                                              gTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAagag                                           >  1:5160335/1‑93 (MQ=255)

                                                ccAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGtt                                         <  2:2908744/93‑1 (MQ=255)

                                                ccAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGtt                                         <  2:2482536/93‑1 (MQ=255)

                                                ccAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGtt                                         <  1:572764/93‑1 (MQ=255)

                                                 cagcagGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTg                                        >  2:3530627/1‑93 (MQ=255)

                                                  agcagGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTgc                                       >  2:5356522/1‑93 (MQ=255)

                                                  agcagGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTgc                                       <  1:1818803/93‑1 (MQ=255)

                                                      ggCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCgt                                   >  2:8691351/1‑93 (MQ=255)

                                                      ggCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCgt                                   >  2:7700133/1‑93 (MQ=255)

                                                        cGTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCgtta                                 <  2:3766093/93‑1 (MQ=255)

                                                         gTGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAGAGAGATCCGTTGGATGATAGAAGAGTTGCGCgttag                                <  2:2709927/93‑1 (MQ=255)

                                                          tGGCAGCAATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCgttagt                               <  2:4754868/93‑1 (MQ=255)

                                                            gcagcaATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCgttagtta                             <  2:1402063/93‑1 (MQ=255)

                                                              agcaATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACt                           >  1:6058286/1‑93 (MQ=255)

                                                               gcaATGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACtt                          >  2:9685084/1‑93 (MQ=255)

                                                                 aaTGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCg                        >  1:2231462/1‑93 (MQ=255)

                                                                 aaTGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCg                        >  2:3917432/1‑93 (MQ=255)

                                                                 aaTGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCg                        >  1:6798486/1‑93 (MQ=255)

                                                                  aTGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGc                       >  2:5105698/1‑93 (MQ=255)

                                                                   tGGATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCt                      >  1:7504557/1‑93 (MQ=255)

                                                                     gATCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTca                    <  2:2975949/93‑1 (MQ=255)

                                                                       tCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTcaac                  >  2:9425001/1‑93 (MQ=255)

                                                                       tCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTcaac                  >  2:2060822/1‑93 (MQ=255)

                                                                       tCAACAAACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTcaac                  <  1:259951/93‑1 (MQ=255)

                                                                          acaaACTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCaacaac               <  1:9395977/93‑1 (MQ=255)

                                                                            aaaCTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCAACAACtg             >  1:1545582/1‑92 (MQ=255)

                                                                             aaCTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCAACAACTTg            <  2:8772243/93‑1 (MQ=255)

                                                                              aCTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCAACAACTTgg           <  2:4412009/93‑1 (MQ=255)

                                                                               cTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCAACAACTTGGt          >  2:7219322/1‑93 (MQ=255)

                                                                               cTGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCAACAACTTGGt          >  1:10044158/1‑93 (MQ=255)

                                                                                tGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCAACAACTTGGTa         <  2:5445456/93‑1 (MQ=255)

                                                                                tGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCAACAACTTGGTa         <  2:9733925/93‑1 (MQ=255)

                                                                                tGCCGCCCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCAACAACTTGGTa         <  1:4116202/93‑1 (MQ=255)

                                                                                 gccgccCGCGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCAACAACTTGGTAc        <  1:5351690/93‑1 (MQ=255)

                                                                                    gcccgcGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCAACAACTTGGTACGcc     >  2:9892469/1‑93 (MQ=255)

                                                                                     cccgcGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCACCAACTAGCTACGCCt    >  2:1073371/1‑93 (MQ=255)

                                                                                     cccgcGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCAACAACTTGGTACGCCt    >  2:1764328/1‑93 (MQ=255)

                                                                                     cccgcGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCAACAACTTGGTACGCCt    >  2:7558816/1‑93 (MQ=255)

                                                                                      ccgcGCGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCAACAACTTGGTACGCCtt   >  2:32400/1‑93 (MQ=255)

                                                                                       cgcgcgTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCAACAACTTGGTACGCCTTa  >  2:1099390/1‑93 (MQ=255)

                                                                                       cgcgcgTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCAACAACTTGGTACGCCTTa  >  2:5005547/1‑93 (MQ=255)

GGTTGATCCACATCGCGACCGTGCGCAGATGCTGAAAGTCGAAAACGTCCAGCAGGCGTGGCAGCAATGGATCAACAAACTGCCGCCCGCACGTCGTGAGGATGAAGACGTGAAAGAGATCCGTTGGATGATAGAAGAGTTGCGCGTTAGTTACTTCGCTCAACAACTTGGTACGCCTTA  >  NC_000913/1484756‑1484935

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