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.6%	913.8 / ‑3.6	259	coding (2014/2028 nt)	b2974	putative endoglucanase
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (117/141); total (117/142)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

Only 100 of 280 total aligned reads displayed.

GTTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑CTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTCC  >  NC_000913/3115023‑3115202

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

  ttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cc                                                                                       <  2:5810963/93‑2 (MQ=255)

  ttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cc                                                                                       <  2:3932329/93‑2 (MQ=255)

  ttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cc                                                                                       <  1:978347/93‑2 (MQ=255)

    gCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cctg                                                                                     <  2:3655306/93‑4 (MQ=255)

     cTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTT‑cctgc                                                                                    >  1:2912853/1‑89 (MQ=255)

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

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

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

         tATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGAt                                                                                 <  1:2207776/93‑1 (MQ=255)

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

           tCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCg                                                                               <  1:3135467/93‑1 (MQ=255)

           tCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCg                                                                               >  1:2572798/1‑93 (MQ=255)

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

             ggCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCgag                                                                             >  1:6201692/2‑93 (MQ=255)

              gCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCgaga                                                                            >  2:4201547/1‑93 (MQ=255)

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

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

                 cGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCg                                                                         >  1:3537990/1‑93 (MQ=255)

                  gCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGc                                                                        <  1:718829/93‑1 (MQ=255)

                    aGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCa                                                                      >  1:2163421/1‑93 (MQ=255)

                     gCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCAt                                                                     <  1:6236556/93‑1 (MQ=255)

                       gggATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAg                                                                   >  1:5432100/1‑93 (MQ=255)

                         gATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCt                                                                 <  2:4547782/93‑1 (MQ=255)

                          aTTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGACGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTa                                                                >  2:4529052/1‑93 (MQ=255)

                           ttGCATAAGGATCGTTACCCCCCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTAc                                                               <  1:5848510/93‑1 (MQ=255)

                           ttGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTAc                                                               <  2:2260124/93‑1 (MQ=255)

                           ttGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTAc                                                               <  2:6410014/93‑1 (MQ=255)

                           ttGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTAc                                                               >  2:892087/1‑93 (MQ=255)

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

                                 aaGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCAACTTAAAt                                                         <  1:4439343/93‑1 (MQ=255)

                                  aGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATg                                                        >  1:4115956/1‑93 (MQ=255)

                                  aGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATg                                                        >  2:3026422/1‑93 (MQ=255)

                                   ggATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGc                                                       >  1:4083511/1‑93 (MQ=255)

                                    gATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCt                                                      <  2:1632996/93‑1 (MQ=255)

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

                                     aTCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTc                                                     <  2:4594429/93‑1 (MQ=255)

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

                                        gTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGCGATCGCCATAGCTACTTAAATGCTCAAc                                                  >  2:1949893/1‑93 (MQ=255)

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

                                           aCCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACaga                                               >  1:3498431/1‑93 (MQ=255)

                                            cccACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACagag                                              >  1:2597446/1‑93 (MQ=255)

                                              cACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGAt                                            <  2:2435092/93‑1 (MQ=255)

                                              cACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGAt                                            >  2:6117834/1‑93 (MQ=255)

                                              cACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGAt                                            <  1:2775546/93‑1 (MQ=255)

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

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

                                                 cTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGcc                                         <  1:1648029/93‑1 (MQ=255)

                                                 cTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGcc                                         <  2:3349772/93‑1 (MQ=255)

                                                  tGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGccc                                        <  2:3144475/93‑1 (MQ=255)

                                                  tGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGccc                                        <  1:1565282/93‑1 (MQ=255)

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

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

                                                     gTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCa                                     >  1:1453727/1‑93 (MQ=255)

                                                       cacaTATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaaa                                   >  1:3835720/1‑93 (MQ=255)

                                                        acaTATTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaaaa                                  >  1:3077774/1‑93 (MQ=255)

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

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

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

                                                           tatTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACGGAGATGCCCGCAAAATcc                               <  1:3060955/93‑1 (MQ=255)

                                                           tatTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATcc                               <  2:114844/93‑1 (MQ=255)

                                                            atTCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCg                              >  1:1241404/1‑93 (MQ=255)

                                                              tCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGa                            <  2:5585670/93‑1 (MQ=255)

                                                              tCAAAGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGa                            <  2:1796641/93‑1 (MQ=255)

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

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

                                                                aaaGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATg                          <  2:1054218/93‑1 (MQ=255)

                                                                aaaGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATg                          <  2:4086910/93‑1 (MQ=255)

                                                                 aaGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATgg                         <  1:1140271/93‑1 (MQ=255)

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

                                                                  aGCCGTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGa                        <  1:858704/93‑1 (MQ=255)

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

                                                                      gTTAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTc                    >  2:5709101/1‑93 (MQ=255)

                                                                       ttAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCg                   >  2:42274/1‑93 (MQ=255)

                                                                        tAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGa                  >  2:786325/1‑93 (MQ=255)

                                                                        tAAGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGa                  <  2:3315747/93‑1 (MQ=255)

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

                                                                          aGGATCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAAc                <  2:2941844/93‑1 (MQ=255)

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

                                                                             aTCAGCAGCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCa             <  2:1519784/93‑1 (MQ=255)

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

                                                                               cagcagCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgc           <  1:2431591/93‑1 (MQ=255)

                                                                               cagcagCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgc           >  1:4957601/1‑93 (MQ=255)

                                                                               cagcagCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgc           <  1:1056555/93‑1 (MQ=255)

                                                                               cagcagCGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAga           >  2:5144977/1‑92 (MQ=255)

                                                                                 gcagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgcgc         >  2:5092780/1‑93 (MQ=255)

                                                                                 gcagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgcgc         >  2:5297340/1‑93 (MQ=255)

                                                                                 gcagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTAAAATGCTCAACAGAGATGCCCGCAACAGCCGGATGGAAGTCGAACGCAgcgc         <  1:450525/93‑1 (MQ=255)

                                                                                  cagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCt        <  2:3862413/93‑1 (MQ=255)

                                                                                  cagcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCt        >  1:6281196/1‑93 (MQ=255)

                                                                                    gcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGt      >  1:2337239/1‑93 (MQ=255)

                                                                                    gcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGt      >  2:1254478/1‑93 (MQ=255)

                                                                                    gcGGCATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGt      >  1:5522354/1‑93 (MQ=255)

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

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

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

                                                                                        cATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTcc  >  1:4390538/1‑93 (MQ=255)

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

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

                                                                                        cATTTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTcc  <  2:1606009/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` `< 24 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`