BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,106,960	T→C	D216G (GAC→GGC)	rfbX ←	putative O‑antigen transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,106,960	0	T	C	99.5%	699.7 / NA	218	D216G (GAC→GGC)	rfbX	putative O‑antigen transporter
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (118/98); total (119/99)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

ACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCCA  >  NC_000913/2106871‑2107049

 cGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGc                                                                                       >  2:5180822/1‑93 (MQ=255)

 cGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGc                                                                                       <  1:3181282/93‑1 (MQ=255)

  gTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCt                                                                                      >  1:2040850/1‑93 (MQ=255)

   ttGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTa                                                                                     >  2:4146125/1‑93 (MQ=255)

   ttGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTa                                                                                     >  1:4721516/1‑93 (MQ=255)

      gTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAg                                                                                  >  1:1338513/1‑93 (MQ=255)

        ccAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGAt                                                                                <  1:919968/93‑1 (MQ=255)

        ccAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGAt                                                                                <  1:3715041/93‑1 (MQ=255)

         cAGATATAAATCCCAGGATAATAGCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAATCCGCCTGCTAAAGATc                                                                               <  2:5737366/93‑1 (MQ=255)

         cAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATc                                                                               <  1:1016862/93‑1 (MQ=255)

         cAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATc                                                                               <  1:5318593/93‑1 (MQ=255)

           gATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGa                                                                             <  1:4105390/93‑1 (MQ=255)

           gATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGa                                                                             <  1:2731750/93‑1 (MQ=255)

           gATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGa                                                                             >  1:279143/1‑93 (MQ=255)

           gATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGa                                                                             <  1:6072735/93‑1 (MQ=255)

            atatAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGAc                                                                            <  2:3188707/93‑1 (MQ=255)

              ataAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCAACGc                                                                          >  2:209694/1‑93 (MQ=255)

               taAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGca                                                                         <  2:126596/93‑1 (MQ=255)

               taAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGca                                                                         >  1:4693165/1‑93 (MQ=255)

               taAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGca                                                                         <  1:4786065/93‑1 (MQ=255)

                aaaTCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGcac                                                                        >  1:2396957/1‑93 (MQ=255)

                aaaTCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGcac                                                                        >  2:1900463/1‑93 (MQ=255)

                 aaTCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGcaca                                                                       >  1:3363758/1‑93 (MQ=255)

                  aTCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCTCAt                                                                      <  2:2911031/93‑1 (MQ=255)

                  aTCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACAt                                                                      <  2:5775073/93‑1 (MQ=255)

                   tCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACAtt                                                                     >  2:4399103/1‑93 (MQ=255)

                      cAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATg                                                                  <  2:4366782/93‑1 (MQ=255)

                       aGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGt                                                                 >  2:6295813/1‑93 (MQ=255)

                        ggATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTa                                                                <  1:2606478/93‑1 (MQ=255)

                         gATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTaa                                                               >  1:801715/1‑93 (MQ=255)

                         gATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTaa                                                               <  2:5580819/93‑1 (MQ=255)

                          ataataACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAAt                                                              >  2:5502168/1‑93 (MQ=255)

                          ataataACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAAt                                                              >  2:5210170/1‑93 (MQ=255)

                           taataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATg                                                             >  2:1602135/1‑93 (MQ=255)

                            aataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGa                                                            >  2:6032539/1‑93 (MQ=255)

                            aataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGa                                                            >  2:5292249/1‑93 (MQ=255)

                             ataaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGAt                                                           <  1:1544136/93‑1 (MQ=255)

                              taaCTATTCCCGTAGAATATCAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATa                                                          >  1:148031/1‑93 (MQ=255)

                              taaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATa                                                          >  2:4569759/1‑93 (MQ=255)

                              taaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATa                                                          <  1:4985440/93‑1 (MQ=255)

                               aaCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAg                                                         <  1:3921281/93‑1 (MQ=255)

                                aCTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGc                                                        <  1:4313180/93‑1 (MQ=255)

                                 cTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCg                                                       >  2:2589626/1‑93 (MQ=255)

                                 cTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCg                                                       <  1:1830976/93‑1 (MQ=255)

                                 cTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCg                                                       <  2:4862352/93‑1 (MQ=255)

                                  tATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGt                                                      <  1:2381195/93‑1 (MQ=255)

                                      cccGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACt                                                  <  1:3596988/93‑1 (MQ=255)

                                       ccGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACtt                                                 <  1:3390315/93‑1 (MQ=255)

                                       ccGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACtt                                                 <  1:3473079/93‑1 (MQ=255)

                                       ccGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACtt                                                 >  2:5984087/1‑93 (MQ=255)

                                        cGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACttt                                                >  2:4273583/1‑93 (MQ=255)

                                        cGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACttt                                                >  2:2387547/1‑93 (MQ=255)

                                        cGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACttt                                                >  2:883710/1‑93 (MQ=255)

                                         gTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTa                                               >  2:1383818/1‑93 (MQ=255)

                                            gAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCa                                            >  2:3678390/1‑93 (MQ=255)

                                             aaTATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCaa                                           <  2:403474/93‑1 (MQ=255)

                                               tataAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATc                                         >  1:5549807/1‑93 (MQ=255)

                                                ataAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATcc                                        >  2:5679727/1‑93 (MQ=255)

                                                 taAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCa                                       >  1:5896837/1‑93 (MQ=255)

                                                 taAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCa                                       <  2:852643/93‑1 (MQ=255)

                                                  aaaCTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCaa                                      >  1:2952015/1‑93 (MQ=255)

                                                  aaaCTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCaa                                      >  2:2958031/1‑93 (MQ=255)

                                                    aCTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAAcc                                    >  1:1789740/1‑93 (MQ=255)

                                                     cTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCt                                   <  1:719263/93‑1 (MQ=255)

                                                     cTAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCt                                   <  1:1856174/93‑1 (MQ=255)

                                                      tAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCtt                                  >  2:2060986/1‑93 (MQ=255)

                                                      tAATAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCtt                                  >  2:3074061/1‑93 (MQ=255)

                                                       aaTAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTc                                 >  2:103595/1‑93 (MQ=255)

                                                       aaTAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTc                                 >  2:1537234/1‑93 (MQ=255)

                                                        aTAGCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCa                                >  1:5628769/1‑93 (MQ=255)

                                                           gCAGATGTGGAAATAAAAACATGAAACCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGa                             <  1:4012312/93‑1 (MQ=255)

                                                           gCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGa                             >  1:4550878/1‑93 (MQ=255)

                                                           gCAGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGa                             >  2:3759370/1‑93 (MQ=255)

                                                             aGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaac                           >  2:1125844/1‑93 (MQ=255)

                                                             aGATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaac                           <  2:329143/93‑1 (MQ=255)

                                                              gATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCACCCTTCATGcaca                          >  1:485791/1‑93 (MQ=255)

                                                              gATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaaca                          <  2:479279/93‑1 (MQ=255)

                                                              gATGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaaca                          <  1:187462/93‑1 (MQ=255)

                                                               aTGTGGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGaacaa                         <  1:1248926/93‑1 (MQ=255)

                                                                 gtgGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACg                       >  1:1766013/1‑93 (MQ=255)

                                                                  tgGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGa                      >  1:4642490/1‑93 (MQ=255)

                                                                  tgGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGa                      <  2:3934831/93‑1 (MQ=255)

                                                                  tgGAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGa                      <  1:2109196/93‑1 (MQ=255)

                                                                   ggAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGAt                     >  2:2739061/1‑93 (MQ=255)

                                                                   ggAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGAt                     >  1:4466866/1‑93 (MQ=255)

                                                                   ggAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGAt                     <  1:4165866/93‑1 (MQ=255)

                                                                   ggAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGAt                     >  2:1665437/1‑93 (MQ=255)

                                                                    gAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATa                    >  1:5587011/1‑93 (MQ=255)

                                                                    gAAATAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATa                    <  1:2295695/93‑1 (MQ=255)

                                                                     aaaTAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAg                   <  2:4955492/93‑1 (MQ=255)

                                                                        tAAAAACATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTa                >  2:3891468/1‑93 (MQ=255)

                                                                          aaaaCATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTGCCAATCCAACCTTCATGAACAACGATAGCTAGt              <  2:4699852/93‑1 (MQ=255)

                                                                            aaCATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGc            <  1:1111373/93‑1 (MQ=255)

                                                                            aaCATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGc            <  1:5136606/93‑1 (MQ=255)

                                                                            aaCATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGc            <  1:5240368/93‑1 (MQ=255)

                                                                             agaTGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCa           >  2:2560975/3‑93 (MQ=255)

                                                                              cATGAAAACCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCaa          >  2:1337339/1‑93 (MQ=255)

                                                                                    aaCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTc    <  1:32091/93‑1 (MQ=255)

                                                                                    aaCCGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTc    >  2:5075840/1‑93 (MQ=255)

                                                                                      ccGCCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCCa  >  2:1299925/1‑93 (MQ=255)

ACGTTGGTCCAGATATAAATCCCAGGATAATAACTATTCCCGTAGAATATAAACTAATAGCAGATGTGGAAATAAAAACATGAAAACCGTCTGCTAAAGATCGACGCACATTATGTAATGATAGCGTAACTTTACCAATCCAACCTTCATGAACAACGATAGCTAGTGCAATAATTCCA  >  NC_000913/2106871‑2107049

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