BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,544,360	+G	coding (166/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,358	1	.	G	99.1%	307.8 / ‑3.8	112	coding (164/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.92e-01

Only 100 of 112 total aligned reads displayed.

GGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTC‑G‑‑C‑GGA‑CG‑TCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCA  >  NC_000913/3544300‑3544419

ggATGAATGCAGGGAGTGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGcc                                                 <  2:1612098/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                                               <  2:3473899/76‑1 (MQ=255)

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

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

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

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

    tAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGc                                             <  1:4063532/76‑1 (MQ=255)

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

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

      aTGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGccc                                           <  1:4367839/76‑1 (MQ=255)

      aTGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGccc                                           >  1:2868821/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:4676268/76‑1 (MQ=255)

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

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

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

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

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

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

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

            ggAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCa                                     >  1:3776752/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)

                tGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGGAGCCCTGCCCTGATCATggg                                 >  2:3815991/1‑76 (MQ=255)

                  ccGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGtt                               <  1:539961/76‑1 (MQ=255)

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

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

                    gCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGcttt                             >  1:1906616/1‑76 (MQ=255)

                      tGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTgg                           >  1:5076309/1‑76 (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:3459738/1‑76 (MQ=255)

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

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

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

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

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

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

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

                               ttATTGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCt                  >  2:1117334/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:2659045/76‑1 (MQ=255)

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

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

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

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

                                  ttGGCTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTggg               >  1:4789431/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:4370591/1‑76 (MQ=255)

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

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

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

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

                                      cTCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaaa           >  2:3403039/1‑76 (MQ=255)

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

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

                                       tCCATCAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCaaaa          >  2:2338542/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:5744834/76‑1 (MQ=255)

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

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

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

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

                                            cAAAGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTg     <  2:5423374/76‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                             aaaGCCGGTGTCGG‑‑CGGGA‑CGCTCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTgg    >  2:2991023/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  <  2:2179279/76‑1 (MQ=255)

GGATTAATGCAAGGAATGCCGCTGGATAAAGTTATTGGCTCCATCAAAGCCGGTGTC‑G‑‑C‑GGA‑CG‑TCGGTAGCCTTGCCCTGATCATGGGTTTTGGCGCAATGCTGGGCAAAATGCTGGCA  >  NC_000913/3544300‑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: `‑`