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	99.6%	783.6 / NA	225	R177C (CGC→TGC)	pqqL	putative peptidase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (104/120); total (104/121)
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 226 total aligned reads displayed.

CCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCGCCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAC  >  NC_000913/1572608‑1572787

ccGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCAcc                                                                                         >  1:6955014/1‑93 (MQ=255)

ccGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCAcc                                                                                         <  2:2508106/93‑1 (MQ=255)

  gTGGCGACAGTATCAATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCAt                                                                                       >  2:7406897/1‑93 (MQ=255)

       gACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGc                                                                                  >  1:7556624/1‑93 (MQ=255)

       gACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGCGCACCATTTCGc                                                                                  >  1:7842995/1‑93 (MQ=255)

        aCAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCa                                                                                 <  1:4635332/93‑1 (MQ=255)

        aCAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCa                                                                                 <  1:408811/93‑1 (MQ=255)

        aCAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCa                                                                                 <  2:6346925/93‑1 (MQ=255)

         cAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCAt                                                                                <  1:6117500/93‑1 (MQ=255)

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

           gTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATcc                                                                              <  1:4899609/93‑1 (MQ=255)

            tATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGTGAGGTGCACCATTTCGCATCCt                                                                             <  2:4983312/93‑1 (MQ=255)

             aTCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTg                                                                            >  2:616787/1‑93 (MQ=255)

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

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

                 aTCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATgg                                                                        <  1:6558799/93‑1 (MQ=255)

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

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

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

                          gtAGGTTCACGGTCTAAATTACGTGTATTTGCCAGCAGGAAATGGCTGCGCTCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCTAt                                                               <  2:6839180/92‑1 (MQ=255)

                          aTAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCAt                                                               >  1:1031229/1‑93 (MQ=255)

                          aTAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCAt                                                               >  2:4326880/1‑93 (MQ=255)

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

                            aGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTc                                                             >  1:723661/1‑93 (MQ=255)

                            aGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTc                                                             >  1:6237877/1‑93 (MQ=255)

                              gTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCt                                                           >  2:7670183/1‑93 (MQ=255)

                              gTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCt                                                           <  2:755750/93‑1 (MQ=255)

                              gTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCt                                                           >  1:2647064/1‑93 (MQ=255)

                              gTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCt                                                           >  1:1648037/1‑93 (MQ=255)

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

                                tCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCa                                                         <  1:311425/93‑1 (MQ=255)

                                 cACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAg                                                        >  2:1405306/1‑93 (MQ=255)

                                 cACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAg                                                        <  1:5956560/93‑1 (MQ=255)

                                 cACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAg                                                        <  1:1019248/93‑1 (MQ=255)

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

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

                                   cGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTa                                                      >  2:5502306/1‑93 (MQ=255)

                                   cGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTa                                                      >  2:7068491/1‑93 (MQ=255)

                                    ggTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGCGCACGCCATTCCTCAGTaa                                                     >  1:6323619/1‑93 (MQ=255)

                                     gTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAAt                                                    <  2:3151920/93‑1 (MQ=255)

                                      tCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAAtt                                                   >  1:434617/1‑93 (MQ=255)

                                       cTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTa                                                  >  2:6889025/1‑93 (MQ=255)

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

                                          aaTTACGGGTATTTGCCAGCAGGAAAGGGGGGCGCGCATGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGc                                               >  2:4837799/1‑93 (MQ=255)

                                          aaTTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGc                                               >  2:2827780/1‑93 (MQ=255)

                                          aaTTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGc                                               <  2:2281376/93‑1 (MQ=255)

                                           aTTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGcc                                              >  2:7852520/1‑93 (MQ=255)

                                           aTTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGcc                                              <  1:740010/93‑1 (MQ=255)

                                            ttACGGGTATTTGCCAGCAGCAAAGGGCGGCGCGCCTGAGACGTGCACCATTTCTCATCCTGATGGGCACGCCATTCCTCAGTAATTACGccc                                             >  2:3221259/1‑92 (MQ=255)

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

                                              aCGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACg                                           <  2:2379901/93‑1 (MQ=255)

                                               cGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTACTTACGCCACGt                                          <  2:3841106/93‑1 (MQ=255)

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

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

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

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

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

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

                                                    aTTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGc                                     <  1:2392805/93‑1 (MQ=255)

                                                    aTTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTAAGc                                     >  1:6259952/1‑93 (MQ=255)

                                                     tttGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGTCACGCCATTCCTCAGTAATTACGCCACGTTCAGCg                                    <  2:7830183/93‑1 (MQ=255)

                                                     tttGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCg                                    >  2:1281393/1‑93 (MQ=255)

                                                        gCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCt                                 <  2:3563567/93‑1 (MQ=255)

                                                         ccAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTa                                >  2:4930170/1‑93 (MQ=255)

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

                                                          cagcagGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTAc                               <  1:4347931/93‑1 (MQ=255)

                                                          cagcagGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTAc                               <  2:4706224/93‑1 (MQ=255)

                                                          cagcagGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTAc                               <  1:3893638/93‑1 (MQ=255)

                                                          cagcagGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTAc                               <  1:7735135/93‑1 (MQ=255)

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

                                                               ggATAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGa                          <  2:6623263/93‑1 (MQ=255)

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

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

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

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

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

                                                                  aaGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCCGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAggt                       >  2:1287641/1‑91 (MQ=255)

                                                                   aGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGttt                      >  1:6166468/1‑93 (MQ=255)

                                                                   aGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGttt                      >  1:531392/1‑93 (MQ=255)

                                                                    gggCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGtttt                     >  2:7522121/1‑93 (MQ=255)

                                                                    gggCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGtttt                     >  2:6113391/1‑93 (MQ=255)

                                                                     ggcggcGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGTGTCTACTTCGAGttttt                    <  1:4781654/93‑1 (MQ=255)

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

                                                                        ggcgCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCaa                 <  1:15234/93‑1 (MQ=255)

                                                                        ggcgCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCaa                 <  2:43707/93‑1 (MQ=255)

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

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

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

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

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

                                                                              ccTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTg           <  2:1871673/93‑1 (MQ=255)

                                                                              ccTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTg           <  2:3899145/93‑1 (MQ=255)

                                                                               cTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGc          <  1:7093387/93‑1 (MQ=255)

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

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

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

                                                                                   gagGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCa      >  1:4364884/1‑93 (MQ=255)

                                                                                   gagGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCa      >  1:2928697/1‑93 (MQ=255)

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

                                                                                     ggTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCAtt    <  1:3394124/93‑1 (MQ=255)

CCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCGCCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAC  >  NC_000913/1572608‑1572787

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