BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,115,116	+C	coding (2013/2028 nt)	b2974 ←	putative endoglucanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,115,115	1	.	C	99.7%	1116.4 / ‑4.1	319	coding (2014/2028 nt)	b2974	putative endoglucanase
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (163/155); total (164/155)
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.52e-01

Only 100 of 363 total aligned reads displayed.

GTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑CTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTCC  >  NC_000913/3115023‑3115202

gTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAttt                                                                                          >  2:4181773/1‑93 (MQ=255)

gTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAttt                                                                                          >  1:8690591/1‑93 (MQ=255)

gTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAttt                                                                                          <  2:6759822/93‑1 (MQ=255)

 tttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑c                                                                                        >  1:1042260/1‑93 (MQ=255)

  ttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cc                                                                                       >  1:8096092/1‑92 (MQ=255)

  ttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cc                                                                                       >  1:2046927/1‑92 (MQ=255)

   tGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cct                                                                                      >  1:4552266/1‑91 (MQ=255)

    gCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cctg                                                                                     >  2:1909684/1‑90 (MQ=255)

     cTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cctgc                                                                                    <  2:3370717/93‑5 (MQ=255)

       ggTATCTGAACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCgg                                                                                   <  1:4816418/93‑1 (MQ=255)

        gTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGa                                                                                  <  2:2652270/93‑1 (MQ=255)

          aTCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGAtg                                                                                >  1:2599328/1‑92 (MQ=255)

          aTCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATc                                                                                <  2:7440584/93‑1 (MQ=255)

           tCTGCACGCAGCGGGATTGCATACGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCg                                                                               <  2:4747001/93‑1 (MQ=255)

            cTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCga                                                                              <  2:2193317/93‑1 (MQ=255)

            cTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCga                                                                              <  1:90696/93‑1 (MQ=255)

            cTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCga                                                                              <  1:7195234/93‑1 (MQ=255)

             tGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCgag                                                                             >  2:3758256/1‑93 (MQ=255)

              gCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCgaga                                                                            <  2:7940471/93‑1 (MQ=255)

               cACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGAt                                                                           <  2:8290920/93‑1 (MQ=255)

               cACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGAt                                                                           <  1:657093/93‑1 (MQ=255)

                aCGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATc                                                                          <  1:4974748/93‑1 (MQ=255)

                    aGCGGGATTGCATAAAGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCa                                                                      >  2:7329933/1‑93 (MQ=255)

                     gCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCAt                                                                     >  1:7947962/1‑93 (MQ=255)

                      cGGGATTTCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATa                                                                    <  2:3251271/93‑1 (MQ=255)

                      cGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATa                                                                    <  1:58247/93‑1 (MQ=255)

                       gggAGTGCATAAGGATCGTTACCCCCCTGAGTCACATATTCAAAGCCGTTAAGGAGCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAg                                                                   <  1:6773351/93‑1 (MQ=255)

                        ggATTGCATAAGGATCGTTACCCACCTGAGTTACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGc                                                                  >  1:2471104/1‑93 (MQ=255)

                        ggATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGc                                                                  >  2:2265834/1‑93 (MQ=255)

                        ggATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGc                                                                  <  1:3669137/93‑1 (MQ=255)

                           ttGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGAGGATCGAGATTGCCATAGCTAc                                                               >  1:3419118/1‑93 (MQ=255)

                               aTAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTaa                                                           >  2:508940/1‑93 (MQ=255)

                               aTAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTaa                                                           >  2:4999928/1‑93 (MQ=255)

                                tAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTaaa                                                          <  2:6049391/93‑1 (MQ=255)

                                 aaGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAAt                                                         <  2:4620236/93‑1 (MQ=255)

                                 aaGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAAt                                                         <  1:618442/93‑1 (MQ=255)

                                  aGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATg                                                        <  1:2285475/93‑1 (MQ=255)

                                   ggATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGc                                                       <  2:2880089/93‑1 (MQ=255)

                                   ggATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGc                                                       <  2:6253854/93‑1 (MQ=255)

                                    gATCGTTACCCACCTGAGTCTCATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCt                                                      <  2:7557150/93‑1 (MQ=255)

                                    gATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCt                                                      >  2:5613736/1‑93 (MQ=255)

                                      tCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCa                                                    >  2:4037219/1‑93 (MQ=255)

                                       cGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCaa                                                   <  1:5264612/93‑1 (MQ=255)

                                       cGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCaa                                                   >  2:426748/1‑93 (MQ=255)

                                        gTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAAc                                                  <  1:2115803/93‑1 (MQ=255)

                                        gTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAAc                                                  <  2:4526326/93‑1 (MQ=255)

                                        gTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAAc                                                  >  2:6138728/1‑93 (MQ=255)

                                          tACCCACCTGAGTCACGTATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACag                                                <  1:4658075/93‑1 (MQ=255)

                                          tACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACag                                                >  2:2524583/1‑93 (MQ=255)

                                           aCCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACaga                                               <  1:6715766/93‑1 (MQ=255)

                                             ccACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAgaga                                             <  1:299714/93‑1 (MQ=255)

                                             ccACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAgaga                                             >  2:6409182/1‑93 (MQ=255)

                                             ccACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAgaga                                             <  2:5513740/93‑1 (MQ=255)

                                             ccACCTGAGACACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAgaga                                             <  2:8037877/93‑1 (MQ=255)

                                               aCCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATg                                           >  1:7866188/1‑93 (MQ=255)

                                               aCCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATg                                           >  1:5983187/1‑93 (MQ=255)

                                                ccTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGGTGc                                          >  2:6474396/1‑93 (MQ=255)

                                                ccTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGc                                          >  1:1662365/1‑93 (MQ=255)

                                                  tGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATcccc                                        >  1:7015246/1‑93 (MQ=255)

                                                   gAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTCCTTAAATGCTCAACAGAGATGCCCg                                       <  1:8390898/93‑1 (MQ=255)

                                                   gAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCg                                       >  1:6339416/1‑93 (MQ=255)

                                                    aGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGc                                      >  1:3568758/1‑93 (MQ=255)

                                                      tCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaa                                    <  2:524087/93‑1 (MQ=255)

                                                      tCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaa                                    <  1:5688381/93‑1 (MQ=255)

                                                         caTATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAAt                                 >  2:1043200/1‑93 (MQ=255)

                                                          atatTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTTCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATc                                >  2:4939246/1‑93 (MQ=255)

                                                          atatTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATc                                >  2:5805522/1‑93 (MQ=255)

                                                               cAAAGCCGTTAAGGATCAGCATCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGAt                           <  2:5694073/93‑1 (MQ=255)

                                                               cAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGAt                           <  2:3184134/93‑1 (MQ=255)

                                                               cAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGAt                           <  2:2813027/93‑1 (MQ=255)

                                                                aaaGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATg                          >  2:4417532/1‑93 (MQ=255)

                                                                 aaGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATgg                         >  1:3264272/1‑93 (MQ=255)

                                                                 aaGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATgg                         >  1:149114/1‑93 (MQ=255)

                                                                    ccGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAg                      >  2:6948641/1‑93 (MQ=255)

                                                                    ccGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAg                      >  1:5525156/1‑93 (MQ=255)

                                                                     cGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGt                     >  1:402940/1‑93 (MQ=255)

                                                                      gTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTc                    <  2:8441926/93‑1 (MQ=255)

                                                                         aaGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGaa                 <  1:2342269/93‑1 (MQ=255)

                                                                         aaGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGaa                 >  2:782103/1‑93 (MQ=255)

                                                                          aGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAAc                >  2:1407956/1‑93 (MQ=255)

                                                                           ggATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACg               >  1:2164996/1‑93 (MQ=255)

                                                                            gATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGc              >  1:7380127/1‑93 (MQ=255)

                                                                            gATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGc              <  1:443548/93‑1 (MQ=255)

                                                                              tCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAg            >  1:133279/1‑93 (MQ=255)

                                                                              tCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAg            >  2:3527007/1‑93 (MQ=255)

                                                                              tCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAg            <  2:2046359/93‑1 (MQ=255)

                                                                                agcagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgcg          >  1:565664/1‑93 (MQ=255)

                                                                                 gcagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgcgc         <  1:4269474/93‑1 (MQ=255)

                                                                                   agcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTg       >  2:689939/1‑93 (MQ=255)

                                                                                      ggCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACGGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGttt    <  1:8595522/93‑1 (MQ=255)

                                                                                      ggCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGttt    >  2:5377006/1‑93 (MQ=255)

                                                                                      ggCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGttt    <  2:4870850/93‑1 (MQ=255)

                                                                                      ggCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGttt    >  1:1925853/1‑93 (MQ=255)

                                                                                       gCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTc   >  1:7477182/1‑93 (MQ=255)

                                                                                       gCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTc   >  1:4488927/1‑93 (MQ=255)

                                                                                       gCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTc   >  2:8616715/1‑93 (MQ=255)

                                                                                        cATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTcc  <  2:166746/93‑1 (MQ=255)

                                                                                        cATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTcc  <  1:6814928/93‑1 (MQ=255)

                                                                                        cATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTcc  <  1:5166674/93‑1 (MQ=255)

                                                                                        cATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTcc  <  1:374915/93‑1 (MQ=255)

GTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑CTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTCC  >  NC_000913/3115023‑3115202

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