BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	657,814	G→A	A113T (GCC→ACC)	ybeM →	putative amidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	657,814	0	G	A	91.5%	757.8 / 34.3	260	A113T (GCC→ACC)	ybeM	putative amidase
Reads supporting (aligned to +/- strand): ref base G (9/13); new base A (120/118); total (129/131)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-04

Only 100 of 288 total aligned reads displayed.

GCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCC  >  NC_000913/657724‑657903

gcACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCCACGTTGCTTAcc                                                                                         >  2:4139256/1‑93 (MQ=255)

gcACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTAcc                                                                                         <  1:3872586/93‑1 (MQ=255)

 cacaGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCg                                                                                        >  1:5573059/1‑93 (MQ=255)

 cacaGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCg                                                                                        >  1:7917057/1‑93 (MQ=255)

   caGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCt                                                                                      >  1:3618122/1‑93 (MQ=255)

   caGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCt                                                                                      >  2:4733947/1‑93 (MQ=255)

   caGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCt                                                                                      >  2:6791310/1‑93 (MQ=255)

    aGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTc                                                                                     >  2:3822394/1‑93 (MQ=255)

     ggCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCg                                                                                    >  1:6851564/1‑93 (MQ=255)

     agCATGACAGGGAGCTGAAATCCTGGTCCTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCg                                                                                    <  1:8902109/92‑1 (MQ=255)

      gCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGt                                                                                   <  2:3119699/93‑1 (MQ=255)

      gCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGt                                                                                   >  2:7138532/1‑93 (MQ=255)

        aTTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTgc                                                                                 >  2:406451/1‑93 (MQ=255)

         ttACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGTCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTgcg                                                                                >  1:9090773/1‑93 (MQ=255)

         ttACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTgcg                                                                                >  1:7732771/1‑93 (MQ=255)

         ttACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTgcc                                                                                >  1:4942867/1‑92 (MQ=255)

         ttACAGGGAGCTGAAATCCGGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTgcg                                                                                <  1:1982165/93‑1 (MQ=255)

          tACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTgcgc                                                                               <  2:3023591/93‑1 (MQ=255)

          tACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTgcgc                                                                               >  1:1087100/1‑93 (MQ=255)

            cAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTg                                                                             <  1:8295783/93‑1 (MQ=255)

              gggAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGa                                                                           <  2:823193/93‑1 (MQ=255)

              gggAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGa                                                                           <  2:9159206/93‑1 (MQ=255)

              gggAGCTGAAATCCCGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGa                                                                           >  1:6008186/1‑93 (MQ=255)

                 aGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATac                                                                        <  2:7683321/93‑1 (MQ=255)

                  gCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATacc                                                                       >  2:3483930/1‑93 (MQ=255)

                    tGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATaccac                                                                     <  1:8680102/93‑1 (MQ=255)

                        aTCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGt                                                                 <  2:3287809/93‑1 (MQ=255)

                         tCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGtt                                                                <  1:1147405/93‑1 (MQ=255)

                            tGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAGGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTtata                                                             <  1:3331626/93‑1 (MQ=255)

                             ggTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTatat                                                            >  2:4315293/1‑93 (MQ=255)

                             ggTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTatat                                                            >  1:5158246/1‑93 (MQ=255)

                              gTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATAtg                                                           >  1:2435168/1‑93 (MQ=255)

                               tACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATAtgg                                                          <  1:7643408/93‑1 (MQ=255)

                                tcTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATAtggt                                                         <  1:2879936/92‑1 (MQ=255)

                                aCTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATAtggt                                                         <  1:1041076/93‑1 (MQ=255)

                                  ttCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATAtggtgg                                                       >  1:1643513/1‑93 (MQ=255)

                                   tCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTggc                                                      <  2:5903320/93‑1 (MQ=255)

                                   tCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATAtggtggg                                                      >  2:8439049/1‑92 (MQ=255)

                                   tCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTggc                                                      >  1:9210063/1‑93 (MQ=255)

                                   tCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTggc                                                      <  1:1492876/93‑1 (MQ=255)

                                    ccTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTggcg                                                     >  2:4680755/1‑93 (MQ=255)

                                      tGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTggcggc                                                   >  2:3687761/1‑93 (MQ=255)

                                      tGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTggcggc                                                   >  2:6925887/1‑93 (MQ=255)

                                      tGCCGCCTGGGTGCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTggcggc                                                   <  2:8523755/93‑1 (MQ=255)

                                          gccTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCgggg                                               <  1:3017284/93‑1 (MQ=255)

                                           ccTGGTTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCggggg                                              <  1:3682800/93‑1 (MQ=255)

                                           ccTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCggggg                                              >  1:4028266/1‑93 (MQ=255)

                                           ccTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCggggg                                              >  1:4104285/1‑93 (MQ=255)

                                           ccTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCggggg                                              >  1:9870789/1‑93 (MQ=255)

                                            cTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCgggggg                                             >  1:83607/1‑92 (MQ=255)

                                            cTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCgggggg                                             >  2:4962157/1‑92 (MQ=255)

                                             tGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAg                                            >  1:1587583/1‑93 (MQ=255)

                                             tGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAg                                            >  1:5391211/1‑93 (MQ=255)

                                             tGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAg                                            <  1:3645684/93‑1 (MQ=255)

                                             tGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAg                                            <  1:126421/93‑1 (MQ=255)

                                               ggTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTg                                          <  1:9299553/93‑1 (MQ=255)

                                               ggTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTg                                          <  2:3618122/93‑1 (MQ=255)

                                                gTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAgggg                                         >  2:1104724/1‑90 (MQ=255)

                                                      gggCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAAc                                   <  2:889271/93‑1 (MQ=255)

                                                      gggCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGGGGGCGGGAAc                                   >  2:518834/1‑93 (MQ=255)

                                                      gggCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAAc                                   <  2:5756707/93‑1 (MQ=255)

                                                       ggCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACa                                  <  2:1357939/93‑1 (MQ=255)

                                                       ggCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACa                                  >  1:4457071/1‑93 (MQ=255)

                                                        gCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGTGGGCGGGAACaa                                 >  1:693011/1‑93 (MQ=255)

                                                           gCTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACaaaaa                              <  1:9756417/93‑1 (MQ=255)

                                                            cTCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAat                             <  2:722134/93‑1 (MQ=255)

                                                             tCAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAata                            <  1:2043502/93‑1 (MQ=255)

                                                              cAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAatat                           >  2:214674/1‑93 (MQ=255)

                                                              cAAAGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAatat                           <  2:4397373/93‑1 (MQ=255)

                                                               aaaGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATc                          <  2:8364078/93‑1 (MQ=255)

                                                               aaaGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATc                          <  1:6522560/93‑1 (MQ=255)

                                                                aaGAGCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCg                         >  2:5254387/1‑93 (MQ=255)

                                                                 agagCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCgg                        <  1:6428799/93‑1 (MQ=255)

                                                                  gagCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGt                       <  2:5507/93‑1 (MQ=255)

                                                                    gCATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCa                     <  2:6614922/93‑1 (MQ=255)

                                                                     cATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCaa                    <  2:247103/93‑1 (MQ=255)

                                                                     cATCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCaa                    >  2:9304516/1‑93 (MQ=255)

                                                                      aTCACTGGTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCaaa                   <  2:2828229/93‑1 (MQ=255)

                                                                            ggTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGa             >  1:7115098/1‑93 (MQ=255)

                                                                            ggTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAAGCAAAATATCCGTCACAGCCGGa             >  1:4271865/1‑93 (MQ=255)

                                                                            ggTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGa             <  2:5711170/93‑1 (MQ=255)

                                                                            ggTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGa             <  1:5965642/93‑1 (MQ=255)

                                                                            ggTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGa             <  1:1789448/93‑1 (MQ=255)

                                                                             gTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGAt            >  1:14297/1‑93 (MQ=255)

                                                                             gTCAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGAt            >  1:4661330/1‑93 (MQ=255)

                                                                              tgaccGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGCCGGCGGGGGAGTGCGGAAAC‑AAAATATCGGTCAAAGCCGGATta          <  1:5010314/89‑1 (MQ=255)

                                                                               cAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATta          <  2:6234456/93‑1 (MQ=255)

                                                                               cAACGTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATta          <  2:171614/93‑1 (MQ=255)

                                                                                   gTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGa      >  1:7415849/1‑93 (MQ=255)

                                                                                   gTTGCTTACCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGAGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGa      <  2:504721/93‑1 (MQ=255)

                                                                                    ttgcttacCGCTCGTGCGCTGGATACCACCTGTTCTATGGTGGCGGCGAGGGGGAGCGGGGACAAAAATATCGGTCAAAGCCGGCTTATAGAt     >  2:5258414/1‑93 (MQ=255)

                                                                                    ttgcttacCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGAt     >  2:7562/1‑93 (MQ=255)

                                                                                    ttgcttacCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGAt     <  2:9239731/93‑1 (MQ=255)

                                                                                    ttgcttacCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGAt     >  1:940940/1‑93 (MQ=255)

                                                                                    ttgcttacCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGAt     >  1:8029945/1‑93 (MQ=255)

                                                                                     tgcttacCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATc    >  2:1799749/1‑93 (MQ=255)

                                                                                      gcttacCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATcc   >  2:9537585/1‑93 (MQ=255)

                                                                                       cttccCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGGGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATccc  >  2:9839270/1‑93 (MQ=255)

                                                                                       cttacCGCTCGTGCGCTGGATACCACTTGTTATATGGTGGCGGCGGGTGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATccc  <  1:2649315/93‑1 (MQ=255)

                                                                                       cttacCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATcct  >  2:7038698/1‑92 (MQ=255)

GCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCC  >  NC_000913/657724‑657903

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