BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,896,324	+G	coding (698/711 nt)	b3718 ←	putative isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,896,324	1	.	G	99.6%	1652.5 / ‑3.1	509	coding (698/711 nt)	b3718	putative isomerase
Reads supporting (aligned to +/- strand): ref base . (2/0); new base G (213/294); total (215/294)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 562 total aligned reads displayed.

TTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CA‑CGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTCGGTAACCGGACCTTGCAGAACA  >  NC_000913/3896232‑3896410

  cGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑Cagc                                                                                        <  2:4800256/93‑3 (MQ=255)

  cGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑Cagc                                                                                        >  1:3514901/1‑91 (MQ=255)

    tAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑Cagcgg                                                                                      <  2:3552607/93‑5 (MQ=255)

    tAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑Cagcgg                                                                                      >  1:8079872/1‑89 (MQ=255)

    tAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAAACGAGTGCTAATTCTG‑CGG‑Cagcgg                                                                                      >  2:265404/1‑89 (MQ=255)

      gTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCt                                                                                    <  1:7441880/93‑1 (MQ=255)

       ttATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCtt                                                                                   >  2:8654144/1‑93 (MQ=255)

       ttATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCtt                                                                                   <  2:4162095/93‑1 (MQ=255)

       ttATATAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCtt                                                                                   <  1:9097203/93‑1 (MQ=255)

         aTTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTa                                                                                 <  1:1489755/93‑1 (MQ=255)

         aTTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTa                                                                                 <  2:3206493/93‑1 (MQ=255)

         aTTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTa                                                                                 <  2:8146795/93‑1 (MQ=255)

           ttAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATc                                                                               <  2:9043315/93‑1 (MQ=255)

            tAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCa                                                                              <  2:4331973/93‑1 (MQ=255)

             aaaTAATGAGAACAGGCCGGAGCGTAATTCACACAGCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAg                                                                             <  2:5148172/93‑1 (MQ=255)

              aataatGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGc                                                                            <  2:6563257/93‑1 (MQ=255)

              aataatGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGc                                                                            >  2:442000/1‑93 (MQ=255)

               ataatGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATGCTG‑CGG‑CAGCGGCTTTATCAGCg                                                                           >  2:5550649/1‑93 (MQ=255)

                taatGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CTG‑CAGCGGCTTTATCAGcta                                                                          >  1:4311917/1‑91 (MQ=255)

                taatGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGa                                                                          <  1:8624656/93‑1 (MQ=255)

                 aatGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGAt                                                                         >  1:3285608/1‑93 (MQ=255)

                 aatGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGAt                                                                         >  2:4315902/1‑93 (MQ=255)

                 aatGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGAt                                                                         <  2:3874054/93‑1 (MQ=255)

                  atGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGAtt                                                                        >  1:7798963/1‑93 (MQ=255)

                  atGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGAtt                                                                        <  2:1392333/93‑1 (MQ=255)

                  atGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGAtt                                                                        >  2:333907/1‑93 (MQ=255)

                   tGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTa                                                                       <  1:5095525/93‑1 (MQ=255)

                    gagaACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTAc                                                                      <  1:4397854/93‑1 (MQ=255)

                      gaACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCa                                                                    <  2:3090751/93‑1 (MQ=255)

                       aaCAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCAt                                                                   <  2:6967023/93‑1 (MQ=255)

                        aCAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATc                                                                  <  1:5366807/93‑1 (MQ=255)

                        aCAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATc                                                                  <  1:6042973/93‑1 (MQ=255)

                         cAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCa                                                                 >  2:1732078/1‑93 (MQ=255)

                         cAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCa                                                                 <  2:10026636/93‑1 (MQ=255)

                              cGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGaa                                                            >  1:6363066/1‑93 (MQ=255)

                              cGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGaa                                                            <  1:2568155/93‑1 (MQ=255)

                               ggAGCGTAATTCACACATCCGGCCTTATTTCTTAGGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAg                                                           >  1:6687225/1‑93 (MQ=255)

                                 aGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGa                                                         <  2:6078434/93‑1 (MQ=255)

                                 aGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGa                                                         >  2:6382276/1‑93 (MQ=255)

                                   cGTAATTCACACGTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATg                                                       <  1:9578330/93‑1 (MQ=255)

                                    gTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGc                                                      <  1:943764/93‑1 (MQ=255)

                                    gTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGc                                                      <  1:9751969/93‑1 (MQ=255)

                                        ttCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCt                                                  <  1:7120398/93‑1 (MQ=255)

                                         tCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTg                                                 >  2:2588291/1‑93 (MQ=255)

                                          cacacaTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGc                                                <  2:9353459/93‑1 (MQ=255)

                                           acacaTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCa                                               <  1:367607/93‑1 (MQ=255)

                                             acaTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAac                                             >  1:9932232/1‑93 (MQ=255)

                                              caTCCGGGCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAaca                                            <  1:1850115/93‑1 (MQ=255)

                                              caTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAaca                                            <  1:8380046/93‑1 (MQ=255)

                                               aTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAacac                                           >  1:7980915/1‑93 (MQ=255)

                                               aTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAacac                                           >  1:2012997/1‑93 (MQ=255)

                                                tCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACg                                          <  2:7590538/93‑1 (MQ=255)

                                                tCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACg                                          <  1:1202881/93‑1 (MQ=255)

                                                   ggCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGaa                                       >  1:9224219/1‑93 (MQ=255)

                                                   ggCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGaa                                       <  2:5887591/93‑1 (MQ=255)

                                                   ggCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGaa                                       <  2:7331885/93‑1 (MQ=255)

                                                    gCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAg                                      >  2:1150126/1‑93 (MQ=255)

                                                    gCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAg                                      <  1:575392/93‑1 (MQ=255)

                                                    gCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAg                                      <  1:2369058/93‑1 (MQ=255)

                                                     ccTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGc                                     <  2:16075/93‑1 (MQ=255)

                                                       ttATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTg                                   <  1:8873557/93‑1 (MQ=255)

                                                       ttATTTCCTAAGCTTAACCGAGTGCTAATTCCG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTg                                   <  2:2026418/93‑1 (MQ=255)

                                                        tATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGACGATgg                                  >  1:3811106/1‑93 (MQ=255)

                                                        tATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTgg                                  >  2:8585124/1‑93 (MQ=255)

                                                        tATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTgg                                  >  2:8879054/1‑93 (MQ=255)

                                                        tATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTgg                                  <  1:8505310/93‑1 (MQ=255)

                                                          tttCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGaa                                <  1:8175883/93‑1 (MQ=255)

                                                          tttCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGaa                                >  1:4148913/1‑93 (MQ=255)

                                                          tttCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGaa                                >  1:5193631/1‑93 (MQ=255)

                                                          tttCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGaa                                <  1:4941579/93‑1 (MQ=255)

                                                           ttCTTAAGCTTAACCGAGTGCTAATTCTGAGGG‑CG‑CGGGTTTATCAGAGATTACCCTCAGAGCAGGGTCGAGATGCAACCCGGCAGCTGGaag                               >  1:695600/1‑92 (MQ=255)

                                                             cTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑t                            >  2:3681978/1‑93 (MQ=255)

                                                              ttAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGGAGCTGGAACA‑Tc                           >  2:6904392/1‑93 (MQ=255)

                                                              ttAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑Tc                           >  2:9602480/1‑93 (MQ=255)

                                                              ttAAGCTTAACCGAGTGCTAATTCCG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑Tc                           <  1:6005910/93‑1 (MQ=255)

                                                               tAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCt                          <  1:485427/93‑1 (MQ=255)

                                                               tAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCt                          <  1:3650938/93‑1 (MQ=255)

                                                                aaGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCtt                         <  1:10028388/93‑1 (MQ=255)

                                                                aaGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCtt                         <  2:798985/93‑1 (MQ=255)

                                                                 aGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTc                        <  2:5262700/93‑1 (MQ=255)

                                                                 aGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTc                        <  2:826934/93‑1 (MQ=255)

                                                                 aGCTTAACCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTc                        <  1:2879352/93‑1 (MQ=255)

                                                                  gCTTAACCGACTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACAAGGAAGCTGGAACA‑TCTTCg                       >  2:7142878/1‑93 (MQ=255)

                                                                       aCCGAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTCGGTAAc                  <  1:6813893/93‑1 (MQ=255)

                                                                        ccGAGTGCTAATTCTG‑CGGTCA‑TGGCTTTATCATCGATTACCATCCGCGAAGGCCGCACCTGCAACACCGACG‑TGGGACACTCTTCGGTATcc                 >  2:8512679/1‑93 (MQ=255)

                                                                          gAGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTCGGTAACCgg               >  2:397945/1‑93 (MQ=255)

                                                                           aGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTCGGTAACCGGa              >  2:2320718/1‑93 (MQ=255)

                                                                           aGTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTCGGTAACCGGa              >  1:4601492/1‑93 (MQ=255)

                                                                            gTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTCGGTAACCGGAc             <  1:732226/93‑1 (MQ=255)

                                                                            gTGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTCGGTAACCGGAc             <  2:4549488/93‑1 (MQ=255)

                                                                             tGCTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTCGGTAACCGGAcc            <  2:9213927/93‑1 (MQ=255)

                                                                               cTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTCGGTAACCGGACCtt          >  1:3990793/1‑93 (MQ=255)

                                                                               cTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTCGGTAACCGGACCtt          <  1:5526967/93‑1 (MQ=255)

                                                                               cTAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTCGGTAACCGGACCtt          <  2:2966321/93‑1 (MQ=255)

                                                                                tAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAATA‑TCTTCGGTAACCGGACCTTg         >  2:5719323/1‑93 (MQ=255)

                                                                                tAATTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTCGGTAACCGGACCTTg         <  1:3163078/93‑1 (MQ=255)

                                                                                 aaTTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTCGGTAACCGGACCTTGc        <  1:7650313/93‑1 (MQ=255)

                                                                                  aTTCTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTCGGTAACCGGACCTTGCa       <  1:7180389/93‑1 (MQ=255)

                                                                                     cTG‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTCGGTAACCGGACCTTGCAGaa    <  1:2631532/93‑1 (MQ=255)

                                                                                       g‑CGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTCGGTAACCGGACCTTGCAGAACa  <  1:1330527/93‑1 (MQ=255)

TTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTG‑CGG‑CA‑CGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACA‑TCTTCGGTAACCGGACCTTGCAGAACA  >  NC_000913/3896232‑3896410

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