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.5%	699.9 / NA	195	R177C (CGC→TGC)	pqqL	putative peptidase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (95/98); total (95/100)
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 199 total aligned reads displayed.

CCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCGCCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAC  >  NC_000913/1572608‑1572787

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

  gTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCAt                                                                                       >  2:6022446/1‑93 (MQ=255)

  gTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCAt                                                                                       >  1:759327/1‑93 (MQ=255)

  gTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGAGGCGCGCCTGAGAGGTGCACCAt                                                                                       >  1:2051990/1‑93 (MQ=255)

    ggCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCCCCAttt                                                                                     <  1:5721932/93‑1 (MQ=255)

    ggCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCAttt                                                                                     <  2:901079/93‑1 (MQ=255)

    ggCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCAttt                                                                                     >  2:5770106/1‑93 (MQ=255)

    ggCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCAttt                                                                                     <  1:5908497/93‑1 (MQ=255)

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

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

        aCAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCTAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCa                                                                                 >  2:4344758/1‑93 (MQ=255)

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

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

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

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

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

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

                  tCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGATGTGCACCATTTCGCATCCTGATgtg                                                                       >  2:5223497/1‑91 (MQ=255)

                  tCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATggg                                                                       <  2:3490479/93‑1 (MQ=255)

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

                    aGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCa                                                                     <  2:4537422/93‑1 (MQ=255)

                    aGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCa                                                                     >  1:4679554/1‑93 (MQ=255)

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

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

                        cGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTAGCATCCTGATGGGCACGcc                                                                 >  2:2412920/1‑93 (MQ=255)

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

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

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

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

                           tAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCAtt                                                              <  1:2379630/93‑1 (MQ=255)

                             ggTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTcc                                                            >  2:4034010/1‑93 (MQ=255)

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

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

                               ttCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTc                                                          >  1:1601628/1‑93 (MQ=255)

                               ttCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTc                                                          >  2:2916249/1‑93 (MQ=255)

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

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

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

                                 cACGGTCTAAATGACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCGCCATTTCGCATCCTGATGGGCACGCCATTCCTCAg                                                        <  2:1458479/93‑1 (MQ=255)

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

                                    ggTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTaa                                                     <  1:216245/93‑1 (MQ=255)

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

                                     gTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAAt                                                    >  1:4819102/1‑93 (MQ=255)

                                       cTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTa                                                  <  1:4027002/93‑1 (MQ=255)

                                       cTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTa                                                  <  2:629142/93‑1 (MQ=255)

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

                                         aaaTTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACg                                                <  1:1064965/93‑1 (MQ=255)

                                          aaTTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGCAATTACGc                                               <  2:6162598/93‑1 (MQ=255)

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

                                            ttACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCa                                             <  2:1271399/93‑1 (MQ=255)

                                            ttACGGGTATGGGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGGGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAAGTACGCCa                                             <  2:6254891/93‑1 (MQ=255)

                                               cGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGt                                          >  1:4197440/1‑93 (MQ=255)

                                                 ggTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTc                                        <  2:1018711/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                                        gCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCt                                 >  2:6047211/1‑93 (MQ=255)

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

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

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

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

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

                                                           agcagGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACt                              <  1:875651/93‑1 (MQ=255)

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

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

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

                                                               ggAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCTCGTTCAGCGTCTACTTCGa                          <  2:6118635/93‑1 (MQ=255)

                                                               ggAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGa                          <  2:2566435/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                   aGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGttt                      <  2:3551380/93‑1 (MQ=255)

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

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

                                                                    aggcggcGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGtttt                     <  2:2481421/92‑1 (MQ=255)

                                                                      gcggcgCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTc                   <  2:1065577/93‑1 (MQ=255)

                                                                       cggcgCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCa                  >  2:1724565/1‑93 (MQ=255)

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

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

                                                                           gcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAgg              >  1:5547862/1‑93 (MQ=255)

                                                                           gcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACACCACGTTCAGCGTCTACTTCGAGTTTTTCAAAgg              >  2:2818801/1‑93 (MQ=255)

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

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

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

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

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

                                                                                tGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCg         <  2:4231746/93‑1 (MQ=255)

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

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

                                                                                      gTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTa   <  1:3241937/93‑1 (MQ=255)

                                                                                      gTGCACCATTTCGCATCCTGATGGGAACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTa   <  1:6063715/93‑1 (MQ=255)

                                                                                       tGCACCATTGCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAc  <  2:4439441/93‑1 (MQ=255)

CCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCGCCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAC  >  NC_000913/1572608‑1572787

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