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.5%	701.2 / ‑3.9	201	coding (698/711 nt)	b3718	putative isomerase
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (111/89); total (112/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 218 total aligned reads displayed.

TTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CA‑CGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACA  >  NC_000913/3896232‑3896410

ttCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑Ca                                                                                         >  1:4895541/1‑93 (MQ=255)

 tCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑Cag                                                                                        >  1:3745153/1‑92 (MQ=255)

 tCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑Cag                                                                                        <  1:2451988/93‑2 (MQ=255)

  cGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑Cagc                                                                                       >  2:2224830/1‑91 (MQ=255)

       ttATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCtt                                                                                  >  1:661731/1‑93 (MQ=255)

         aTTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTa                                                                                <  1:6161230/93‑1 (MQ=255)

          tttAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTAt                                                                               >  2:3291644/1‑93 (MQ=255)

          tttAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTAt                                                                               <  2:5059261/93‑1 (MQ=255)

           ttAAATAATGAGAACAGGCCGTAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATc                                                                              <  2:2584450/93‑1 (MQ=255)

           ttAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATc                                                                              <  1:2184557/93‑1 (MQ=255)

           ttAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATc                                                                              <  1:1480702/93‑1 (MQ=255)

             aaaTAATGAGAACAGGCCGGTGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAg                                                                            >  1:5251616/1‑93 (MQ=255)

             aaaTAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAg                                                                            >  1:2865640/1‑93 (MQ=255)

             aaaTAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAg                                                                            >  2:3210736/1‑93 (MQ=255)

              aataatGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGc                                                                           >  2:1828949/1‑93 (MQ=255)

              aataatGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGc                                                                           <  1:3086295/93‑1 (MQ=255)

               ataatGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCg                                                                          >  1:4138574/1‑93 (MQ=255)

                taatGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGa                                                                         <  2:3777478/93‑1 (MQ=255)

                 aatGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGAt                                                                        <  2:2690557/93‑1 (MQ=255)

                  atGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGAtt                                                                       <  2:452896/93‑1 (MQ=255)

                     agaACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTAcc                                                                    >  1:1946671/1‑93 (MQ=255)

                      gaACAGGCCGGAGCGTAATTCACGCATCCGGCCTTATTTCTTAAGCTTAGCCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCa                                                                   <  2:6328772/93‑1 (MQ=255)

                      gaACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCa                                                                   >  1:1454708/1‑93 (MQ=255)

                      gaACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCa                                                                   <  1:5738585/93‑1 (MQ=255)

                      gaACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCa                                                                   >  1:3036770/1‑93 (MQ=255)

                       aaCAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCAt                                                                  >  1:3540820/1‑93 (MQ=255)

                       aaCAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCAt                                                                  >  2:2833347/1‑93 (MQ=255)

                        aCAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATc                                                                 >  1:3831996/1‑93 (MQ=255)

                        aCAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATc                                                                 <  1:4558142/93‑1 (MQ=255)

                         cAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCa                                                                <  2:2097905/93‑1 (MQ=255)

                          aGGCCGGAGCGTAATTCCCACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAg                                                               <  1:3700437/93‑1 (MQ=255)

                           ggCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGc                                                              <  2:1512520/93‑1 (MQ=255)

                            gCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCg                                                             >  1:83253/1‑93 (MQ=255)

                            gCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCg                                                             <  2:5769386/93‑1 (MQ=255)

                              cGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGaa                                                           >  2:657608/1‑93 (MQ=255)

                              cGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGaa                                                           <  1:564604/93‑1 (MQ=255)

                               ggAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAg                                                          >  1:2362057/1‑93 (MQ=255)

                               ggAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAg                                                          <  2:3723402/93‑1 (MQ=255)

                               ggAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAg                                                          <  1:5388123/93‑1 (MQ=255)

                               ggAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAg                                                          >  2:1365852/1‑93 (MQ=255)

                               ggAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAg                                                          >  1:3273616/1‑93 (MQ=255)

                                 aGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGa                                                        <  2:4073089/93‑1 (MQ=255)

                                  gCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGAt                                                       <  2:4391184/93‑1 (MQ=255)

                                  gCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGAt                                                       <  1:337116/93‑1 (MQ=255)

                                    gTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGc                                                     >  2:4279295/1‑93 (MQ=255)

                                    gTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGc                                                     <  1:4479522/93‑1 (MQ=255)

                                     tAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCa                                                    <  2:817695/93‑1 (MQ=255)

                                      aaTTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAg                                                   <  2:3558641/93‑1 (MQ=255)

                                      aaTTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAg                                                   <  2:1650113/93‑1 (MQ=255)

                                      aaTTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAg                                                   >  2:3865126/1‑93 (MQ=255)

                                          cacacaTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGc                                               >  2:5657956/1‑93 (MQ=255)

                                          cacacaTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGc                                               >  2:1463984/1‑93 (MQ=255)

                                            cacaTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCaa                                             <  1:2968760/93‑1 (MQ=255)

                                              caTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAaca                                           >  2:3111105/1‑93 (MQ=255)

                                              caTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAaca                                           <  1:25437/93‑1 (MQ=255)

                                                tCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGTGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACg                                         >  1:1167295/1‑93 (MQ=255)

                                                tCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACg                                         >  1:5897144/1‑93 (MQ=255)

                                                  cGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGa                                       >  1:3500122/1‑93 (MQ=255)

                                                    gCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAg                                     >  2:1198324/1‑93 (MQ=255)

                                                      cTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCt                                   >  2:3131756/1‑93 (MQ=255)

                                                       ttATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTg                                  >  1:5191078/1‑93 (MQ=255)

                                                       ttATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTg                                  <  2:4996049/93‑1 (MQ=255)

                                                        tATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTgg                                 >  2:3191665/1‑93 (MQ=255)

                                                         aTTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGa                                >  1:5047819/1‑93 (MQ=255)

                                                         aTTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGa                                >  2:4920510/1‑93 (MQ=255)

                                                          tttCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGaa                               <  2:5448851/93‑1 (MQ=255)

                                                           ttCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAAc                              >  2:1979528/1‑93 (MQ=255)

                                                            tCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACa                             >  1:6362004/1‑93 (MQ=255)

                                                            tCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACa                             >  1:3569702/1‑93 (MQ=255)

                                                             cTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACAt                            <  1:5055970/93‑1 (MQ=255)

                                                              ttAAGCTTAACCGAGTGCTAATTCTGCGGCCA‑CGGCTTTATCAGCGATTACCATCAGCGTAGGATGCAGCTGCAACACGGAAGCTGGAACATc                           >  1:59606/1‑93 (MQ=255)

                                                              ttAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATc                           >  2:6408426/1‑93 (MQ=255)

                                                              ttAAGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATc                           >  2:4266768/1‑93 (MQ=255)

                                                                 aGCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTc                        >  1:2707149/1‑93 (MQ=255)

                                                                 aGCTTAACCGAGTGCTAATCCTGCGG‑CAGCGGCTTTATCAGCGATTCCCATGAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTc                        <  2:1416235/93‑1 (MQ=255)

                                                                  gCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCg                       >  2:5197387/1‑93 (MQ=255)

                                                                  gCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCg                       <  2:1101643/93‑1 (MQ=255)

                                                                  gCTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCg                       <  2:3713096/93‑1 (MQ=255)

                                                                   cTTAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCgg                      >  1:1281804/1‑93 (MQ=255)

                                                                    ttAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGt                     >  1:4969050/1‑93 (MQ=255)

                                                                     tAACCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTa                    >  1:1676591/1‑93 (MQ=255)

                                                                       aCCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAAc                  >  1:5121099/1‑93 (MQ=255)

                                                                       aCCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAAc                  >  2:2321906/1‑93 (MQ=255)

                                                                       aCCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAAc                  >  1:991857/1‑93 (MQ=255)

                                                                       aCCGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAAc                  >  1:4294030/1‑93 (MQ=255)

                                                                        ccGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAAcc                 >  2:4903528/1‑93 (MQ=255)

                                                                        ccGAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAAcc                 >  1:4167048/1‑93 (MQ=255)

                                                                          gAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAAcctg               >  2:2382034/1‑91 (MQ=255)

                                                                          gAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCgg               >  2:4702412/1‑93 (MQ=255)

                                                                          gAGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCgg               >  2:4317368/1‑93 (MQ=255)

                                                                           aGTGCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGa              >  1:6172870/1‑93 (MQ=255)

                                                                              gCTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCt           <  2:6447081/93‑1 (MQ=255)

                                                                               cTAATTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCtt          <  1:4237376/93‑1 (MQ=255)

                                                                                 aaTTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGc        <  2:3450578/93‑1 (MQ=255)

                                                                                 aaTTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGc        >  2:1849323/1‑93 (MQ=255)

                                                                                 aaTTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGc        >  2:6209096/1‑93 (MQ=255)

                                                                                  aTTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCa       >  1:5280549/1‑93 (MQ=255)

                                                                                  aTTCTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCa       <  2:25193/93‑1 (MQ=255)

                                                                                     cTGCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGaa    >  1:927853/1‑93 (MQ=255)

                                                                                       gCGG‑CAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACa  >  2:5258686/1‑93 (MQ=255)

TTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGG‑CA‑CGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACA  >  NC_000913/3896232‑3896410

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