BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,544,357	+G	coding (163/1314 nt)	gntT →	gluconate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,544,356	1	.	G	99.1%	304.6 / ‑3.7	112	coding (162/1314 nt)	gntT	gluconate permease
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (57/54); total (58/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 121 total aligned reads displayed.

CCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTC‑G‑‑C‑GGA‑CG‑TCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCA  >  NC_000913/3544281‑3544419

ccTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTc                                                                       >  1:3824517/1‑76 (MQ=255)

ccTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTc                                                                       >  1:4985139/1‑76 (MQ=255)

 cTCGTGTCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTc‑g                                                                     <  1:2243582/76‑1 (MQ=255)

 cTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTc‑g                                                                     <  2:3460270/76‑1 (MQ=255)

 cTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTc‑g                                                                     <  1:409165/76‑1 (MQ=255)

 cTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTc‑g                                                                     <  1:5267509/76‑1 (MQ=255)

   cGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTc‑ggc                                                                   >  1:1000142/1‑74 (MQ=255)

    gTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTc‑ggcg                                                                  >  1:4239359/1‑73 (MQ=255)

                   ggATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGcc                                                 <  2:2497904/76‑1 (MQ=255)

                     aTTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCtt                                               <  1:139256/76‑1 (MQ=255)

                     aTTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCtt                                               <  2:4282857/76‑1 (MQ=255)

                     aTTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCtt                                               <  2:3473899/76‑1 (MQ=255)

                     aTTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCtt                                               >  1:4190895/1‑76 (MQ=255)

                     aTTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCtt                                               <  2:1866624/76‑1 (MQ=255)

                     aTTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCtt                                               <  1:4607170/76‑1 (MQ=255)

                     aTTAACGCAAGGAATGCCGCTGGATAAAGGTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCtt                                               >  1:4035020/1‑76 (MQ=255)

                      ttAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTg                                              <  2:2768587/76‑1 (MQ=255)

                        aaTGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGcc                                            <  1:3145897/76‑1 (MQ=255)

                         aTGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGccc                                           >  1:2868821/1‑76 (MQ=255)

                         aTGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGccc                                           >  1:2521441/1‑76 (MQ=255)

                         aTGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGGAGCCTTGtcc                                           >  1:3223225/1‑76 (MQ=255)

                          tGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCt                                          >  1:2295828/1‑76 (MQ=255)

                           gCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTg                                         >  1:4240413/1‑76 (MQ=255)

                           gCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTg                                         >  1:3667360/1‑76 (MQ=255)

                           gCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTg                                         <  2:1586554/76‑1 (MQ=255)

                            cAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGa                                        <  1:4625560/76‑1 (MQ=255)

                             aaGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGAt                                       <  2:2116731/76‑1 (MQ=255)

                             aaGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGAt                                       >  1:5587924/1‑76 (MQ=255)

                             aaGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGAt                                       <  1:251881/76‑1 (MQ=255)

                             aaGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGAt                                       <  1:1938261/76‑1 (MQ=255)

                              aGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATc                                      <  1:4948689/76‑1 (MQ=255)

                               ggAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGACCC‑CCGGTAGCCCTGCCCtcatca                                     >  2:5480547/1‑76 (MQ=255)

                               ggAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCCGATCa                                     >  2:4452916/1‑76 (MQ=255)

                                gAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCAt                                    <  2:4440872/76‑1 (MQ=255)

                                 aaTGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATg                                   <  2:2631144/76‑1 (MQ=255)

                                   tGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGGAGCCTTGCCCTGATCATggg                                 >  1:5423375/1‑76 (MQ=255)

                                     ccGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCCTGCCCTGATCATGGGtt                               >  2:4780317/1‑76 (MQ=255)

                                       tgtGGATAATGATATTGGTGCCGTCAGAGCCGGGGTCGG‑‑CGGGA‑GGCTCGGTATCGTTGCCCTGATCATGGGtttt                             <  2:212245/74‑1 (MQ=255)

                                       gCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGtttt                             <  2:630269/76‑1 (MQ=255)

                                        cTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTg                            <  1:519394/76‑1 (MQ=255)

                                          ggATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGgc                          <  2:4771712/76‑1 (MQ=255)

                                          ggATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGgc                          <  2:4814330/76‑1 (MQ=255)

                                           gATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGgcg                         <  1:5618031/76‑1 (MQ=255)

                                            gtAAAGTGGTGGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGG‑CGTTCGGTAGCCTTGCCCTGATCATGGGTTTTTgcgg                        <  2:579426/75‑2 (MQ=255)

                                            aTAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGgcgc                        >  1:5072361/1‑76 (MQ=255)

                                            aTAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGgcgc                        >  2:673719/1‑76 (MQ=255)

                                            aTAAAGTTATTGGCTCCATCAAAGCCGGTGTC‑GGAG‑GGA‑CGCTCGGTAGCCTTGCCCTGATCCTGGGGTTTGgcgc                        >  2:1329153/1‑76 (MQ=255)

                                             tAAAGTTATTGGCTCCCTCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCa                       >  2:4004585/1‑76 (MQ=255)

                                             tAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGCTCATGGGTTTTGGCGCa                       >  2:3958780/1‑76 (MQ=255)

                                             tAAAGTTATTGGCGCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCa                       <  2:3918918/76‑1 (MQ=255)

                                             tAAAGTTATGGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGTGTGTTGgcgaa                       <  2:1532820/76‑3 (MQ=255)

                                              aaaGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCaa                      >  1:1924857/1‑76 (MQ=255)

                                              aaaGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCaa                      >  1:3459738/1‑76 (MQ=255)

                                               aaGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAAt                     <  1:2985884/76‑1 (MQ=255)

                                                aGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATg                    <  1:933246/76‑1 (MQ=255)

                                                 gtgaTTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGc                   <  2:328401/73‑1 (MQ=255)

                                                 gTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGc                   >  2:4163251/1‑76 (MQ=255)

                                                 gTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGc                   <  1:2159679/76‑1 (MQ=255)

                                                 gTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGc                   <  1:3559093/76‑1 (MQ=255)

                                                  ttATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCt                  >  2:1753366/1‑76 (MQ=255)

                                                  ttATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCt                  >  1:909486/1‑76 (MQ=255)

                                                  ttATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGGTTTGGCGCAATGCt                  >  1:1638130/1‑76 (MQ=255)

                                                  ttATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCCGATCATGGGGTTTTGCGCAATGCt                  >  2:321815/1‑76 (MQ=255)

                                                   tATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTg                 <  1:3884097/76‑1 (MQ=255)

                                                   tATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTg                 <  2:5115693/76‑1 (MQ=255)

                                                   tATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTg                 >  1:5146332/1‑76 (MQ=255)

                                                   tATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTg                 <  1:2659045/76‑1 (MQ=255)

                                                    aTTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTgg                <  1:1889284/76‑1 (MQ=255)

                                                     ttGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTggg               <  2:3572832/76‑1 (MQ=255)

                                                     ttGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTggg               >  1:2413846/1‑76 (MQ=255)

                                                      tGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGc              >  2:2887020/1‑76 (MQ=255)

                                                       ggCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCa             >  1:613195/1‑76 (MQ=255)

                                                        gCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaa            >  1:534100/1‑76 (MQ=255)

                                                        gCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaa            >  2:3447709/1‑76 (MQ=255)

                                                        gCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaa            >  1:4370591/1‑76 (MQ=255)

                                                        gCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCCGATCATGGGTTTTGGCGCAATGCTGGGCaa            >  2:3183513/1‑76 (MQ=255)

                                                         cTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaaa           <  1:1468813/76‑1 (MQ=255)

                                                          tCCATCAGAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaaaa          <  2:1718865/76‑1 (MQ=255)

                                                          tCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaaaa          >  2:2338542/1‑76 (MQ=255)

                                                          tCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaaaa          >  1:4907051/1‑76 (MQ=255)

                                                           ccATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTTCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAAt         <  2:1022652/76‑1 (MQ=255)

                                                           ccATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAAt         <  2:5076308/76‑1 (MQ=255)

                                                            cATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATg        <  2:3864426/76‑1 (MQ=255)

                                                            cATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATg        >  2:4605687/1‑76 (MQ=255)

                                                             aTCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGc       >  1:3901027/1‑76 (MQ=255)

                                                             aTCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGc       >  1:425972/1‑76 (MQ=255)

                                                             aTCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGc       <  1:5744834/76‑1 (MQ=255)

                                                              tCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCCTGGGTTTTGGCGCAATGCTGGGCAAAATGCt      >  2:3133455/1‑76 (MQ=255)

                                                              tCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCt      <  1:4656250/76‑1 (MQ=255)

                                                               cAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTg     >  1:524487/1‑76 (MQ=255)

                                                                aaaGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTgg    >  1:2106846/1‑76 (MQ=255)

                                                                aaaGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTgg    >  2:2991023/1‑76 (MQ=255)

                                                                aaaGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTgg    >  1:2980210/1‑76 (MQ=255)

                                                                aaaGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTgg    <  1:3998541/76‑1 (MQ=255)

                                                                aaaGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTgg    >  2:1290999/1‑76 (MQ=255)

                                                                aaaGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTgg    >  1:5410459/1‑76 (MQ=255)

                                                                 aaGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGc   >  2:4444845/1‑76 (MQ=255)

                                                                  aGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCCATGCTGGGCAAAATGCTGGCa  >  2:1014804/1‑76 (MQ=255)

                                                                  aGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCa  >  1:4375502/1‑76 (MQ=255)

                                                                  aGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCa  <  1:2173756/76‑1 (MQ=255)

CCTCGTGGCGCTTGCTGTTGGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTC‑G‑‑C‑GGA‑CG‑TCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCA  >  NC_000913/3544281‑3544419

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 8 ≤` `ATCG/ATCG` `< 18 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`