BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,750,491	T→C	Y455C (TAC→TGC)	b1673 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,750,491	0	T	C	99.6%	725.0 / ‑4.8	224	Y455C (TAC→TGC)	b1673	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (117/106); total (117/107)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

Only 100 of 229 total aligned reads displayed.

ACAATGGAACCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGTACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACGGTAG  >  NC_000913/1750401‑1750580

aCAATGGAACCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCAc                                                                                         <  2:4399540/93‑1 (MQ=255)

  aaTGGAACCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCActc                                                                                       >  1:3948833/1‑93 (MQ=255)

  aaTGGAACCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCActc                                                                                       >  2:2167181/1‑93 (MQ=255)

  aaTGGAACCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCActc                                                                                       <  2:315300/93‑1 (MQ=255)

  aaTGGAACCGACCTGCGCCGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCActc                                                                                       <  1:2748002/93‑1 (MQ=255)

   aTGGAACCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCActct                                                                                      >  2:5418198/1‑93 (MQ=255)

   aTGGAACCGACCTGCGCTGACCCTTAATTGGCAAGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCActct                                                                                      <  2:4429408/93‑1 (MQ=255)

     ggAACCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTc                                                                                    <  2:630870/93‑1 (MQ=255)

     ggAACCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTc                                                                                    <  2:3343526/93‑1 (MQ=255)

      gAACCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCa                                                                                   >  2:6657897/1‑93 (MQ=255)

       aaCCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑ATTGTTTTTCGCGTAGCCCCAGCACTCTTCAc                                                                                  >  1:4599667/1‑93 (MQ=255)

        aCCGACCTGCGCCGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCAcc                                                                                 >  1:6587761/1‑93 (MQ=255)

          cGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCa                                                                               >  1:1802565/1‑93 (MQ=255)

            aCCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCGCTCTTCACCCaaa                                                                             >  2:5873578/1‑93 (MQ=255)

            aCCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCaaa                                                                             >  2:3437197/1‑93 (MQ=255)

            aCCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCaaa                                                                             <  2:5259419/93‑1 (MQ=255)

            aCCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCaaa                                                                             >  2:810790/1‑93 (MQ=255)

            aCCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCaaa                                                                             >  1:1445947/1‑93 (MQ=255)

             ccTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAAt                                                                            <  2:1390481/93‑1 (MQ=255)

               tGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTc                                                                          >  1:5359820/1‑93 (MQ=255)

                   cTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGc                                                                      >  1:498798/1‑93 (MQ=255)

                   cTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGc                                                                      <  2:6468500/93‑1 (MQ=255)

                   cTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TGTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGc                                                                      <  2:3555864/93‑1 (MQ=255)

                    tGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCg                                                                     >  2:620049/1‑93 (MQ=255)

                       cGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTc                                                                  <  2:4294121/93‑1 (MQ=255)

                       cGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTc                                                                  <  2:4213743/93‑1 (MQ=255)

                       cGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTc                                                                  >  2:6631059/1‑93 (MQ=255)

                        gCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCt                                                                 <  1:4422345/93‑1 (MQ=255)

                        gCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCt                                                                 >  2:2294382/1‑93 (MQ=255)

                        gCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCt                                                                 >  1:307726/1‑93 (MQ=255)

                          ttCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTgg                                                               >  2:4958066/1‑92 (MQ=255)

                           tCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCg                                                              >  1:4121707/1‑93 (MQ=255)

                             aTTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGAt                                                            >  2:190323/1‑93 (MQ=255)

                              ttagCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATg                                                           >  2:2719552/4‑93 (MQ=255)

                              ttGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTCCTGCGATg                                                           <  1:1723393/93‑1 (MQ=255)

                              ttGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGGTCTGCGATg                                                           >  1:2931972/1‑93 (MQ=255)

                               tGGCATGGTGAATCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATgg                                                          <  2:5082108/93‑1 (MQ=255)

                               tGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATgg                                                          >  1:1027712/1‑93 (MQ=255)

                               tGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATgg                                                          <  1:5698497/93‑1 (MQ=255)

                                 gCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCa                                                        >  2:1514250/1‑93 (MQ=255)

                                 gCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCa                                                        <  2:4630209/93‑1 (MQ=255)

                                   aTGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGtata                                                      <  1:5550143/93‑1 (MQ=255)

                                   aTGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCata                                                      >  1:4971830/1‑93 (MQ=255)

                                     ggTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATg                                                    >  2:1581804/1‑93 (MQ=255)

                                     ggTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATg                                                    <  1:5114681/93‑1 (MQ=255)

                                      gTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGa                                                   <  1:2033376/93‑1 (MQ=255)

                                        gAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAAc                                                 <  1:1164548/93‑1 (MQ=255)

                                         aaCCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAAcc                                                >  2:6355773/1‑93 (MQ=255)

                                            cGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATc                                             <  2:319647/93‑1 (MQ=255)

                                            cGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATc                                             >  1:3815161/1‑93 (MQ=255)

                                             ggATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCa                                            <  2:1845324/93‑1 (MQ=255)

                                             ggATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCa                                            <  1:4164942/93‑1 (MQ=255)

                                              gATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAg                                           >  1:2403436/1‑93 (MQ=255)

                                               aTAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGc                                          <  2:2246867/93‑1 (MQ=255)

                                                tAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCCCTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGcc                                         <  1:3207147/93‑1 (MQ=255)

                                                tAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGcc                                         <  2:4686739/93‑1 (MQ=255)

                                                tAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGcc                                         <  1:3972949/93‑1 (MQ=255)

                                                 tagcAAACGGCCGCCCGAGTTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGcc                                         <  1:6153736/90‑1 (MQ=255)

                                                 aGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCa                                        <  1:3416989/93‑1 (MQ=255)

                                                 aGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCa                                        >  2:3892048/1‑93 (MQ=255)

                                                  gCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAg                                       <  1:240496/93‑1 (MQ=255)

                                                  gCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAg                                       >  1:6531632/1‑93 (MQ=255)

                                                   ccAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAgg                                      >  1:2160780/1‑93 (MQ=255)

                                                    cAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGt                                     <  1:2146158/93‑1 (MQ=255)

                                                     aaaCGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTg                                    >  2:5512960/1‑93 (MQ=255)

                                                     aaaCGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTg                                    >  2:4758731/1‑93 (MQ=255)

                                                       aCGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGAc                                  >  1:2248256/1‑93 (MQ=255)

                                                          gCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCa                               <  2:3672932/93‑1 (MQ=255)

                                                            gACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATTGCATATGAACCATCAGCCAGGTGACTCAGc                             >  2:2943671/1‑93 (MQ=255)

                                                            gACCAGAG‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGc                             >  2:4521330/1‑93 (MQ=255)

                                                                agag‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCAc                         >  1:4759628/1‑93 (MQ=255)

                                                                 gag‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCAcc                        >  1:5826812/1‑93 (MQ=255)

                                                                 gag‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCAcc                        >  2:2384351/1‑93 (MQ=255)

                                                                   g‑TTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCca                      >  2:5221944/1‑93 (MQ=255)

                                                                     tttATTTTTCGCGTAGCCCCAGTACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCcac                     <  1:1093103/93‑1 (MQ=255)

                                                                      ttATTTTTCGCGTAGCCGCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCcaca                    <  1:1994607/93‑1 (MQ=255)

                                                                      ttATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCcaca                    >  1:3301936/1‑93 (MQ=255)

                                                                       tATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCacac                   >  1:1368123/1‑93 (MQ=255)

                                                                        aTTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACg                  >  1:5386978/1‑93 (MQ=255)

                                                                         tttttCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACgc                 >  1:1526931/1‑92 (MQ=255)

                                                                         tttttCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGa                 <  1:3539367/93‑1 (MQ=255)

                                                                         tttttCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGa                 <  1:230422/93‑1 (MQ=255)

                                                                          ttttCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGAt                <  1:6016837/93‑1 (MQ=255)

                                                                          ttttCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGAt                <  1:2841199/93‑1 (MQ=255)

                                                                               gcgTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGcc           <  1:6426085/93‑1 (MQ=255)

                                                                               gcgTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGcc           >  2:5025073/1‑93 (MQ=255)

                                                                               gcgTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGcc           >  2:695858/1‑93 (MQ=255)

                                                                                 gTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAg         >  1:1747296/1‑93 (MQ=255)

                                                                                 gTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAg         >  1:3900840/1‑93 (MQ=255)

                                                                                 gTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAg         >  1:6659046/1‑93 (MQ=255)

                                                                                   aGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGAc       <  1:3772730/93‑1 (MQ=255)

                                                                                   aGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGAc       <  1:4341980/93‑1 (MQ=255)

                                                                                   aGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGAc       <  1:1889494/93‑1 (MQ=255)

                                                                                   aGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGAc       >  1:6269963/1‑93 (MQ=255)

                                                                                   aGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGAc       <  2:2060440/93‑1 (MQ=255)

                                                                                    gCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACg      <  2:6389201/93‑1 (MQ=255)

                                                                                     ccccAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACgg     <  1:906860/93‑1 (MQ=255)

                                                                                      cccAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACACAGCTCACCCACACGATGCGCCAGACGgt    <  2:3073618/93‑1 (MQ=255)

                                                                                       ccAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACGgta   >  2:6026694/1‑93 (MQ=255)

                                                                                       ccAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACGgta   <  2:3976133/93‑1 (MQ=255)

ACAATGGAACCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAG‑TTTATTTTTCGCGTAGCCCCAGTACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACGGTAG  >  NC_000913/1750401‑1750580

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