BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,830,970	T→C	R173R (CGT→CGC)	xthA →	exonuclease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,830,970	0	T	C	99.5%	749.4 / ‑5.6	212	R173R (CGT→CGC)	xthA	exonuclease III
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (97/114); total (98/114)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 219 total aligned reads displayed.

TAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGTACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAA  >  NC_000913/1830880‑1831059

tAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTCAGCGCTGGCTGCGCAc                                                                                         >  1:4914052/1‑93 (MQ=255)

tAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCAc                                                                                         <  2:2718535/93‑1 (MQ=255)

 aaTCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCAcc                                                                                        <  1:5984125/93‑1 (MQ=255)

   tCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCgg                                                                                      >  2:2172359/1‑93 (MQ=255)

    ccGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGt                                                                                     <  1:2205188/93‑1 (MQ=255)

     cGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTa                                                                                    <  1:53892/93‑1 (MQ=255)

      ggTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTaa                                                                                   >  1:106840/1‑93 (MQ=255)

       gTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTaaa                                                                                  <  1:200896/93‑1 (MQ=255)

        tACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAAt                                                                                 <  2:4626636/93‑1 (MQ=255)

        tACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAAt                                                                                 <  1:5187373/93‑1 (MQ=255)

        tACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGCAAAt                                                                                 >  1:4405836/1‑93 (MQ=255)

           tGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGct                                                                              <  1:1573699/93‑1 (MQ=255)

           tGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGct                                                                              >  2:5977691/1‑93 (MQ=255)

           tGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGct                                                                              >  2:1892156/1‑93 (MQ=255)

            gATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGctt                                                                             >  2:4911619/1‑92 (MQ=255)

            gATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGctc                                                                             >  1:3950330/1‑93 (MQ=255)

              ttATGGGCGATATGAATATCAGCCCTACAGATCTGTATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCtt                                                                           >  1:96211/1‑93 (MQ=255)

              ttATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCtt                                                                           >  2:1141146/1‑93 (MQ=255)

              ttATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCtt                                                                           >  2:167010/1‑93 (MQ=255)

              ttATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCtt                                                                           >  2:1397/1‑93 (MQ=255)

               tATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCttt                                                                          >  1:5757880/1‑93 (MQ=255)

                aTGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTc                                                                         <  2:1982930/93‑1 (MQ=255)

                 tGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTcc                                                                        <  2:4559047/93‑1 (MQ=255)

                 tGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTcc                                                                        >  1:487045/1‑93 (MQ=255)

                 tGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTcc                                                                        <  2:1611076/93‑1 (MQ=255)

                  gggCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCt                                                                       <  2:488884/93‑1 (MQ=255)

                  gggCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCt                                                                       <  2:3986409/93‑1 (MQ=255)

                    gCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGc                                                                     >  2:4389302/1‑93 (MQ=255)

                     cGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGcc                                                                    <  1:6088747/93‑1 (MQ=255)

                      gATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCg                                                                   <  1:931956/93‑1 (MQ=255)

                      gATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCg                                                                   >  1:4346804/1‑93 (MQ=255)

                        tatGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGa                                                                 >  2:183253/1‑93 (MQ=255)

                        tatGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGa                                                                 >  1:3393964/1‑93 (MQ=255)

                        tatGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGa                                                                 >  1:29155/1‑93 (MQ=255)

                         atGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGaa                                                                >  1:5727783/1‑93 (MQ=255)

                          tGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAag                                                               >  2:1076965/1‑93 (MQ=255)

                          tGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAag                                                               <  2:3880538/93‑1 (MQ=255)

                           gAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAaga                                                              >  1:5658000/1‑93 (MQ=255)

                            aaTATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGCagag                                                             >  2:71429/1‑93 (MQ=255)

                               atCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAgcgc                                                          <  2:5731928/93‑1 (MQ=255)

                                tCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGcgcg                                                         >  1:4525466/1‑93 (MQ=255)

                                 cAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGa                                                        <  2:853476/93‑1 (MQ=255)

                                   gCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGATCTTTCCTGCCGGAAGAGCGCGAat                                                      >  1:2812476/1‑93 (MQ=255)

                                     ccTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAatgg                                                    <  2:880552/93‑1 (MQ=255)

                                      cTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAatgga                                                   >  1:3805081/1‑93 (MQ=255)

                                       tACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAatggat                                                  >  2:613926/1‑93 (MQ=255)

                                       tACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAatggat                                                  <  1:1414717/93‑1 (MQ=255)

                                       tACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAatggat                                                  >  2:3969597/1‑93 (MQ=255)

                                       tACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAatggat                                                  <  2:2971445/93‑1 (MQ=255)

                                       tACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAatggat                                                  <  1:2513227/93‑1 (MQ=255)

                                       tACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAatggat                                                  >  2:264451/1‑93 (MQ=255)

                                       tACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAatggat                                                  >  1:2565140/1‑93 (MQ=255)

                                        aCAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAatggatg                                                 >  1:4957737/1‑93 (MQ=255)

                                            aTCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACa                                             >  2:433892/1‑93 (MQ=255)

                                            aTCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACa                                             <  1:459431/93‑1 (MQ=255)

                                             tCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAg                                            <  2:546030/93‑1 (MQ=255)

                                             tCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAg                                            >  2:1323496/1‑93 (MQ=255)

                                               tGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGc                                          <  1:5078190/93‑1 (MQ=255)

                                                ggATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCt                                         >  1:128014/1‑93 (MQ=255)

                                                   tatCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCtgat                                      >  1:317101/1‑93 (MQ=255)

                                                   tatCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCtgat                                      >  1:3102020/1‑93 (MQ=255)

                                                   tatCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCtgat                                      >  2:468330/1‑93 (MQ=255)

                                                   tatCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCtgat                                      <  2:4523788/93‑1 (MQ=255)

                                                    atCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCtgatg                                     >  1:1156140/1‑93 (MQ=255)

                                                    atCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCtgatg                                     <  1:747739/93‑1 (MQ=255)

                                                     tCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCtgatga                                    <  1:1534801/93‑1 (MQ=255)

                                                       ggCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGc                                  <  1:2986831/93‑1 (MQ=255)

                                                        gCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCt                                 >  1:3912423/1‑93 (MQ=255)

                                                         cATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTg                                <  2:3517410/93‑1 (MQ=255)

                                                          aTTGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTgg                               <  1:4633295/93‑1 (MQ=255)

                                                           ttGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTggg                              >  1:1029952/1‑93 (MQ=255)

                                                           ttGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGAAGAGCgggg                              >  1:2637437/1‑93 (MQ=255)

                                                            tGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTgggg                             <  1:351671/93‑1 (MQ=255)

                                                            tGGCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTgggg                             <  1:1371059/93‑1 (MQ=255)

                                                              gCGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGt                           >  2:5074570/1‑93 (MQ=255)

                                                               cGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGtt                          <  1:1121897/93‑1 (MQ=255)

                                                               cGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGtt                          <  1:5904094/93‑1 (MQ=255)

                                                               cGAAGAAAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGtt                          <  2:1864441/93‑1 (MQ=255)

                                                                gaagaaAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTg                         <  2:5130274/93‑1 (MQ=255)

                                                                gaagaaAACCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGGTg                         >  1:4181375/1‑93 (MQ=255)

                                                                      acccGTAAGCGCGGGCTGCACACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGAt                   <  1:5465324/91‑1 (MQ=255)

                                                                      aaCCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGAt                   >  2:1575310/1‑93 (MQ=255)

                                                                      aaCCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGAt                   <  2:173938/93‑1 (MQ=255)

                                                                      aaCCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGAt                   <  1:1298122/93‑1 (MQ=255)

                                                                      aaCCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCACGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGAt                   >  1:174619/1‑93 (MQ=255)

                                                                       aCCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATa                  <  1:1524009/93‑1 (MQ=255)

                                                                       aCCGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATa                  <  2:5339388/93‑1 (MQ=255)

                                                                         cGTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATAcc                >  2:562099/1‑93 (MQ=255)

                                                                          ttAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGGCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCt               <  1:81370/92‑1 (MQ=255)

                                                                          gTAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCt               >  2:224377/1‑93 (MQ=255)

                                                                          ctAAGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCt               <  1:4793164/92‑1 (MQ=255)

                                                                            aaGCGCTGGCTGCGTACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTc             >  2:5191252/1‑93 (MQ=255)

                                                                            aaGCGCTGGCTGCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTc             <  1:2376897/93‑1 (MQ=255)

                                                                             aGCGCTGGCTGCGCACAGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTcc            <  2:4746797/93‑1 (MQ=255)

                                                                                 ctggctgCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCa        <  2:2416117/93‑1 (MQ=255)

                                                                                  tggctgCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCAt       <  2:1218411/93‑1 (MQ=255)

                                                                                  tggctgCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCAt       >  1:1818344/1‑93 (MQ=255)

                                                                                  tggctgCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCAt       <  2:3741356/93‑1 (MQ=255)

                                                                                  tggctgCGCACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCAt       <  2:2904429/93‑1 (MQ=255)

                                                                                       gcgcACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGaa  <  2:3187442/93‑1 (MQ=255)

TAATCCGGTACTGATTATGGGCGATATGAATATCAGCCCTACAGATCTGGATATCGGCATTGGCGAAGAAAACCGTAAGCGCTGGCTGCGTACCGGTAAATGCTCTTTCCTGCCGGAAGAGCGCGAATGGATGGACAGGCTGATGAGCTGGGGGTTGGTCGATACCTTCCGCCATGCGAA  >  NC_000913/1830880‑1831059

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