BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,449,738	+G	coding (1106/1254 nt)	ytfR →	YtfR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,449,737	1	.	G	97.2%	379.2 / 4.8	145	coding (1105/1254 nt)	ytfR	YtfR
Reads supporting (aligned to +/- strand): ref base . (2/2); new base G (69/72); total (71/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 154 total aligned reads displayed.

AGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCT‑G‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCATTCGC  >  NC_000913/4449662‑4449807

aGGAAAACCGGTGGCGTCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGgct                                                                           <  2:1903787/76‑1 (MQ=255)

 ggAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGgct‑g                                                                         >  2:3246413/1‑76 (MQ=255)

 ggAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGgct‑g                                                                         >  1:1913969/1‑76 (MQ=255)

 gaaaaaGCGGTGGCATCACCGCGTCCGCCTCGGGGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGgct‑g                                                                         <  2:5725682/74‑1 (MQ=255)

  gAAAACCGATGGCATCATCGCGGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGgct‑gg                                                                        <  2:5544014/76‑2 (MQ=255)

   aaaaCCGATGGCGTCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGgct‑ggc                                                                       <  2:1361181/76‑3 (MQ=255)

       ccGATGGCATCATCGCTGCCGCCTCGGCGCGGGAAAATATCATCCTCGCTCCCCACGCCCCGCGCCGCT‑GACTACg                                                                   >  2:1656852/1‑76 (MQ=255)

        cGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGt                                                                  <  2:4024860/76‑1 (MQ=255)

            ggCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGcc‑ccc                                                             >  2:1798397/1‑75 (MQ=255)

             gCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCACGCCCCGCGCGGCTGG‑CTACGGC‑CCAt                                                            >  2:751539/1‑76 (MQ=255)

              catcatCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCAtt                                                           <  2:1046461/76‑1 (MQ=255)

                tcatcGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCGGCGCGGCTGG‑CTACGTC‑CCATTTc                                                         <  2:1288204/76‑1 (MQ=255)

                tcatcGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CCACGGC‑CCCTTTc                                                         >  2:4938243/1‑76 (MQ=255)

                  ctcGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTccc                                                       <  1:961411/75‑1 (MQ=255)

                  atcGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCCTTTccc                                                       >  2:3607657/1‑76 (MQ=255)

                  atcGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTccc                                                       <  1:3319174/76‑1 (MQ=255)

                    cGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGc                                                     >  1:163178/1‑76 (MQ=255)

                     gcggccgccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCa                                                    <  2:5239403/73‑1 (MQ=255)

                     gCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCa                                                    >  1:1013684/1‑76 (MQ=255)

                     gCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCa                                                    <  1:1586355/76‑1 (MQ=255)

                     gCTGCCGCCTCGGTGCGGGAAAATATCAGACTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTA‑CCATTTTCCGCa                                                    >  1:3267460/1‑76 (MQ=255)

                       tGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCcaa                                                  >  2:2603249/1‑76 (MQ=255)

                       ggccgccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCaaa                                                  <  2:2749429/75‑1 (MQ=255)

                       cgccgcgTCGGGGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCaaa                                                  <  2:575168/75‑1 (MQ=255)

                        gccgccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAg                                                 >  2:623551/1‑76 (MQ=255)

                         ccgccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTCCCCGCAAAGa                                                >  2:5089724/1‑76 (MQ=255)

                          cgccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGaa                                               <  1:4858106/76‑1 (MQ=255)

                          cgccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGaa                                               >  2:1242934/1‑76 (MQ=255)

                          cgccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGaa                                               <  2:4717909/76‑1 (MQ=255)

                          cgccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGaa                                               <  1:4040895/76‑1 (MQ=255)

                           gccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGaaa                                              >  1:5773737/1‑75 (MQ=255)

                           gccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAAc                                              <  2:5398995/76‑1 (MQ=255)

                           gccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAAc                                              <  1:3675968/76‑1 (MQ=255)

                           gccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAAc                                              <  2:2347302/76‑1 (MQ=255)

                           gacTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCCAAGAAc                                              <  2:3929508/74‑1 (MQ=255)

                            ccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACa                                             >  2:483860/1‑76 (MQ=255)

                            ccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACa                                             <  2:1313770/76‑1 (MQ=255)

                            ccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACa                                             >  1:5788147/1‑76 (MQ=255)

                            ccTCGGTGCGGGAAAATATCATCCTCGCTATCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACa                                             >  2:2037735/1‑76 (MQ=255)

                              tCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGc                                           >  2:821708/1‑76 (MQ=255)

                              tCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGc                                           <  1:2980604/76‑1 (MQ=255)

                               cGGTGCGGGAAAATATCATCCTCGGTCGCCAGGCCCAGCGCGGATGG‑CGACGGT‑GCATGTCCCCCAAAGAacaaca                                          >  1:472723/1‑76 (MQ=255)

                                ggTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTCGCCATTT‑CCGCTAAGAACAGCaa                                         >  1:1134624/1‑76 (MQ=255)

                                ggTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCCCAGAACAGCac                                         >  2:2939170/1‑75 (MQ=255)

                                ggTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCaa                                         <  1:862388/76‑1 (MQ=255)

                                 gTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAag                                        <  1:1834083/76‑1 (MQ=255)

                                 gTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAag                                        <  1:4338691/76‑1 (MQ=255)

                                  tGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAaga                                       >  1:2298655/1‑76 (MQ=255)

                                  tGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAaga                                       >  1:4850229/1‑76 (MQ=255)

                                    cGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTNCCCGCAAAGAACAGCAAgaga                                     <  1:2948740/76‑1 (MQ=255)

                                     gggAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAgagan                                    <  2:3350243/76‑2 (MQ=255)

                                      ggAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGAtt                                   <  2:1476861/76‑1 (MQ=255)

                                      ggAAAATATCATACTAGCTCTCCAGGCCCAGCGCGGCTGG‑ATACGTC‑CCGATTCCCGCAAAGAACAGCAAGACAta                                   >  1:3382028/1‑75 (MQ=255)

                                        aaaaTATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGc                                 >  1:3784631/1‑76 (MQ=255)

                                         aaaTATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGcc                                <  1:1144592/76‑1 (MQ=255)

                                         aaaTATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGcc                                >  1:4655803/1‑76 (MQ=255)

                                           atagCATCCTCGCTCGCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGa                              <  1:822263/76‑1 (MQ=255)

                                            tatCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGaa                             >  2:3417841/1‑76 (MQ=255)

                                            tatCATCCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGaa                             <  2:743225/76‑1 (MQ=255)

                                               catcCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGc                          <  2:3727889/76‑1 (MQ=255)

                                                atcCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCt                         >  2:4304800/1‑76 (MQ=255)

                                                 tcCTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCtt                        <  1:3159805/76‑1 (MQ=255)

                                                 accTCGCTCGCCAGGCCCAGCGCTGCT‑GTCTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGGTTGCCGAACGCtt                        <  2:2532640/75‑1 (MQ=255)

                                                  ccTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGCACGCgtt                       >  2:2442922/1‑76 (MQ=255)

                                                  ccTCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCttt                       <  1:732836/76‑1 (MQ=255)

                                                    tCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTAt                     >  1:4380982/1‑76 (MQ=255)

                                                    tCGCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTAt                     >  1:2765801/1‑76 (MQ=255)

                                                      gCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATcc                   >  1:3087249/1‑76 (MQ=255)

                                                      gCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATcc                   <  1:3454283/76‑1 (MQ=255)

                                                      gCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATcc                   <  1:2088515/76‑1 (MQ=255)

                                                      gCTCTCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATcc                   >  2:3261295/1‑76 (MQ=255)

                                                      cctctCCAGGCCCAGCGCGGCTGG‑CTACTTC‑CCATTTCCCGCAAAGAACAGCACGAGATTGCCGCACGCTTTATcc                   <  1:2803243/75‑1 (MQ=255)

                                                       ctctCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATccc                  >  2:1199658/1‑75 (MQ=255)

                                                       ctctCCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCg                  <  2:1589083/76‑1 (MQ=255)

                                                        tctcCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGc                 >  1:1692931/1‑76 (MQ=255)

                                                        tctcCAGGCCCAGCGAGGCTGG‑CTACGTC‑CCATTTACCGAAAAGAACAGCAAGAGATTGCCGCACGCTTTATTCGc                 >  1:2188006/1‑76 (MQ=255)

                                                        tcgccAGGCCCAGCGCGGCTGG‑CTAAGTC‑CCATTTCCCGCAAAGAACAGCAAGAGAGTGCCGAACGCTTTAGCCGc                 >  2:3203705/4‑76 (MQ=255)

                                                        cctcCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGc                 <  2:2976496/75‑1 (MQ=255)

                                                         ctcCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGcc                <  1:2760983/76‑1 (MQ=255)

                                                          tcCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCa               >  1:2054774/1‑76 (MQ=255)

                                                          tcCAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCa               >  1:1874859/1‑76 (MQ=255)

                                                           ccAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAg              <  1:5846922/76‑1 (MQ=255)

                                                           ccAGGCCCAGCGCGGCTGG‑CTACGGC‑CCATTTCCCGCAAAGAACAGCAAGAGAGTGCCGAACGCTTTATCCGCCAg              >  1:3824054/1‑76 (MQ=255)

                                                            cAGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGc             <  2:3274982/76‑1 (MQ=255)

                                                             aGGCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCt            >  1:5360720/1‑76 (MQ=255)

                                                             aGGCCCAGCGCGGCTGG‑CTACGGC‑CCACTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTCTCCGCCAGCt            <  2:603200/76‑1 (MQ=255)

                                                               gCCCAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTg          >  2:2758461/1‑76 (MQ=255)

                                                                cccAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTgg         <  1:1418410/76‑1 (MQ=255)

                                                                 ccAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTCGc        >  2:2181165/1‑76 (MQ=255)

                                                                 ccAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCCTGGc        >  2:3696294/1‑76 (MQ=255)

                                                                  cAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCa       >  2:4279985/1‑76 (MQ=255)

                                                                  cAGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCCTGGCa       >  2:3916866/1‑76 (MQ=255)

                                                                   aGCTCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCAt      <  1:2273009/76‑1 (MQ=255)

                                                                   aGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCAt      <  1:4030047/76‑1 (MQ=255)

                                                                   aGCGCGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCAt      <  1:2810547/76‑1 (MQ=255)

                                                                   aGCGCGGCT‑GCCTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCCGCTTGGcct      >  2:3348208/1‑74 (MQ=255)

                                                                    gcgcGGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCAtt     <  2:1706651/76‑1 (MQ=255)

                                                                     cgcgGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCATTc    <  2:2636311/76‑1 (MQ=255)

                                                                      gcgGCTGG‑CTACGTC‑CCATTTCCCGCAAAGCACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCATTCg   <  2:4027539/76‑1 (MQ=255)

                                                                      gcgGCTGG‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCATTCg   >  2:34957/1‑76 (MQ=255)

AGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCT‑G‑CTACGTC‑CCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCATTCGC  >  NC_000913/4449662‑4449807

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 6 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`