BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	981,066	A→T	D266V (GAC→GTC)	ycbB →	putative amidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	981,066	0	A	T	100.0%	834.4 / NA	221	D266V (GAC→GTC)	ycbB	putative amidase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base T (109/112); total (109/112)

Only 100 of 228 total aligned reads displayed.

AACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAA  >  NC_000913/980976‑981155

aaCTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCg                                                                                         >  1:1244618/1‑93 (MQ=255)

aaCTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCg                                                                                         >  2:2738007/1‑93 (MQ=255)

    gACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCggggg                                                                                     >  1:4460073/1‑93 (MQ=255)

    gACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCggggg                                                                                     <  2:222929/93‑1 (MQ=255)

    gACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCggggg                                                                                     <  2:3848804/93‑1 (MQ=255)

          cAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGaa                                                                               >  1:4331379/1‑93 (MQ=255)

          cAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGaa                                                                               >  2:1354156/1‑93 (MQ=255)

           aaaGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGcaa                                                                              >  1:5600545/1‑93 (MQ=255)

            aaGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCgcaca                                                                             >  2:2578354/1‑89 (MQ=255)

            aaGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGaaaa                                                                             <  2:1570001/93‑1 (MQ=255)

             aGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAAt                                                                            >  1:190735/1‑93 (MQ=255)

             aGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAAt                                                                            <  1:5461091/93‑1 (MQ=255)

               cAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTa                                                                          <  2:3823379/93‑1 (MQ=255)

                aaCGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTAc                                                                         <  1:4414657/93‑1 (MQ=255)

                 aCGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTAct                                                                        >  2:2401749/1‑93 (MQ=255)

                 aCGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTAct                                                                        >  2:6429121/1‑93 (MQ=255)

                  cGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTActc                                                                       <  2:5168706/93‑1 (MQ=255)

                  cGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTActc                                                                       <  1:4905291/93‑1 (MQ=255)

                   gTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTActct                                                                      >  2:1406129/1‑93 (MQ=255)

                    ttGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTa                                                                     >  1:2979734/1‑93 (MQ=255)

                    ttGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTa                                                                     >  2:3707787/1‑93 (MQ=255)

                     tGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTAc                                                                    >  1:338720/1‑93 (MQ=255)

                      gcgcCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTAcc                                                                   <  2:4690401/93‑1 (MQ=255)

                      gcgcCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTAcc                                                                   <  1:4880861/93‑1 (MQ=255)

                      gcgcCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTAcc                                                                   <  2:568910/93‑1 (MQ=255)

                       cgcCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCt                                                                  <  1:4788851/93‑1 (MQ=255)

                       cgcCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCt                                                                  <  2:1011298/93‑1 (MQ=255)

                           cAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCg                                                              >  2:1069061/1‑93 (MQ=255)

                           cAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCg                                                              >  2:54279/1‑93 (MQ=255)

                           cAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCg                                                              >  1:4245493/1‑93 (MQ=255)

                            aGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGa                                                             >  1:1800497/1‑93 (MQ=255)

                             gggCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGAt                                                            >  2:2280498/1‑93 (MQ=255)

                             gggCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGAt                                                            >  1:5342710/1‑93 (MQ=255)

                              ggCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATg                                                           <  1:4603354/93‑1 (MQ=255)

                              ggCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATg                                                           >  1:3154599/1‑93 (MQ=255)

                              ggCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATg                                                           <  2:4296496/93‑1 (MQ=255)

                              ggCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATg                                                           >  2:1509392/1‑93 (MQ=255)

                               gCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGa                                                          >  1:3716525/1‑93 (MQ=255)

                                cAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGac                                                         <  1:3757560/93‑1 (MQ=255)

                                 aGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGaca                                                        <  2:2784/93‑1 (MQ=255)

                                 aGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGaca                                                        <  2:3583408/93‑1 (MQ=255)

                                  gTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGacac                                                       >  1:3032865/1‑93 (MQ=255)

                                  gTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGacac                                                       >  2:1526800/1‑93 (MQ=255)

                                    ggAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGc                                                     <  1:1405147/93‑1 (MQ=255)

                                     gAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGcc                                                    <  1:4977556/93‑1 (MQ=255)

                                       gTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCaa                                                  <  2:2092478/93‑1 (MQ=255)

                                       gTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCaa                                                  <  1:2088933/93‑1 (MQ=255)

                                       gTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCaa                                                  >  2:5951523/1‑93 (MQ=255)

                                        tAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAAc                                                 >  2:7555/1‑93 (MQ=255)

                                          acgacgTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTg                                               <  2:1447306/93‑1 (MQ=255)

                                            gacgTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGAc                                             <  1:5291516/93‑1 (MQ=255)

                                            gacgTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGAc                                             >  2:5857906/1‑93 (MQ=255)

                                            gacgTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGAc                                             >  2:4177170/1‑93 (MQ=255)

                                             acgTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGAcg                                            >  1:3469374/1‑93 (MQ=255)

                                               gTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGAcgcg                                          >  2:600329/1‑93 (MQ=255)

                                                tACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCgg                                         >  2:5592282/1‑93 (MQ=255)

                                                 aCCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGt                                        <  1:1009273/93‑1 (MQ=255)

                                                 aCCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGt                                        >  1:5689229/1‑93 (MQ=255)

                                                  ccGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTa                                       <  2:4912279/93‑1 (MQ=255)

                                                  ccGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTa                                       >  1:6436117/1‑93 (MQ=255)

                                                  ccGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTa                                       <  1:2420049/93‑1 (MQ=255)

                                                   cGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAg                                      >  2:2982851/1‑93 (MQ=255)

                                                    ggCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGt                                     >  2:2323001/1‑93 (MQ=255)

                                                     gCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTc                                    >  1:6156373/1‑93 (MQ=255)

                                                     gCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTc                                    <  1:4387678/93‑1 (MQ=255)

                                                      cGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCa                                   >  1:2641585/1‑93 (MQ=255)

                                                       gTTGCGCGAAATATTGGAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAg                                  >  1:2902154/1‑93 (MQ=255)

                                                       gTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAg                                  >  2:6457176/1‑93 (MQ=255)

                                                       gTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAg                                  <  2:1363973/93‑1 (MQ=255)

                                                         tGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGcc                                >  2:4526866/1‑93 (MQ=255)

                                                            gcgAAATATTGCAACGCACAGGCATGTTGGTTGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCAt                             <  2:1091919/93‑1 (MQ=255)

                                                            gcgAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCAt                             >  1:3619295/1‑93 (MQ=255)

                                                             cgAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATc                            <  2:5052277/93‑1 (MQ=255)

                                                             cgAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATc                            <  2:566227/93‑1 (MQ=255)

                                                             cgAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATc                            <  2:2840529/93‑1 (MQ=255)

                                                              gAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATcc                           <  2:3913663/93‑1 (MQ=255)

                                                              gAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATcc                           <  2:3127251/93‑1 (MQ=255)

                                                               aaaTATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCg                          >  1:1492016/1‑93 (MQ=255)

                                                               aaaTATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCg                          >  2:1651916/1‑93 (MQ=255)

                                                                 atatTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCt                        <  1:2480475/93‑1 (MQ=255)

                                                                 atatTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCt                        >  2:5015551/1‑93 (MQ=255)

                                                                  tatTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCtt                       >  1:1139689/1‑92 (MQ=255)

                                                                   atTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGt                      <  1:4066686/93‑1 (MQ=255)

                                                                    ttGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGtt                     <  2:3631480/93‑1 (MQ=255)

                                                                      gCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTAc                   <  2:4513390/93‑1 (MQ=255)

                                                                         aCGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGt                <  1:5784992/93‑1 (MQ=255)

                                                                         aCGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCGGTTACTGt                <  1:6301365/93‑1 (MQ=255)

                                                                           gCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTg              >  1:3417485/1‑93 (MQ=255)

                                                                           gCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTg              <  2:6279104/93‑1 (MQ=255)

                                                                            cacaGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGa             >  1:5970/1‑93 (MQ=255)

                                                                             acaGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGaa            >  1:750453/1‑93 (MQ=255)

                                                                               aGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAc          >  2:3160121/1‑93 (MQ=255)

                                                                                ggCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAca         <  2:5769358/93‑1 (MQ=255)

                                                                                ggCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAca         >  1:1356699/1‑93 (MQ=255)

                                                                                  cATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAcagc       >  1:1570459/1‑93 (MQ=255)

                                                                                   aTGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAcagca      >  1:2265171/1‑93 (MQ=255)

                                                                                   aTGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAcagca      >  1:3322520/1‑93 (MQ=255)

                                                                                   aTGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAcagca      <  2:1819361/93‑1 (MQ=255)

                                                                                     gTTGGTCGCGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGa    <  2:4067309/93‑1 (MQ=255)

                                                                                      ttGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGaa   >  2:2164357/1‑93 (MQ=255)

AACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAA  >  NC_000913/980976‑981155

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