BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,823,984	A→G	Q383R (CAG→CGG)	recG →	DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,823,984	0	A	G	84.9%	767.8 / 132.4	332	Q383R (CAG→CGG)	recG	DNA helicase
Reads supporting (aligned to +/- strand): ref base A (23/27); new base G (119/163); total (142/190)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.59e-01

Only 100 of 336 total aligned reads displayed.

GGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTCTGGCATTGTGGGAGA  >  NC_000913/3823893‑3824073

 gTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCggg                                                                                         >  2:2881603/1‑93 (MQ=255)

 gTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCggg                                                                                         >  2:4423579/1‑93 (MQ=255)

  tAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGa                                                                                        <  1:4813531/93‑1 (MQ=255)

    aGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAAc                                                                                      <  1:5850783/93‑1 (MQ=255)

     gCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACa                                                                                     >  1:3497871/1‑93 (MQ=255)

     gCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACa                                                                                     <  2:6582083/93‑1 (MQ=255)

     gCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACa                                                                                     <  2:5588728/93‑1 (MQ=255)

     gCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACa                                                                                     >  2:4141699/1‑93 (MQ=255)

      cAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAg                                                                                    <  2:6847205/93‑1 (MQ=255)

       aGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAgg                                                                                   <  1:1544494/93‑1 (MQ=255)

       aGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAgg                                                                                   <  2:3670255/93‑1 (MQ=255)

        gAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGt                                                                                  >  1:4519788/1‑93 (MQ=255)

          aaGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGc                                                                                >  1:2395293/1‑93 (MQ=255)

            ggTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAg                                                                              <  2:5087319/93‑1 (MQ=255)

            ggTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAg                                                                              <  1:3260595/93‑1 (MQ=255)

             gTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGt                                                                             <  1:6030492/93‑1 (MQ=255)

             gTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGt                                                                             >  2:2172985/1‑93 (MQ=255)

              tAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGCGCAGtt                                                                            >  1:3209688/1‑93 (MQ=255)

                 aGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTaa                                                                         <  1:4607666/93‑1 (MQ=255)

                  gCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAAc                                                                        <  2:5820844/93‑1 (MQ=255)

                  gCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAAc                                                                        >  1:980360/1‑93 (MQ=255)

                   cACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACg                                                                       <  2:3294927/93‑1 (MQ=255)

                   cACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACg                                                                       <  2:270990/93‑1 (MQ=255)

                    aCGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACgg                                                                      <  2:1436138/93‑1 (MQ=255)

                    aCGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACgg                                                                      <  2:5122195/93‑1 (MQ=255)

                    aCGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACgg                                                                      <  1:2127281/93‑1 (MQ=255)

                      ggCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGcc                                                                    <  1:6669428/93‑1 (MQ=255)

                      ggCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGcc                                                                    >  2:4791717/1‑93 (MQ=255)

                       gCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACGGGTGCAGTTTAACGGCCt                                                                   <  1:1075988/93‑1 (MQ=255)

                       gCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCt                                                                   <  1:2787537/93‑1 (MQ=255)

                       gCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCt                                                                   <  1:2677648/93‑1 (MQ=255)

                       gCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCt                                                                   >  1:3103101/1‑93 (MQ=255)

                         tGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTgg                                                                 <  2:5538472/93‑1 (MQ=255)

                         tGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTgg                                                                 >  2:1222907/1‑93 (MQ=255)

                         tGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTgg                                                                 >  2:266850/1‑93 (MQ=255)

                          ggCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGgc                                                                <  1:4423494/93‑1 (MQ=255)

                          ggCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGgc                                                                >  2:3193187/1‑93 (MQ=255)

                           gCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGgcg                                                               >  2:4875849/1‑93 (MQ=255)

                              cagcagGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTg                                                            <  1:7530112/93‑1 (MQ=255)

                               agcagGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTgg                                                           <  1:5623497/93‑1 (MQ=255)

                                  agGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGa                                                        >  1:5344468/1‑93 (MQ=255)

                                  agGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGa                                                        <  1:1026016/93‑1 (MQ=255)

                                  agGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGa                                                        <  1:7223608/93‑1 (MQ=255)

                                   ggAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACGGGTGCAGTTTAACGGCCTGGCGCTGGTGAt                                                       <  1:4709005/93‑1 (MQ=255)

                                   ggAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGAt                                                       <  2:7216744/93‑1 (MQ=255)

                                     aaGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCCGGCGCTGGTGATTa                                                     <  2:6449077/93‑1 (MQ=255)

                                     aaGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTa                                                     >  2:1564346/1‑93 (MQ=255)

                                      aGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTAt                                                    <  1:847888/93‑1 (MQ=255)

                                       gCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATc                                                   <  2:5727988/93‑1 (MQ=255)

                                         cATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATcga                                                 <  2:2255688/93‑1 (MQ=255)

                                         cATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATcga                                                 <  1:5259035/93‑1 (MQ=255)

                                          aTCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATcgac                                                <  1:4042667/93‑1 (MQ=255)

                                          aTCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATcgac                                                <  1:293833/93‑1 (MQ=255)

                                           tCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATcgacg                                               >  2:5495290/1‑93 (MQ=255)

                                             gCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGaa                                             <  2:7304680/93‑1 (MQ=255)

                                             gCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGaa                                             >  2:74890/1‑93 (MQ=255)

                                             gCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGaa                                             >  1:509229/1‑93 (MQ=255)

                                              ccAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAAc                                            <  1:6474357/93‑1 (MQ=255)

                                              ccAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAAc                                            <  2:3538723/93‑1 (MQ=255)

                                               cAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATCATCGACGAACa                                           <  1:1769082/93‑1 (MQ=255)

                                               cAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACa                                           <  2:7420937/93‑1 (MQ=255)

                                                 gCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGc                                         >  2:186900/1‑93 (MQ=255)

                                                   ggTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCAt                                       <  2:1721101/93‑1 (MQ=255)

                                                    gTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATc                                      <  2:4015954/93‑1 (MQ=255)

                                                    gTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATc                                      <  1:30264/93‑1 (MQ=255)

                                                      cAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGt                                    <  2:3054009/93‑1 (MQ=255)

                                                      cAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGt                                    <  1:1695802/93‑1 (MQ=255)

                                                       aGGTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGtt                                   <  1:5970614/93‑1 (MQ=255)

                                                        ggTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGttt                                  >  2:5221898/1‑93 (MQ=255)

                                                        ggTGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGttt                                  >  1:319070/1‑93 (MQ=255)

                                                         gTGCAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGtttt                                 <  2:7060345/93‑1 (MQ=255)

                                                            cAGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGc                              <  1:2339224/93‑1 (MQ=255)

                                                             aGATGATTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCg                             <  1:2079099/93‑1 (MQ=255)

                                                              gatgatTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGt                            >  1:1371452/1‑93 (MQ=255)

                                                              gatgatTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGt                            <  1:6594509/93‑1 (MQ=255)

                                                                 gatTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCa                         <  1:7085736/93‑1 (MQ=255)

                                                                 gatTGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCa                         <  2:4927953/93‑1 (MQ=255)

                                                                   ttGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATc                       >  1:4380724/1‑93 (MQ=255)

                                                                   ttGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATc                       <  2:6719536/93‑1 (MQ=255)

                                                                    tGTCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCa                      <  2:3786321/93‑1 (MQ=255)

                                                                      tCGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGc                    >  2:3094307/1‑93 (MQ=255)

                                                                       cGGTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGCGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCg                   <  1:6207892/93‑1 (MQ=255)

                                                                         gTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTc                 <  1:5758856/93‑1 (MQ=255)

                                                                         gTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTc                 <  2:2440064/93‑1 (MQ=255)

                                                                         gTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTc                 <  1:6864859/93‑1 (MQ=255)

                                                                         gTACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATAGTTTTGGCGTGCATCAGCGTc                 >  1:2577979/1‑93 (MQ=255)

                                                                          tACACACGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTCt                <  2:4237349/93‑1 (MQ=255)

                                                                           acacacGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTCTg               <  1:3909999/93‑1 (MQ=255)

                                                                             acacGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTCTGGc             <  2:5975093/93‑1 (MQ=255)

                                                                              cacGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTCTGGCa            >  2:4536903/1‑93 (MQ=255)

                                                                               acGCCATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTCTGGCAt           <  1:2467021/93‑1 (MQ=255)

                                                                                  ccATCTTGCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGGGATTATCGCCGAACAGCATCGTTTTGCCGTGCCTCAGAGTATCACATtgt        >  1:1967612/1‑93 (MQ=255)

                                                                                  ccATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTCTGGCATtgt        >  1:5395658/1‑93 (MQ=255)

                                                                                   cATCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTCTGGCATtgtg       >  1:7327108/1‑93 (MQ=255)

                                                                                    aTCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTCTGGCATTGTgg      >  1:5500530/1‑93 (MQ=255)

                                                                                    aTCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTCTGGCATTGTgg      >  1:7442692/1‑93 (MQ=255)

                                                                                     tCTTCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTCTGGCATTGTggg     >  1:6326103/1‑93 (MQ=255)

                                                                                     tCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTCTGGCATTGTggg     >  2:299650/1‑93 (MQ=255)

                                                                                       ttCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTCTGGCATTGTGGgag   >  1:3360398/1‑93 (MQ=255)

                                                                                        tCCGGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTCTGGCATTGTGGgaga  >  2:6521389/1‑93 (MQ=255)

GGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTTTTGGCGTGCATCAGCGTCTGGCATTGTGGGAGA  >  NC_000913/3823893‑3824073

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