BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,572,698	G→A	R177C (CGC→TGC)	pqqL ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,572,698	0	G	A	98.7%	522.1 / ‑3.6	150	R177C (CGC→TGC)	pqqL	putative peptidase
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (69/79); total (70/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 159 total aligned reads displayed.

GTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCGCCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAC  >  NC_000913/1572610‑1572787

gTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTTAGAGGTGCACCAt                                                                                       >  2:3605349/1‑93 (MQ=255)

   gCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTc                                                                                    <  2:1577287/93‑1 (MQ=255)

   gCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTc                                                                                    >  2:2204697/1‑93 (MQ=255)

      aCAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCa                                                                                 >  2:3933626/1‑93 (MQ=255)

      aCAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCa                                                                                 >  2:1765915/1‑93 (MQ=255)

        aGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATc                                                                               <  2:709999/93‑1 (MQ=255)

         gTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATcc                                                                              >  1:2101910/1‑93 (MQ=255)

          tATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGCATTTGCCAGCAGGAAAGGGTGGCGCGCCTGAGATGTGCACCATTTCGCATCCt                                                                             >  2:2444958/1‑93 (MQ=255)

            tCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGa                                                                           <  2:97290/93‑1 (MQ=255)

            tCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGa                                                                           >  2:4391826/1‑93 (MQ=255)

            tCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGa                                                                           >  1:3827680/1‑93 (MQ=255)

             ccATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGag                                                                          >  2:1705772/1‑92 (MQ=255)

             ccATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGAt                                                                          <  1:260341/93‑1 (MQ=255)

             ccATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGAt                                                                          <  2:3057574/93‑1 (MQ=255)

             ccATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGAt                                                                          <  2:387374/93‑1 (MQ=255)

              cATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATg                                                                         <  2:421665/93‑1 (MQ=255)

              cATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATg                                                                         >  1:1094456/1‑93 (MQ=255)

               ctCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATgg                                                                        <  1:3604096/92‑1 (MQ=255)

               aTCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATgg                                                                        >  1:3965829/1‑93 (MQ=255)

               aTCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATgg                                                                        >  1:372816/1‑93 (MQ=255)

                tCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATggg                                                                       >  2:4152128/1‑93 (MQ=255)

                 cAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGc                                                                      <  2:2699799/93‑1 (MQ=255)

                   ggCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGGGCGCCTGAGAGGTGCACCATTTCGCATTCTGATGGGCAc                                                                    >  2:1293475/1‑93 (MQ=255)

                   ggCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCAc                                                                    >  2:1737184/1‑93 (MQ=255)

                   ggCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCAc                                                                    >  1:108617/1‑93 (MQ=255)

                   ggCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCAc                                                                    <  1:3718638/93‑1 (MQ=255)

                       gATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCa                                                                >  2:799346/1‑93 (MQ=255)

                        aTAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCAt                                                               <  1:1712824/93‑1 (MQ=255)

                         tAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCAtt                                                              >  1:414246/1‑93 (MQ=255)

                         tAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCAtt                                                              >  1:616341/1‑93 (MQ=255)

                             ttCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTc                                                          <  2:2251575/93‑1 (MQ=255)

                              tCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCa                                                         >  2:4308530/1‑93 (MQ=255)

                              tCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCa                                                         >  1:789466/1‑93 (MQ=255)

                                aCGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGt                                                       >  2:599655/1‑93 (MQ=255)

                                  ggTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTaa                                                     >  2:1373977/1‑93 (MQ=255)

                                    tCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAAtt                                                   <  1:4483434/93‑1 (MQ=255)

                                      tAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTAc                                                 >  2:835205/1‑93 (MQ=255)

                                      tAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTAc                                                 <  1:656060/93‑1 (MQ=255)

                                        aaTTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTATGc                                               <  2:3437717/93‑1 (MQ=255)

                                         aTTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCGCCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGcc                                              <  2:4040395/93‑1 (MQ=255)

                                          ttACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCa                                             >  1:1338002/1‑93 (MQ=255)

                                          atACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCa                                             <  2:3842052/92‑1 (MQ=255)

                                           tACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCAc                                            <  1:1890329/93‑1 (MQ=255)

                                            aCGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACg                                           >  2:2687403/1‑93 (MQ=255)

                                            aCGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACg                                           >  1:1853480/1‑93 (MQ=255)

                                             cGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGt                                          <  2:1063517/93‑1 (MQ=255)

                                             cGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGt                                          <  2:2122653/93‑1 (MQ=255)

                                             cGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGt                                          <  1:2427389/93‑1 (MQ=255)

                                              gggTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGtt                                         <  2:389199/93‑1 (MQ=255)

                                              gggTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGtt                                         <  2:4693269/93‑1 (MQ=255)

                                               ggTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTc                                        >  1:547784/1‑93 (MQ=255)

                                                gTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCTCGTTCa                                       >  2:4224787/1‑93 (MQ=255)

                                                gTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCa                                       >  2:203781/1‑93 (MQ=255)

                                                gTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCa                                       <  1:2910015/93‑1 (MQ=255)

                                                gTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCa                                       <  1:3070679/93‑1 (MQ=255)

                                                gTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCa                                       <  1:2910016/93‑1 (MQ=255)

                                                 tATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAg                                      >  2:815436/1‑93 (MQ=255)

                                                 tATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAg                                      <  1:1780219/93‑1 (MQ=255)

                                                     tGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTc                                  <  1:910425/93‑1 (MQ=255)

                                                      gCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCCCATCCTCATGGGCACGCCATTCCTAGGCAATTACGCCACATTCCGCGTCt                                 >  2:1222242/1‑93 (MQ=255)

                                                        cagcagGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTAc                               >  1:360073/1‑93 (MQ=255)

                                                         agcagGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACt                              >  1:540361/1‑93 (MQ=255)

                                                          gcagGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACtt                             <  1:1210013/93‑1 (MQ=255)

                                                           cagGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTc                            <  2:2838082/93‑1 (MQ=255)

                                                           cagGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTc                            <  1:1601807/93‑1 (MQ=255)

                                                           cagGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTc                            >  2:167943/1‑93 (MQ=255)

                                                            agGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCg                           >  2:1301402/1‑93 (MQ=255)

                                                             ggAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGa                          >  2:2329501/1‑93 (MQ=255)

                                                              gAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAg                         >  2:320556/1‑93 (MQ=255)

                                                               aaaGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGt                        <  2:1868026/93‑1 (MQ=255)

                                                               aaaGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGt                        >  2:2616647/1‑93 (MQ=255)

                                                                aaGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGtt                       <  2:3524317/93‑1 (MQ=255)

                                                                aaGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGtt                       >  2:3416214/1‑93 (MQ=255)

                                                                   ggcggcGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGttttt                    <  1:477604/93‑1 (MQ=255)

                                                                     cggcgCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCGCGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCa                  <  1:417159/93‑1 (MQ=255)

                                                                      ggcgCGCCTGAGAGGTGCACCACTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTCCGCCACGTTCAGCGTCTACTTCGAGTTTTTCaa                 <  1:3676369/93‑1 (MQ=255)

                                                                       gcgcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCaaa                <  1:756071/93‑1 (MQ=255)

                                                                       gcgcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCaaa                <  2:1810025/93‑1 (MQ=255)

                                                                       gcgcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCaaa                <  1:4434499/93‑1 (MQ=255)

                                                                       gcgcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCaaa                >  2:1124312/1‑93 (MQ=255)

                                                                        ggcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAg               <  2:859388/92‑1 (MQ=255)

                                                                        cgcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAg               <  1:2799989/93‑1 (MQ=255)

                                                                        cgcgcCTGAGAGGTGCACCATTTCGAATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAg               <  1:1772930/93‑1 (MQ=255)

                                                                         gcgcCTTAGATGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACTTGCAGCGTCTACTTCGAGTTTTTCAAAgg              <  1:1468921/93‑1 (MQ=255)

                                                                          cgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGt             <  1:3849200/93‑1 (MQ=255)

                                                                           gcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGtt            <  2:3620587/93‑1 (MQ=255)

                                                                           gcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGtt            >  1:1696302/1‑93 (MQ=255)

                                                                           gcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGtt            <  2:1155298/93‑1 (MQ=255)

                                                                           gcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGtt            >  1:1184175/1‑93 (MQ=255)

                                                                           gcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGtt            <  1:1974914/93‑1 (MQ=255)

                                                                           cccTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGtt            <  2:1510503/92‑1 (MQ=255)

                                                                             cTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGc          >  2:435021/1‑93 (MQ=255)

                                                                              tGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCg         >  1:2242864/1‑93 (MQ=255)

                                                                              tGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCg         >  2:2606406/1‑93 (MQ=255)

                                                                                agagGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGc       >  1:2740107/1‑93 (MQ=255)

                                                                                 gagGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCCAAGGTTGCGGCa      >  1:310205/1‑93 (MQ=255)

                                                                                  agGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCAt     <  1:585700/93‑1 (MQ=255)

                                                                                   ggTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCAtt    >  1:2694687/1‑93 (MQ=255)

                                                                                     tGCGCCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAc  <  1:318006/93‑1 (MQ=255)

                                                                                     tGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAc  <  2:3738878/93‑1 (MQ=255)

GTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCGCCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAC  >  NC_000913/1572610‑1572787

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