BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,224,926	+C	coding (1052/1434 nt)	ygjI →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,224,924	1	.	C	99.7%	1085.1 / ‑4.7	341	coding (1050/1434 nt)	ygjI	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (150/190); total (150/191)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 363 total aligned reads displayed.

TTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTAT‑CCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTATTC  >  NC_000913/3224832‑3225011

  tGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTAtacc                                                                                       >  2:4496392/1‑90 (MQ=255)

   ggTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTAT‑cccg                                                                                     >  1:4013604/1‑91 (MQ=255)

    gTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTAT‑cccga                                                                                    <  1:2341376/93‑4 (MQ=255)

        gAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCt                                                                                 >  2:1449825/1‑93 (MQ=255)

         aaaaCCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTc                                                                                >  2:3051926/1‑93 (MQ=255)

           aaCCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCgg                                                                              >  1:3517955/1‑93 (MQ=255)

            aCCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGt                                                                             >  2:1003729/1‑93 (MQ=255)

             ccGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGtt                                                                            >  1:6388806/1‑93 (MQ=255)

              cGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTc                                                                           <  1:6193662/93‑1 (MQ=255)

               gTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTcc                                                                          <  2:4881487/93‑1 (MQ=255)

               gTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTcc                                                                          <  2:7698823/93‑1 (MQ=255)

               gTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTcc                                                                          >  2:7825793/1‑93 (MQ=255)

                 cGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCaa                                                                        <  1:2842039/93‑1 (MQ=255)

                    aCTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATAc                                                                     <  1:5459730/93‑1 (MQ=255)

                    aCTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATAc                                                                     <  2:2992277/93‑1 (MQ=255)

                    aCTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATAc                                                                     <  1:6163080/93‑1 (MQ=255)

                    aCTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGAAGATTATCCCGATGCTCGGTTCCAATAc                                                                     <  2:7520908/93‑1 (MQ=255)

                     cTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATAcc                                                                    <  1:6026079/93‑1 (MQ=255)

                      agAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTTATGATTATCCCGATGCTCGGTTCCAATACCg                                                                   <  2:6723804/92‑1 (MQ=255)

                       gAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCgg                                                                  >  2:7193846/1‑92 (MQ=255)

                       gAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGt                                                                  >  2:1130777/1‑93 (MQ=255)

                         aCGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGc                                                                >  1:3726279/1‑93 (MQ=255)

                         aCGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGc                                                                <  1:4509534/93‑1 (MQ=255)

                          cGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCa                                                               <  2:5109574/93‑1 (MQ=255)

                           gAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCCTCCCGCTGATGCTTATCCCGATGCTCGCTTCCCATACCGTGCAg                                                              >  1:6971652/1‑93 (MQ=255)

                           gAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAg                                                              <  1:633624/93‑1 (MQ=255)

                           gAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAg                                                              <  1:6064702/93‑1 (MQ=255)

                            aaaaCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAgg                                                             >  1:2527955/1‑93 (MQ=255)

                             aaaCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGa                                                            <  2:3183458/93‑1 (MQ=255)

                              aaCGGCGTTCCGGCGC‑GCGGAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGAt                                                           >  1:7664498/1‑93 (MQ=255)

                              aaCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGAt                                                           >  2:1482497/1‑93 (MQ=255)

                               aCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATc                                                          <  2:4397339/93‑1 (MQ=255)

                               aCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATc                                                          >  2:1297092/1‑93 (MQ=255)

                               aCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATc                                                          <  2:1722104/93‑1 (MQ=255)

                                 ggCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCtg                                                        >  1:100916/1‑93 (MQ=255)

                                 ggCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCtg                                                        <  2:4749222/93‑1 (MQ=255)

                                   cGTGCCGGCGC‑GCTCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTGGGTTCCAATGCCGTGCAGGATCtgat                                                      <  2:279330/93‑1 (MQ=255)

                                    gTTCCGGCGC‑GCGCAGCGTTGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCtgatg                                                     <  1:5647937/93‑1 (MQ=255)

                                    gTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCtgatg                                                     <  2:3941375/93‑1 (MQ=255)

                                    gTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCGCGGTTGCAATACCGTGGAGGATGtgatg                                                     >  1:1729269/1‑93 (MQ=255)

                                     ttCCGGCGC‑GCGCAGCGTGGATCCGGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCtgatga                                                    <  2:3822976/93‑1 (MQ=255)

                                     ttCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCtgatga                                                    >  2:6787022/1‑93 (MQ=255)

                                     ttCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCtgatga                                                    >  2:842897/1‑93 (MQ=255)

                                      tCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGaa                                                   <  1:1104455/93‑1 (MQ=255)

                                      tCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGaa                                                   >  1:3228944/1‑93 (MQ=255)

                                       ccGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAAt                                                  >  2:3286265/1‑93 (MQ=255)

                                        cggcgc‑gcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATa                                                 <  2:4649862/93‑1 (MQ=255)

                                         ggcgc‑gcgcaGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATATCGTGCAGGATCTGATGAATat                                                >  2:5812782/1‑92 (MQ=255)

                                           gcgcgcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACt                                               <  2:6140359/89‑1 (MQ=255)

                                           gcgcgcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACt                                               <  2:256745/89‑1 (MQ=255)

                                            cgcgcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACta                                              <  1:1342418/90‑1 (MQ=255)

                                            cgcgcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACta                                              <  2:5738208/90‑1 (MQ=255)

                                            cgcgcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACta                                              <  1:4953137/90‑1 (MQ=255)

                                            cgcgcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACta                                              <  2:5392446/90‑1 (MQ=255)

                                              cgcgcAGCGTGGATCCAGTGCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACtatt                                            <  1:1932467/92‑1 (MQ=255)

                                                 gcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGCTGCTCGGTTCCCATCCCGTGCAGGATCTGATGAATACTattatt                                         >  1:144080/1‑93 (MQ=255)

                                                 gcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTattatt                                         <  2:5886380/93‑1 (MQ=255)

                                                   aGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTaa                                       <  2:2040892/93‑1 (MQ=255)

                                                   aGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTaa                                       >  2:2381396/1‑93 (MQ=255)

                                                    gCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAat                                      <  1:4118571/93‑1 (MQ=255)

                                                    gCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAat                                      >  2:3523842/1‑93 (MQ=255)

                                                    gCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAat                                      <  1:6593267/93‑1 (MQ=255)

                                                    gCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAat                                      >  2:6468333/1‑93 (MQ=255)

                                                      gTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAatat                                    <  2:5604929/93‑1 (MQ=255)

                                                      gTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAatat                                    <  1:3972055/93‑1 (MQ=255)

                                                       tGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATg                                   >  1:1904076/1‑93 (MQ=255)

                                                         gATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGAc                                 >  1:7436461/1‑93 (MQ=255)

                                                         gATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGAc                                 <  2:2611179/93‑1 (MQ=255)

                                                          aTCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGAcc                                <  1:4265190/93‑1 (MQ=255)

                                                          aTCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGAcc                                <  1:6287582/93‑1 (MQ=255)

                                                           tCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGAccg                               >  2:7322572/1‑93 (MQ=255)

                                                            ccAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGAccgc                              >  2:406770/1‑93 (MQ=255)

                                                             cAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGAccgcc                             <  1:300304/93‑1 (MQ=255)

                                                              aGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGAccgccg                            >  2:6532272/1‑93 (MQ=255)

                                                               gTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACcgccgc                           >  1:904098/1‑93 (MQ=255)

                                                               gTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACcgccgc                           <  1:3778318/93‑1 (MQ=255)

                                                                  ccTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGc                        <  2:4967815/93‑1 (MQ=255)

                                                                  ccTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGc                        <  2:5561240/93‑1 (MQ=255)

                                                                    tGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGt                      >  1:850100/1‑93 (MQ=255)

                                                                       tCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCa                   <  1:4550470/93‑1 (MQ=255)

                                                                       tCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCa                   <  1:501686/93‑1 (MQ=255)

                                                                         gTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATg                 >  2:4199454/1‑93 (MQ=255)

                                                                          tCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGc                >  2:1958889/1‑93 (MQ=255)

                                                                          tCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGc                >  1:4316005/1‑93 (MQ=255)

                                                                            aTCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCtt              >  2:4102001/1‑93 (MQ=255)

                                                                             tCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTc             >  1:2102124/1‑93 (MQ=255)

                                                                             tCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTc             >  1:7283195/1‑93 (MQ=255)

                                                                             tCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTc             <  1:7868325/93‑1 (MQ=255)

                                                                             tCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTc             <  1:1663493/93‑1 (MQ=255)

                                                                              cccGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGGCCGCCGCAGCGTCCATGCTTcc            <  2:603886/93‑1 (MQ=255)

                                                                                 gCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTccgcc         >  2:4286855/1‑93 (MQ=255)

                                                                                  cTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTccgccg        >  1:2761206/1‑93 (MQ=255)

                                                                                   tgatgaTTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGt       <  1:6751004/93‑1 (MQ=255)

                                                                                   tgatgaTTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGt       >  2:3737957/1‑93 (MQ=255)

                                                                                   tgatgaTTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGt       >  2:7370410/1‑93 (MQ=255)

                                                                                   tgatgaTTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGt       >  1:5222934/1‑93 (MQ=255)

                                                                                    gatgatTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGtt      <  1:6824241/93‑1 (MQ=255)

                                                                                      tgatTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTAt    >  1:5835132/1‑93 (MQ=255)

                                                                                      tgatTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTAt    <  2:2147348/93‑1 (MQ=255)

                                                                                        atTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTATtc  >  1:3318783/1‑93 (MQ=255)

TTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTAT‑CCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTATTC  >  NC_000913/3224832‑3225011

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: `‑`