BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,633,590	A→G	G94G (GGT→GGC)	b1547 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,633,590	0	A	G	99.1%	340.5 / ‑4.9	115	G94G (GGT→GGC)	b1547	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (64/50); total (64/51)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 116 total aligned reads displayed.

AAGTTTCGGGTATGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCACCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCGTCGGTCGTATCCCGTCGT  >  NC_000913/1633500‑1633677

aaGTTTCGGGTATGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGcc                                                                                       <  2:3869783/93‑1 (MQ=18)

 aGTTTCGGGTATGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCg                                                                                      >  1:627211/1‑93 (MQ=21)

  gTTTCGGGTATGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGa                                                                                     <  2:713211/93‑1 (MQ=18)

  gTTTCGGGTATGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGa                                                                                     >  2:387885/1‑93 (MQ=18)

   tttCGGGTATGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGAc                                                                                    >  2:1580338/1‑93 (MQ=18)

   tttCGGGTATGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGAc                                                                                    <  2:1352534/93‑1 (MQ=17)

    ttCGGGTATGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACt                                                                                   <  2:103537/93‑1 (MQ=17)

      cGGGTATGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTgg                                                                                 <  1:1801575/93‑1 (MQ=17)

      cGGGTATGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTgg                                                                                 >  1:2736154/1‑93 (MQ=21)

       ggggATGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGCGCGCCGACTGGa                                                                                >  2:4515686/1‑93 (MQ=255)

        ggTATGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGAt                                                                               >  1:4238384/1‑93 (MQ=25)

        ggTATGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGAt                                                                               <  2:3184540/93‑1 (MQ=21)

        ggTATGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGAt                                                                               >  1:2835203/1‑93 (MQ=25)

         gTATGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATa                                                                              <  2:1713888/93‑1 (MQ=25)

            tGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATaaaa                                                                           <  2:948815/93‑1 (MQ=255)

            tGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATaaaa                                                                           <  2:2193019/93‑1 (MQ=255)

            tGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATaaaa                                                                           >  2:3932474/1‑93 (MQ=255)

            tGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATaaaa                                                                           >  1:2924732/1‑93 (MQ=255)

             gCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAAc                                                                          <  1:951969/93‑1 (MQ=255)

             gCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAAc                                                                          <  1:805930/93‑1 (MQ=255)

              cAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAAct                                                                         <  1:1884182/93‑1 (MQ=255)

               aGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAActc                                                                        <  1:4300860/93‑1 (MQ=255)

               aGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAActc                                                                        <  1:2469488/93‑1 (MQ=255)

               aGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAActc                                                                        <  2:2497629/93‑1 (MQ=255)

                gATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAActct                                                                       >  2:593314/1‑93 (MQ=255)

                gATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAActct                                                                       >  1:3743245/1‑93 (MQ=255)

                gATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAActct                                                                       >  2:3578603/1‑93 (MQ=255)

                  tttGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGc                                                                     >  1:2792354/1‑93 (MQ=255)

                   ttGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCa                                                                    >  1:333668/1‑93 (MQ=255)

                      tCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCagga                                                                 <  1:2019071/93‑1 (MQ=255)

                      tCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCagga                                                                 >  1:496638/1‑93 (MQ=255)

                      tCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCagga                                                                 <  1:45912/93‑1 (MQ=255)

                      tCAAAAGTCTGCACCTGCATCAGGGCATAACAAGACGGCACGGTATCTGATGGCGACGGGATTGGTGCGCAGACTGGATACAACGCTGCagga                                                                 <  2:4579007/93‑1 (MQ=14)

                       cAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCaggag                                                                <  1:3996060/93‑1 (MQ=255)

                       cAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCaggag                                                                <  1:1810327/93‑1 (MQ=255)

                       cAAAAGTCTGCCCCGGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCaggag                                                                <  2:548151/93‑1 (MQ=255)

                        aaaaGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCaggagg                                                               >  2:631267/1‑93 (MQ=255)

                          aaGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGAt                                                             >  1:2263053/1‑93 (MQ=255)

                          aaGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGAt                                                             >  2:1298240/1‑93 (MQ=255)

                            gTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGa                                                           >  1:3493518/1‑93 (MQ=255)

                            gTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGa                                                           >  2:1785721/1‑93 (MQ=255)

                             tCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAg                                                          >  2:923485/1‑93 (MQ=255)

                              cTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGc                                                         <  2:2166591/93‑1 (MQ=255)

                               tGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGcc                                                        >  2:4506929/1‑93 (MQ=255)

                                 ccccTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGa                                                      >  2:3332433/1‑93 (MQ=255)

                                 acacTGCATAAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGa                                                      <  1:4228447/90‑1 (MQ=18)

                                  cccTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAg                                                     <  2:2321833/93‑1 (MQ=255)

                                  cccTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAg                                                     >  1:3024331/1‑93 (MQ=255)

                                  cccTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAg                                                     <  2:803199/93‑1 (MQ=255)

                                   cgtGCATCAGGTCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACGGGATAAAACTCTGCAGGAGGATGAGCCGAgg                                                    <  2:506199/91‑1 (MQ=18)

                                   ccTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAgg                                                    >  2:3330188/1‑93 (MQ=255)

                                   ccTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAgg                                                    >  1:3974190/1‑93 (MQ=255)

                                   ccTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAgg                                                    >  1:3700951/1‑93 (MQ=255)

                                    cTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAggg                                                   >  1:1228614/1‑92 (MQ=255)

                                    cTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGgt                                                   >  2:1621594/1‑93 (MQ=255)

                                    cTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGgt                                                   >  2:2076649/1‑93 (MQ=255)

                                      gCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGgtgt                                                 >  1:807151/1‑93 (MQ=255)

                                        aTCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTaa                                               >  1:2914933/1‑93 (MQ=255)

                                        aTCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTaa                                               >  2:4371697/1‑93 (MQ=255)

                                        aTCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTaa                                               <  1:4347914/93‑1 (MQ=255)

                                        aTCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGGGTaa                                               >  2:2029790/1‑93 (MQ=255)

                                        aTCAGGGCATAACCAGACGGAACGGGATCTGATGGCCGCGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTaa                                               <  2:378558/93‑1 (MQ=255)

                                         tCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTaaa                                              <  1:2889245/93‑1 (MQ=255)

                                         tCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTaaa                                              >  1:588066/1‑93 (MQ=255)

                                           aGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGc                                            >  1:4254184/1‑93 (MQ=255)

                                              gCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTgc                                         >  1:3714991/1‑93 (MQ=255)

                                                aTAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTgcgc                                       >  2:3869056/1‑93 (MQ=255)

                                                 tAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGcc                                      >  2:4010255/1‑93 (MQ=255)

                                                 tAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGcc                                      <  1:2735688/93‑1 (MQ=255)

                                                  aaCCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGccc                                     <  1:4067584/93‑1 (MQ=255)

                                                    ccAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCAc                                   <  1:980106/93‑1 (MQ=255)

                                                       gACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGc                                <  2:3145736/93‑1 (MQ=255)

                                                       gACGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGc                                >  2:162339/1‑93 (MQ=255)

                                                         cGGAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGa                              >  2:1038444/1‑93 (MQ=255)

                                                           gAACGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGAcc                            <  2:2630378/93‑1 (MQ=255)

                                                            aaCGGTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCa                           >  1:408420/1‑93 (MQ=255)

                                                             aCGGTATCTGATGGCCACGGGATTGGTCCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAg                          <  2:3829838/93‑1 (MQ=255)

                                                               ggTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGtt                        <  1:3346698/93‑1 (MQ=255)

                                                               ggTATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGtt                        <  2:3827744/93‑1 (MQ=255)

                                                                 tATCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTg                      >  1:3522871/1‑93 (MQ=255)

                                                                  aTCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGc                     >  2:1629313/1‑93 (MQ=255)

                                                                   tCTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCg                    >  2:4514399/1‑93 (MQ=255)

                                                                    cTGATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCgt                   <  2:45996/93‑1 (MQ=255)

                                                                      gATGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGGCGAGGCGTAAAGCTGCGCCCACGGCGACCAGTTTGCgtcg                 >  2:3429170/1‑93 (MQ=255)

                                                                       aTGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCgtcgg                <  1:256054/93‑1 (MQ=255)

                                                                        tGGCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCgtcggt               >  1:3420605/1‑93 (MQ=255)

                                                                         ggCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCgtcggtc              >  2:4119623/1‑93 (MQ=255)

                                                                          gCCACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCgtcggtcg             >  2:4626855/1‑93 (MQ=255)

                                                                           ccACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCGTCGGTCGt            >  2:2912393/1‑93 (MQ=255)

                                                                            cACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCGTCGGTCGTa           <  2:2901924/93‑1 (MQ=255)

                                                                            cACGGGATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCGTCGGTCGTa           >  2:4341440/1‑93 (MQ=255)

                                                                               gggATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCGTCGGTCGTATcc        <  1:1738891/93‑1 (MQ=255)

                                                                               gggATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCGTCGGTCGTATcc        <  1:1397900/93‑1 (MQ=255)

                                                                               gggATTGGTGCACCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCGTCGGTCGTATcc        <  1:4132125/93‑1 (MQ=255)

                                                                                ggATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCGTCGGTCGTATccc       >  2:3431032/1‑93 (MQ=255)

                                                                                 gATTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCGTCGGTCGTATCCcg      >  1:2459189/1‑93 (MQ=255)

                                                                                  aTTGGTGCGCCGACTGGCTAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGCCGACCAGTTTGCGTCGGTCGTATCCcgt     >  2:1367047/1‑93 (MQ=255)

                                                                                  aTTGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCGTCGGTCGTATCCcgt     <  2:541400/93‑1 (MQ=255)

                                                                                   ttGGTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGGGTAAAGCTGCGCCCACGGCGACCAGTTTGCGTCGGTCGTATCCcgtc    >  2:1080735/1‑93 (MQ=255)

                                                                                     ggTGCGCCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCGTCGGTCGTATCCcgtcgt  >  1:2868656/1‑93 (MQ=255)

AAGTTTCGGGTATGCAGATTTGTCAAAAGTCTGCCCCTGCATCAGGGCATAACCAGACGGAACGGTATCTGATGGCCACGGGATTGGTGCACCGACTGGATAAAACTCTGCAGGAGGATGAGCCGAGGTGTAAAGCTGCGCCCACGGCGACCAGTTTGCGTCGGTCGTATCCCGTCGT  >  NC_000913/1633500‑1633677

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