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%	957.3 / ‑3.8	257	coding (2014/2028 nt)	b2974	putative endoglucanase
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (131/125); total (132/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

Only 100 of 271 total aligned reads displayed.

TTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TT‑CTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTCC  >  NC_000913/3115024‑3115202

 ttGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TT‑cc                                                                                       <  2:6337275/93‑2 (MQ=255)

  tGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TT‑cct                                                                                      <  1:5388371/93‑3 (MQ=255)

   gCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TT‑cctg                                                                                     <  1:4170282/93‑4 (MQ=255)

      ggTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCgg                                                                                   >  2:3813415/1‑93 (MQ=255)

        tATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGAt                                                                                 <  1:5846327/93‑1 (MQ=255)

         aTCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATc                                                                                <  1:4973553/93‑1 (MQ=255)

         aTCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATc                                                                                <  1:1703832/93‑1 (MQ=255)

          tCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCg                                                                               <  1:5138637/93‑1 (MQ=255)

          tCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCg                                                                               <  1:6075119/93‑1 (MQ=255)

          tCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCg                                                                               >  2:1372880/1‑93 (MQ=255)

           cTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCga                                                                              <  1:1523576/93‑1 (MQ=255)

           cTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCga                                                                              <  1:2456962/93‑1 (MQ=255)

            tGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCgag                                                                             >  1:5363728/1‑93 (MQ=255)

            tGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGCCACATATTCAAAGCCGTTAAGGATCACCAGCGGCAT‑TTCCTGCGGATCgcg                                                                             >  1:6353817/1‑91 (MQ=255)

             gCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCgaga                                                                            >  2:2667461/1‑93 (MQ=255)

             gCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCgaga                                                                            <  2:4508617/93‑1 (MQ=255)

              cACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGAt                                                                           <  1:1232819/93‑1 (MQ=255)

              cACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGAt                                                                           <  2:5551323/93‑1 (MQ=255)

              cACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGAt                                                                           <  1:852468/93‑1 (MQ=255)

              cACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGAt                                                                           >  2:4462757/1‑93 (MQ=255)

               aCGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATc                                                                          >  1:4225357/1‑93 (MQ=255)

                 gCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGc                                                                        >  1:5768853/1‑93 (MQ=255)

                    gCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCAt                                                                     >  2:882036/1‑93 (MQ=255)

                    gCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCAt                                                                     >  1:2807558/1‑93 (MQ=255)

                     ggggATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATa                                                                    <  2:2528599/92‑1 (MQ=255)

                     cGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATa                                                                    >  2:5750427/1‑93 (MQ=255)

                       ggATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGc                                                                  >  2:2204235/1‑93 (MQ=255)

                        gATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCt                                                                 >  2:1621031/1‑93 (MQ=255)

                         aTTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTa                                                                <  2:4638137/93‑1 (MQ=255)

                          ttGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTat                                                               >  2:5350910/1‑92 (MQ=255)

                           tGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACt                                                              >  1:4390215/1‑93 (MQ=255)

                              aTAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTaa                                                           >  2:393246/1‑93 (MQ=255)

                              aTAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTaa                                                           >  2:4529450/1‑93 (MQ=255)

                              aTAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTaa                                                           >  2:3508188/1‑93 (MQ=255)

                               tAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTaaa                                                          >  2:4139551/1‑93 (MQ=255)

                               tAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTaaa                                                          >  1:1941067/1‑93 (MQ=255)

                               tAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTaaa                                                          >  2:6472628/1‑93 (MQ=255)

                               tAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTaaa                                                          >  2:3311450/1‑93 (MQ=255)

                               tAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTaaa                                                          >  1:3269977/1‑93 (MQ=255)

                                aaGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAAt                                                         <  2:2942726/93‑1 (MQ=255)

                                 aGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATg                                                        >  1:229428/1‑93 (MQ=255)

                                  ggATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGc                                                       <  1:6554137/93‑1 (MQ=255)

                                   gATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCt                                                      >  2:4064570/1‑93 (MQ=255)

                                    aTCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTc                                                     <  1:3315865/93‑1 (MQ=255)

                                     acGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCa                                                    <  2:760417/92‑1 (MQ=255)

                                      cGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCaa                                                   >  1:5617687/1‑93 (MQ=255)

                                      cGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCaa                                                   <  1:2835368/93‑1 (MQ=255)

                                        ttACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACa                                                 >  1:3799420/1‑93 (MQ=255)

                                        ttACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACa                                                 >  2:730077/1‑93 (MQ=255)

                                         tACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACag                                                <  2:4762757/93‑1 (MQ=255)

                                         tACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACag                                                <  1:1104162/93‑1 (MQ=255)

                                         tACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACag                                                <  1:5318126/93‑1 (MQ=255)

                                         tACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGAAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACag                                                >  2:2767700/1‑93 (MQ=255)

                                          aCCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACaga                                               >  1:1434559/1‑93 (MQ=255)

                                          aCCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACaga                                               <  2:1428426/93‑1 (MQ=255)

                                          aCCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACaga                                               >  2:546305/1‑93 (MQ=255)

                                              cccTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATg                                           <  1:5963216/92‑1 (MQ=255)

                                              aCCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATg                                           <  1:6180107/93‑1 (MQ=255)

                                                cTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGcc                                         <  2:131331/93‑1 (MQ=255)

                                                 tGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGccc                                        <  2:6607570/93‑1 (MQ=255)

                                                    gTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCa                                     <  1:3627115/93‑1 (MQ=255)

                                                    gTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCa                                     >  2:4178644/1‑93 (MQ=255)

                                                     tCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaa                                    <  1:496679/93‑1 (MQ=255)

                                                     tCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaa                                    <  1:5442336/93‑1 (MQ=255)

                                                       acaTATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaaaa                                  <  2:6059643/93‑1 (MQ=255)

                                                       acaTATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaaaa                                  <  1:2949919/93‑1 (MQ=255)

                                                       acaTATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaaaa                                  >  2:678557/1‑93 (MQ=255)

                                                       acaTATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaaaa                                  >  1:189357/1‑93 (MQ=255)

                                                        caTATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAAt                                 >  2:6277601/1‑93 (MQ=255)

                                                        caTATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAAt                                 <  2:5162987/93‑1 (MQ=255)

                                                         atatTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATc                                >  1:4748432/1‑93 (MQ=255)

                                                          tatTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATcc                               <  2:3371130/93‑1 (MQ=255)

                                                          tatTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATcc                               <  2:4725407/93‑1 (MQ=255)

                                                            ttCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCgg                             >  2:2092089/1‑93 (MQ=255)

                                                               aaaGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATg                          >  1:705568/1‑93 (MQ=255)

                                                                 aGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGa                        >  1:3934882/1‑93 (MQ=255)

                                                                 aGCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGa                        <  2:3080531/93‑1 (MQ=255)

                                                                  gCCGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGaa                       <  1:1064446/93‑1 (MQ=255)

                                                                   ccGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAg                      <  2:1525045/93‑1 (MQ=255)

                                                                   ccGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAg                      <  2:1245261/93‑1 (MQ=255)

                                                                    cGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGt                     <  2:6303375/93‑1 (MQ=255)

                                                                    cGTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGt                     <  1:2634750/93‑1 (MQ=255)

                                                                     gTTAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTc                    <  1:2261577/93‑1 (MQ=255)

                                                                      ttAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCg                   >  1:2520267/1‑93 (MQ=255)

                                                                      ttAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCg                   >  1:4722812/1‑93 (MQ=255)

                                                                       tgaGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGa                  >  1:5278293/3‑93 (MQ=255)

                                                                       tAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGa                  >  1:4601523/1‑93 (MQ=255)

                                                                       tAAGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGa                  >  2:4029175/1‑93 (MQ=255)

                                                                         aGGATCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAAc                <  2:2589237/93‑1 (MQ=255)

                                                                            aTCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCa             <  2:1170520/93‑1 (MQ=255)

                                                                             tCAGCAGCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAg            <  2:2410197/93‑1 (MQ=255)

                                                                              cagcagCGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgc           <  1:3121413/93‑1 (MQ=255)

                                                                                gcagcGGCATGTC‑CTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgcgc         <  2:320073/93‑1 (MQ=255)

                                                                                  agcGGCAT‑TTCCTGCGGATCGAGGTCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTg       >  1:764201/1‑93 (MQ=255)

                                                                                   gcGGCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGt      >  2:501807/1‑93 (MQ=255)

                                                                                     ggCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGttt    >  1:3527396/1‑93 (MQ=255)

                                                                                     ggCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGttt    <  1:5691762/93‑1 (MQ=255)

                                                                                      gCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTc   >  2:4321586/1‑93 (MQ=255)

                                                                                      gCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTc   >  2:1912688/1‑93 (MQ=255)

                                                                                      gCAT‑TTCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTc   >  2:1066680/1‑93 (MQ=255)

TTTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCAT‑TT‑CTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTCC  >  NC_000913/3115024‑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: `‑`