BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,221,058	A→G	Q285R (CAG→CGG)	dld →	D‑lactate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,221,058	0	A	G	86.8%	457.8 / 62.8	190	Q285R (CAG→CGG)	dld	D‑lactate dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (11/14); new base G (84/81); total (95/95)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.62e-01

Only 100 of 205 total aligned reads displayed.

GTTGCGCCGGGAAGCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCAGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAATATATGCACCGGGATATCT  >  NC_000913/2220968‑2221147

  tGCGCCGGGAAGCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCg                                                                                       <  2:3166094/93‑1 (MQ=255)

   gcgcCGGGAAGCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAAccggccgg                                                                                      >  1:4112543/1‑93 (MQ=255)

    cgcCGGGAAGCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGa                                                                                     <  2:1807577/93‑1 (MQ=255)

    cgcCGGGAAGCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGa                                                                                     <  1:6485001/93‑1 (MQ=255)

    cgcCGGGAAGCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGa                                                                                     >  1:6045735/1‑93 (MQ=255)

      ccGGGAAGCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCAGCCGGAAg                                                                                   <  1:3973876/93‑1 (MQ=255)

       cGGGAAGCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGt                                                                                  <  2:587293/93‑1 (MQ=255)

       cGGGAAGCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGt                                                                                  <  1:6601292/93‑1 (MQ=255)

        gggAAGCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTg                                                                                 >  2:691589/1‑93 (MQ=255)

         ggAAGCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGc                                                                                >  2:7433033/1‑93 (MQ=255)

         ggAAGCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGc                                                                                <  1:5085360/93‑1 (MQ=255)

         ggAAGCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGc                                                                                >  1:4342597/1‑93 (MQ=255)

           aaGCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTg                                                                              <  1:3624235/93‑1 (MQ=255)

           aaGCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTg                                                                              >  1:1820945/1‑93 (MQ=255)

             gCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGAc                                                                            <  2:4009204/93‑1 (MQ=255)

             gCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGAc                                                                            <  1:4427464/93‑1 (MQ=255)

              cTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGAcc                                                                           <  1:4267302/93‑1 (MQ=255)

               tGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCg                                                                          >  1:91897/1‑93 (MQ=255)

               tGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCg                                                                          >  1:1511564/1‑93 (MQ=255)

                ggCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGa                                                                         >  1:2701790/1‑93 (MQ=255)

                  cGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGaaa                                                                       <  1:7205006/93‑1 (MQ=255)

                  cGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGaaa                                                                       <  1:2994423/93‑1 (MQ=255)

                  cGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGaaa                                                                       <  2:3836971/93‑1 (MQ=255)

                   ggTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAAt                                                                      <  1:742782/93‑1 (MQ=255)

                    gTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATc                                                                     >  1:4852183/1‑93 (MQ=255)

                    gTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATc                                                                     >  1:3442126/1‑93 (MQ=255)

                    gTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCAGCCGGAAGTGCTGACCGAAATc                                                                     >  2:1187865/1‑93 (MQ=255)

                     tCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATcc                                                                    >  1:2690136/1‑93 (MQ=255)

                       tttGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATccgc                                                                  <  1:1549404/93‑1 (MQ=255)

                       tttGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATccgc                                                                  <  2:6410662/93‑1 (MQ=255)

                            aGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCa                                                             <  2:2876853/93‑1 (MQ=255)

                             gTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCat                                                            >  1:7117288/1‑93 (MQ=255)

                                cGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATAtt                                                         <  2:5297061/93‑1 (MQ=255)

                                 gTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTc                                                        <  2:2690525/93‑1 (MQ=255)

                                  tCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCt                                                       <  1:6896211/93‑1 (MQ=255)

                                  tCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCt                                                       <  1:2840320/93‑1 (MQ=255)

                                   cTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTg                                                      <  1:1735052/93‑1 (MQ=255)

                                    ttGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTgg                                                     >  2:363955/1‑93 (MQ=255)

                                    ttGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTgg                                                     >  1:3145833/1‑93 (MQ=255)

                                     tGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGc                                                    >  1:3815377/1‑93 (MQ=255)

                                      gATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCt                                                   >  1:2640698/1‑93 (MQ=255)

                                       aTACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTa                                                  >  2:4687178/1‑93 (MQ=255)

                                        tACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTaa                                                 >  2:4896932/1‑93 (MQ=255)

                                           cTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACtt                                              <  1:7363571/93‑1 (MQ=255)

                                            ttCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTc                                             >  1:63858/1‑93 (MQ=255)

                                            ttCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTc                                             <  1:5166349/93‑1 (MQ=255)

                                              cGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCAGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGa                                           <  2:3788553/93‑1 (MQ=255)

                                               gAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGaa                                          >  2:3377/1‑93 (MQ=255)

                                                aaGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGaaa                                         <  1:536359/93‑1 (MQ=255)

                                                 aGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGaaaa                                        <  1:6631211/93‑1 (MQ=255)

                                                 aGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGaaaa                                        <  1:7288508/93‑1 (MQ=255)

                                                   cGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATc                                      <  1:5567309/93‑1 (MQ=255)

                                                    ggAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCt                                     <  1:3027536/93‑1 (MQ=255)

                                                     gAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTg                                    >  2:2754190/1‑93 (MQ=255)

                                                      aaaaaaaTCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGc                                   <  2:2569761/93‑1 (MQ=255)

                                                      aaaaaaaTCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGc                                   >  2:3440188/1‑93 (MQ=255)

                                                        aaaaaTCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCg                                 <  2:6968676/93‑1 (MQ=255)

                                                          aaaTCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGt                               <  2:1941480/93‑1 (MQ=255)

                                                           aaTCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGtt                              <  1:6758325/93‑1 (MQ=255)

                                                            aTCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTg                             <  1:4766898/93‑1 (MQ=255)

                                                            aTCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTg                             >  1:1275143/1‑93 (MQ=255)

                                                            aTCAGCAGGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTg                             <  2:3550286/93‑1 (MQ=255)

                                                               agcagGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCg                          <  1:116786/93‑1 (MQ=255)

                                                                gcagGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCgg                         >  1:167222/1‑93 (MQ=255)

                                                                gcagGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCgg                         >  2:1111132/1‑93 (MQ=255)

                                                                 cagGTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCggg                        <  2:7004253/93‑1 (MQ=255)

                                                                   ggTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGa                      >  2:483883/1‑93 (MQ=255)

                                                                   ggTGTTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGa                      >  1:2375697/1‑93 (MQ=255)

                                                                    gtgtTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGaa                     <  1:5079029/93‑1 (MQ=255)

                                                                    gtgtTTTATATCGGCACCAACCGGCCGGAAGGGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGaa                     >  1:3673049/1‑93 (MQ=255)

                                                                      gtTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAata                   <  2:5621107/93‑1 (MQ=255)

                                                                      gtTTTATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAata                   >  2:1592282/1‑93 (MQ=255)

                                                                      gtTTTATATCGGCACCAACCAGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAata                   <  2:5898009/93‑1 (MQ=255)

                                                                       ttttATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAatat                  >  2:2408091/1‑93 (MQ=255)

                                                                       ttttATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAatat                  >  1:1427831/1‑93 (MQ=255)

                                                                       ttttATATCGGCACCAACCAGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAatat                  <  2:5208833/93‑1 (MQ=255)

                                                                        tttATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAatata                 <  2:6306153/93‑1 (MQ=255)

                                                                        tttATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAatata                 <  2:6065739/93‑1 (MQ=255)

                                                                        tttATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAatata                 <  2:6682480/93‑1 (MQ=255)

                                                                        tttATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAatata                 >  2:4624498/1‑93 (MQ=255)

                                                                        tttATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAatata                 <  1:1840919/93‑1 (MQ=255)

                                                                         ttATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAatatat                >  1:1034773/1‑93 (MQ=255)

                                                                         ttATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAatatat                >  2:5077828/1‑93 (MQ=255)

                                                                         ttATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAatatat                <  1:2133760/93‑1 (MQ=255)

                                                                         ttATATCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAatatat                <  2:7136002/93‑1 (MQ=255)

                                                                          tataTCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAATATATg               >  1:5796179/1‑93 (MQ=255)

                                                                           atatCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAATATATGc              >  2:5793878/1‑93 (MQ=255)

                                                                           atatCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAATATATGc              <  2:283457/93‑1 (MQ=255)

                                                                            tatCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAATATATGCa             >  2:1698071/1‑93 (MQ=255)

                                                                            tatCGGCACCAACCAGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAATATATGCa             >  2:3333802/1‑93 (MQ=255)

                                                                             atCGGCACCAACCGGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAATATATGCAc            <  2:4161248/93‑1 (MQ=255)

                                                                              tCGGCACCAACCAGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAATATATGCAcc           >  2:2303276/1‑93 (MQ=255)

                                                                                ggCACCAACCAGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAATATATGCACCgg         <  2:1392933/93‑1 (MQ=255)

                                                                                    ccaaccgGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAATATATGCACCGGGata     >  1:2780238/1‑93 (MQ=255)

                                                                                    ccaaccgGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAATATATGCACCGGGata     >  2:4822651/1‑93 (MQ=255)

                                                                                     caaccgGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAATATATGCACCGGGatat    >  1:4384122/1‑93 (MQ=255)

                                                                                     caaccgGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAATATATGCACCGGGatat    <  1:5882617/93‑1 (MQ=255)

                                                                                      aaccgGCCGGAAGTGCTGACCGAAATCCGCCGTCGTATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAATATATGCACCGGGATATc   <  2:4464344/93‑1 (MQ=255)

                                                                                      aaccgGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAATATATGCACCGGGATATc   <  1:7127748/93‑1 (MQ=255)

                                                                                       accgGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAATATATGCACCGGGATATCt  <  2:2463666/93‑1 (MQ=255)

GTTGCGCCGGGAAGCTGGCGGTCTTTGCAGTACGTCTTGATACCTTCGAAGCGGAAAAAAATCAGCAGGTGTTTTATATCGGCACCAACCAGCCGGAAGTGCTGACCGAAATCCGCCGTCATATTCTGGCTAACTTCGAAAATCTGCCGGTTGCCGGGGAATATATGCACCGGGATATCT  >  NC_000913/2220968‑2221147

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