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.2%	825.4 / ‑1.5	254	coding (2014/2028 nt)	b2974	putative endoglucanase
Reads supporting (aligned to +/- strand): ref base . (1/1); new base C (132/120); total (133/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 = 7.74e-01

Only 100 of 270 total aligned reads displayed.

TTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑T‑CTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTCC  >  NC_000913/3115025‑3115202

   cTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑T‑cctgc                                                                                    <  1:2340552/93‑5 (MQ=255)

     ggTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCgg                                                                                   >  2:2876291/1‑93 (MQ=255)

      gTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGa                                                                                  <  2:2139715/93‑1 (MQ=255)

      gTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGa                                                                                  >  1:5030239/1‑93 (MQ=255)

       tATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGAt                                                                                 >  2:2451696/1‑93 (MQ=255)

        aTCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATc                                                                                <  2:6020040/93‑1 (MQ=255)

          cTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCga                                                                              <  1:4018374/93‑1 (MQ=255)

          cTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCga                                                                              <  1:1724065/93‑1 (MQ=255)

          cTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCga                                                                              <  2:2030022/93‑1 (MQ=255)

           tGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCgag                                                                             >  2:3606315/1‑93 (MQ=255)

             cACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGAt                                                                           <  1:6116084/93‑1 (MQ=255)

                gCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGc                                                                        >  1:2321184/1‑93 (MQ=255)

                 cAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGcc                                                                       <  2:1769308/93‑1 (MQ=255)

                   gCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCAt                                                                     >  2:2766397/1‑93 (MQ=255)

                   gCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCAt                                                                     >  1:95750/1‑93 (MQ=255)

                   gCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCAt                                                                     >  1:1143202/1‑93 (MQ=255)

                    cGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATa                                                                    <  2:5924182/93‑1 (MQ=255)

                    cGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATa                                                                    <  1:4306317/93‑1 (MQ=255)

                    agggATTGCATAAGGATCGTTACCGACGTGAGTCTCATATTCACAGCCGTTAAGGATCAGCAGCTTCATT‑TCCTGCGGATCGAGATCGCCATa                                                                    <  1:6215416/92‑1 (MQ=255)

                     gggATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAg                                                                   <  2:4653031/93‑1 (MQ=255)

                     gggATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAg                                                                   >  2:2926722/1‑93 (MQ=255)

                      ggATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAATGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGc                                                                  >  1:3375515/1‑93 (MQ=255)

                       gATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCt                                                                 >  2:3714469/1‑93 (MQ=255)

                         ttGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTAc                                                               >  2:5612456/1‑93 (MQ=255)

                         ttGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTAc                                                               >  1:3139620/1‑93 (MQ=255)

                         ttGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTAc                                                               <  1:2952845/93‑1 (MQ=255)

                         ttGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTAc                                                               >  2:1575949/1‑93 (MQ=255)

                            cATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTa                                                            >  2:4842192/1‑93 (MQ=255)

                            cATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTa                                                            >  2:1062581/1‑93 (MQ=255)

                              tAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTaaa                                                          <  1:5542167/93‑1 (MQ=255)

                              tAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTaaa                                                          >  1:4095397/1‑93 (MQ=255)

                              tAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTaaa                                                          >  2:5187434/1‑93 (MQ=255)

                                aGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATg                                                        >  2:5738011/1‑93 (MQ=255)

                                  gATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCt                                                      <  1:3757434/93‑1 (MQ=255)

                                   aTCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATtctc                                                     <  1:1649112/93‑1 (MQ=255)

                                   aTCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTc                                                     <  1:1522283/93‑1 (MQ=255)

                                   aTCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTc                                                     <  1:4841259/93‑1 (MQ=255)

                                   aTCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTc                                                     >  2:5428922/1‑93 (MQ=255)

                                    tCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCa                                                    >  2:407872/1‑93 (MQ=255)

                                     cGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCaa                                                   <  2:5484542/93‑1 (MQ=255)

                                      gTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAAc                                                  <  2:427114/93‑1 (MQ=255)

                                      gTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAAc                                                  >  1:3267034/1‑93 (MQ=255)

                                          cccACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGGTACTTAAATGCTCAACagag                                              <  2:2490055/93‑1 (MQ=255)

                                           cccccTGAGCCACATATTCAAAGCCGTTACGGATCAGCAGCGGAATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTCACTGCTCTACAgaga                                             >  2:463075/4‑93 (MQ=255)

                                           ccACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAgaga                                             >  2:1836821/1‑93 (MQ=255)

                                            cACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGAt                                            >  1:2103919/1‑93 (MQ=255)

                                            cACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGAt                                            <  2:129965/93‑1 (MQ=255)

                                             cccTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATg                                           <  2:2958479/92‑1 (MQ=255)

                                             aCCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATg                                           >  1:5239704/1‑93 (MQ=255)

                                              ccTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGc                                          >  1:5811968/1‑93 (MQ=255)

                                               cTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGcc                                         <  1:3834270/93‑1 (MQ=255)

                                                tGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGccc                                        >  1:3965555/1‑93 (MQ=255)

                                                tGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGccc                                        <  1:2902395/93‑1 (MQ=255)

                                                tGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGccc                                        <  2:2324156/93‑1 (MQ=255)

                                                   gTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCa                                     >  2:6220555/1‑93 (MQ=255)

                                                    tCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaa                                    >  2:751692/1‑93 (MQ=255)

                                                    tCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaa                                    >  1:3861238/1‑93 (MQ=255)

                                                      acaTATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaaaa                                  >  1:2611921/1‑93 (MQ=255)

                                                      acaTATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaaaa                                  >  1:1891425/1‑93 (MQ=255)

                                                      acaTATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCaaaa                                  >  2:4441393/1‑93 (MQ=255)

                                                       caTATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAAt                                 <  2:6045991/93‑1 (MQ=255)

                                                       caTATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAAt                                 >  1:796444/1‑93 (MQ=255)

                                                         tatTGAACGCCGTTAGGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATcc                               <  1:5764661/93‑1 (MQ=255)

                                                         tatTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATcc                               >  1:1672816/1‑93 (MQ=255)

                                                         tatTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATcc                               >  2:1035931/1‑93 (MQ=255)

                                                            tCAAAGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGa                            >  2:3264537/1‑93 (MQ=255)

                                                               aaGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATgg                         <  1:1376276/93‑1 (MQ=255)

                                                               aaGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATgg                         <  1:5512820/93‑1 (MQ=255)

                                                               aaGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATgg                         <  2:2283885/93‑1 (MQ=255)

                                                               aaGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATgg                         >  2:2603656/1‑93 (MQ=255)

                                                                aGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGa                        >  2:4375794/1‑93 (MQ=255)

                                                                aGCCGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGa                        <  1:3507562/93‑1 (MQ=255)

                                                                   cGTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGt                     <  1:4505972/93‑1 (MQ=255)

                                                                    gTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTc                    <  1:344063/93‑1 (MQ=255)

                                                                    gTTAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTc                    <  2:2219406/93‑1 (MQ=255)

                                                                     ttAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCg                   >  1:5782192/1‑93 (MQ=255)

                                                                     ttAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCg                   <  1:4175740/93‑1 (MQ=255)

                                                                     ttAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCg                   >  1:823823/1‑93 (MQ=255)

                                                                     ttAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCg                   >  2:380877/1‑93 (MQ=255)

                                                                      tAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGa                  <  1:4650796/93‑1 (MQ=255)

                                                                      tAAGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGa                  >  2:1350214/1‑93 (MQ=255)

                                                                       aaGGATCAGCAGCGGCATTGC‑CTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAAATGCCCGCAAAATCCGGATGGAAGTCGaa                 <  2:1432553/93‑1 (MQ=255)

                                                                       aaGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGaa                 <  1:5414042/93‑1 (MQ=255)

                                                                       aaGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGaa                 <  2:3671952/93‑1 (MQ=255)

                                                                       aaGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGaa                 <  1:3337303/93‑1 (MQ=255)

                                                                        aGGATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAAc                <  1:2652573/93‑1 (MQ=255)

                                                                         ggATCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACg               >  2:5026647/1‑93 (MQ=255)

                                                                           aTCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCa             >  2:436127/1‑93 (MQ=255)

                                                                           aTCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCa             <  2:120944/93‑1 (MQ=255)

                                                                           aTCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCa             <  1:4300468/93‑1 (MQ=255)

                                                                           aTCAGCAGCGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCa             <  1:2208130/93‑1 (MQ=255)

                                                                              agcagcGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAgcg          >  1:6033602/1‑93 (MQ=255)

                                                                              agcagaGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTAAACAGAGATGCCCGCAAAATCCGGATGGAAGCCGAACGCAgcg          >  1:4347534/1‑93 (MQ=255)

                                                                                  gcGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGt      >  1:5043589/1‑93 (MQ=255)

                                                                                  gcGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGt      >  1:1781796/1‑93 (MQ=255)

                                                                                  gcGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGt      <  1:2489682/93‑1 (MQ=255)

                                                                                  gcGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGt      <  1:610153/93‑1 (MQ=255)

                                                                                  gcGGCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGt      >  2:4102598/1‑93 (MQ=255)

                                                                                    ggCATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGttt    >  2:3135974/1‑93 (MQ=255)

                                                                                      cATT‑TCCTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTcc  <  1:5643511/93‑1 (MQ=255)

TTGCTGGTATCTGCACGCAGCGGGATTGCATAAGGATCGTTACCCACCTGAGTCACATATTCAAAGCCGTTAAGGATCAGCAGCGGCATT‑T‑CTGCGGATCGAGATCGCCATAGCTACTTAAATGCTCAACAGAGATGCCCGCAAAATCCGGATGGAAGTCGAACGCAGCGCTGTTTCC  >  NC_000913/3115025‑3115202

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