BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,044,033	(C)_5→6	coding (1149/7152 nt)	b1978 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,044,028	1	.	C	91.8%	311.3 / 21.7	110	coding (1144/7152 nt)	b1978	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/8); new base C (64/37); total (65/45)
Fisher's exact test for biased strand distribution p-value = 3.21e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

Only 100 of 127 total aligned reads displayed.

TCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑CCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTACCT  >  NC_000913/2043936‑2044115

tCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACtat‑a                                                                                          <  2:2442755/93‑1 (MQ=255)

tCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACtat‑a                                                                                          <  1:1923187/93‑1 (MQ=255)

tCTATGAACAGTATTATGGCGATGAAGTCGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACtat‑a                                                                                          <  1:4240752/93‑1 (MQ=255)

 cTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑c                                                                                        <  2:2650368/93‑1 (MQ=255)

  tATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑cc                                                                                       <  2:3408421/93‑1 (MQ=255)

    tGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑cccc                                                                                     <  1:1054648/93‑1 (MQ=255)

      aaCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑cccccc                                                                                   >  2:4339351/1‑92 (MQ=255)

      aaCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑CCCCCt                                                                                   <  2:2966114/93‑1 (MQ=255)

       aCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGCCGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑cccccct                                                                                  <  2:1525561/93‑3 (MQ=255)

       aCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑cccccct                                                                                  >  1:1554251/1‑91 (MQ=255)

       aCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑cccccct                                                                                  <  1:3422196/93‑3 (MQ=255)

       aCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑cccccct                                                                                  <  2:1099724/93‑3 (MQ=255)

       aCAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑cccccct                                                                                  <  2:1651614/93‑3 (MQ=255)

        cAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑cccccctt                                                                                 >  2:2439615/1‑90 (MQ=255)

         aGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑ccccccttc                                                                                <  1:1530526/93‑5 (MQ=255)

           tattatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTccc                                                                               >  1:4427159/1‑93 (MQ=255)

           tattatGGCGATGAAGTGGCCCGGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGAACTTAACTAT‑ACCCCCCTTccc                                                                               <  1:1749349/93‑1 (MQ=255)

           tatgatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTccc                                                                               <  1:3214557/93‑1 (MQ=255)

            attatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTcccc                                                                              >  2:2146165/1‑92 (MQ=255)

             ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGc                                                                             >  1:2979259/1‑93 (MQ=255)

             ttatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGc                                                                             >  1:4191394/1‑93 (MQ=255)

              tatGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCt                                                                            <  1:3654297/93‑1 (MQ=255)

                tGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCtga                                                                          >  2:1668696/1‑93 (MQ=255)

                     aTGAAGTGCCCGTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGAcc                                                                     <  2:923122/93‑1 (MQ=255)

                      tGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCt                                                                    >  1:303837/1‑93 (MQ=255)

                      tGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCt                                                                    >  1:698911/1‑93 (MQ=255)

                      tGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCt                                                                    <  1:4418275/93‑1 (MQ=255)

                      tGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCt                                                                    >  2:3071643/1‑93 (MQ=255)

                       gAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCtt                                                                   <  2:2568035/93‑1 (MQ=255)

                        aaGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTc                                                                  >  1:681748/1‑93 (MQ=255)

                        aaGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTc                                                                  >  2:1478776/1‑93 (MQ=255)

                        aaGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTc                                                                  <  2:3047369/93‑1 (MQ=255)

                         aGTGGCCCTGTTCGATAAAGACGATCTGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATGC‑CCCCCTTCCCGCTGATGACCTTCa                                                                 <  2:3060914/93‑1 (MQ=255)

                         aGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCa                                                                 >  2:1886272/1‑93 (MQ=255)

                          gTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAg                                                                <  2:293374/93‑1 (MQ=255)

                           gggCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAgc                                                               <  1:2770645/92‑1 (MQ=255)

                               ccTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCgg                                                           >  1:440536/1‑93 (MQ=255)

                               ccTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCgg                                                           >  2:4470694/1‑93 (MQ=255)

                               ccTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCgg                                                           <  1:2422937/93‑1 (MQ=255)

                                cTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGa                                                          >  1:3662500/1‑93 (MQ=255)

                                cTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGa                                                          >  2:3585506/1‑93 (MQ=255)

                                cTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGa                                                          >  1:608238/1‑93 (MQ=255)

                                cTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGa                                                          >  1:891892/1‑93 (MQ=255)

                                cTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGa                                                          <  2:1840621/93‑1 (MQ=255)

                                  gTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGc                                                        >  1:4268505/1‑93 (MQ=255)

                                     tgATAAAGACAATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAAc                                                     <  2:2472548/92‑1 (MQ=255)

                                     cGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTATCTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAAc                                                     <  1:3304289/93‑1 (MQ=255)

                                      gATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACg                                                    >  1:2058447/1‑93 (MQ=255)

                                        tAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGcc                                                  >  2:663007/1‑93 (MQ=255)

                                        tAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGcc                                                  >  2:1649579/1‑93 (MQ=255)

                                          aaGACGATCGGCAAAGTAATCCTCATGCCATAACCTCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAg                                                <  1:4019901/93‑1 (MQ=255)

                                          aaGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAg                                                <  2:566785/93‑1 (MQ=255)

                                          aaGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAAAGCCAg                                                >  2:613643/1‑93 (MQ=255)

                                           aGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAgg                                               >  1:1786356/1‑93 (MQ=255)

                                               gATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTaa                                           >  2:2680751/1‑93 (MQ=255)

                                               caccGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTaa                                           <  2:3087803/90‑1 (MQ=255)

                                                  cGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACa                                        <  2:1010353/93‑1 (MQ=255)

                                                   gtcAAAGTAATCCTCAGGCCATAACCGCTGGACTTAACTATCC‑CCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAg                                       <  2:2787973/91‑1 (MQ=255)

                                                   cgCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAg                                       <  2:3133567/92‑1 (MQ=255)

                                                    gCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATCC‑CCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAgg                                      <  1:953741/93‑1 (MQ=255)

                                                    gCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAgg                                      >  2:2119008/1‑93 (MQ=255)

                                                    gCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAgg                                      >  1:1884368/1‑93 (MQ=255)

                                                    gCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAgg                                      <  1:4668131/93‑1 (MQ=255)

                                                     cAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTATCC‑CCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAggg                                     <  1:31416/93‑1 (MQ=255)

                                                     cAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAggg                                     >  1:1977251/1‑93 (MQ=255)

                                                     cAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAggg                                     <  1:1430526/93‑1 (MQ=255)

                                                      aaaGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCTGGGTAAACAGGGc                                    <  2:4685670/93‑1 (MQ=255)

                                                      aaaGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGc                                    >  1:4599510/1‑93 (MQ=255)

                                                         gTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGaa                                 >  1:3234634/1‑93 (MQ=255)

                                                          tAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGaaa                                >  1:1443010/1‑93 (MQ=255)

                                                          tAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGaaa                                >  1:4634709/1‑93 (MQ=255)

                                                          tAATCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGaaa                                >  2:3070747/1‑93 (MQ=255)

                                                           aaTCCTCATGCCATAACCGCTGGACTTAACTATCC‑CCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGaaaa                               <  2:4306773/93‑1 (MQ=255)

                                                            aTCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAAt                              <  2:3169996/93‑1 (MQ=255)

                                                             tCCTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATg                             <  1:2839584/93‑1 (MQ=255)

                                                              ccTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGa                            <  1:314027/93‑1 (MQ=255)

                                                               cTCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGac                           >  2:1212234/1‑93 (MQ=255)

                                                                tCATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGaca                          >  2:2359117/1‑93 (MQ=255)

                                                                 cATGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGacac                         <  2:718206/93‑1 (MQ=255)

                                                                   tGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACAccc                       <  2:3887796/93‑1 (MQ=255)

                                                                   tGCCATAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACAccc                       >  2:4177931/1‑93 (MQ=255)

                                                                       aTAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGTCAGGGTAAACAGGGCGAAAATGACACCCGttt                   >  2:2023375/1‑93 (MQ=255)

                                                                       aTAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGttt                   >  2:237644/1‑93 (MQ=255)

                                                                        tAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAACCAGGGCGAAAATGACACCCGtctt                  >  2:3891188/1‑93 (MQ=255)

                                                                        tAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGtttt                  >  1:3396856/1‑93 (MQ=255)

                                                                        tAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGtttt                  <  1:2857463/93‑1 (MQ=255)

                                                                        tAACCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGtttt                  >  2:3509061/1‑93 (MQ=255)

                                                                          aCCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGc                <  2:2300469/93‑1 (MQ=255)

                                                                          aCCGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGc                >  2:922071/1‑93 (MQ=255)

                                                                            cGCTGGACTTAACTATCC‑CCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCg              <  1:2586721/93‑1 (MQ=255)

                                                                            cGCTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCg              <  1:2940906/93‑1 (MQ=255)

                                                                              cTGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTc            <  1:651842/93‑1 (MQ=255)

                                                                               tGGACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCg           >  1:329270/1‑93 (MQ=255)

                                                                                 gACTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGAt         >  2:328763/1‑93 (MQ=255)

                                                                                  aCTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGAtt        >  2:2610033/1‑93 (MQ=255)

                                                                                   cTTAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGAttt       <  1:1393736/93‑1 (MQ=255)

                                                                                    ttAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGAtttt      >  2:1417310/1‑93 (MQ=255)

                                                                                     tAACTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTa     >  2:2452996/1‑93 (MQ=255)

                                                                                        cTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTacct  <  2:4318915/93‑1 (MQ=255)

                                                                                        cTAT‑ACCCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTacct  <  2:276539/93‑1 (MQ=255)

TCTATGAACAGTATTATGGCGATGAAGTGGCCCTGTTCGATAAAGACGATCGGCAAAGTAATCCTCATGCCATAACCGCTGGACTTAACTAT‑A‑CCCCCTTCCCGCTGATGACCTTCAGCGCGGAGCAACGCCAGGGTAAACAGGGCGAAAATGACACCCGTTTTGCCGTCGATTTTACCT  >  NC_000913/2043936‑2044115

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