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	84.4%	793.0 / 129.1	262	D266V (GAC→GTC)	ycbB	putative amidase
Reads supporting (aligned to +/- strand): ref base A (23/18); new base T (84/137); total (107/155)
Fisher's exact test for biased strand distribution p-value = 3.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

Only 100 of 266 total aligned reads displayed.

CCAACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAACT  >  NC_000913/980974‑981157

  aaCTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCg                                                                                           <  2:5089162/93‑1 (MQ=255)

  aaCTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCg                                                                                           <  1:4286080/93‑1 (MQ=255)

  aaCTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCg                                                                                           <  1:2904369/93‑1 (MQ=255)

  aaCTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCg                                                                                           <  2:6997647/93‑1 (MQ=255)

   aCTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCgg                                                                                          >  1:3950664/1‑93 (MQ=255)

    cTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCggg                                                                                         <  1:5723004/93‑1 (MQ=255)

       aCCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGACgggggg                                                                                      >  2:2819761/1‑93 (MQ=255)

             aaaGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGaaa                                                                                >  2:3654033/1‑93 (MQ=255)

             aaaGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGaaa                                                                                >  2:6737098/1‑93 (MQ=255)

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

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

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

                  aaCGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTAc                                                                           >  1:853705/1‑93 (MQ=255)

                   aCGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTAct                                                                          >  2:6886017/1‑93 (MQ=255)

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

                    cGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTActc                                                                         >  2:7418470/1‑93 (MQ=255)

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

                     gTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTActct                                                                        <  2:5969597/93‑1 (MQ=255)

                     gTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTActct                                                                        <  1:3336738/93‑1 (MQ=255)

                     gTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTActct                                                                        <  1:5861729/93‑1 (MQ=255)

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

                       tGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTAc                                                                      <  2:5350639/93‑1 (MQ=255)

                         cgcCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTACTCTACCt                                                                    <  1:3835935/93‑1 (MQ=255)

                            ccAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGc                                                                 >  2:469793/1‑93 (MQ=255)

                             cAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCg                                                                <  1:5608543/93‑1 (MQ=255)

                              aGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGa                                                               <  1:4886845/93‑1 (MQ=255)

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

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

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

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

                                    gTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGacac                                                         <  1:2405876/93‑1 (MQ=255)

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

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

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

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

                                           aaCGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACt                                                  >  2:3828647/1‑93 (MQ=255)

                                           aaCGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACt                                                  >  2:4516612/1‑93 (MQ=255)

                                            acgacgTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTg                                                 >  1:4361993/1‑93 (MQ=255)

                                             cgacgTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCAAACTGa                                                <  2:7094153/93‑1 (MQ=255)

                                              gacgTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGAc                                               >  1:1804448/1‑93 (MQ=255)

                                              gacgTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGAc                                               <  1:1899651/93‑1 (MQ=255)

                                               acgTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGAcg                                              <  1:6526036/93‑1 (MQ=255)

                                                cgTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGAcgc                                             <  1:3523852/93‑1 (MQ=255)

                                                cgTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGAcgc                                             <  2:6219326/93‑1 (MQ=255)

                                                 gTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGAcgcg                                            <  2:6353517/93‑1 (MQ=255)

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

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

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

                                                     cGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAg                                        <  1:7108773/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                          ttGCGCGAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGc                                   >  1:231155/1‑93 (MQ=255)

                                                             cgcgAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCa                                >  1:1717971/1‑93 (MQ=255)

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

                                                              gcgAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCAt                               <  1:5308023/93‑1 (MQ=255)

                                                              gcgAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCAt                               <  2:6674383/93‑1 (MQ=255)

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

                                                               cgAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATc                              <  2:7468057/93‑1 (MQ=255)

                                                               cgAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATc                              <  2:2229134/93‑1 (MQ=255)

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

                                                                gAAATATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATcc                             >  2:4363201/1‑93 (MQ=255)

                                                                  aaTATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGc                           >  1:3149405/1‑93 (MQ=255)

                                                                  aaTATTGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGc                           >  1:1625658/1‑93 (MQ=255)

                                                                    tatTGCAACGCACAGGCATGTTGGTCGGGGTGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTg                         <  1:1288613/93‑1 (MQ=255)

                                                                      ttGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGtt                       >  1:3076448/1‑93 (MQ=255)

                                                                       tGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTa                      <  2:1367864/93‑1 (MQ=255)

                                                                       tGCAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTa                      <  1:1599594/93‑1 (MQ=255)

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

                                                                         cAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACt                    <  1:6877628/93‑1 (MQ=255)

                                                                         cAACGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACt                    <  2:6212971/93‑1 (MQ=255)

                                                                         cAACGCACAGGCATGTTGATCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACt                    <  2:3220904/93‑1 (MQ=255)

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

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

                                                                            cGCACAGGCATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGtt                 <  1:275085/93‑1 (MQ=255)

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

                                                                              cacaGGCATGTTGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGa               <  1:6382336/93‑1 (MQ=255)

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

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

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

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

                                                                                    cATGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAcagc         <  1:2568177/93‑1 (MQ=255)

                                                                                    cATGTTGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAcagc         >  1:3403118/1‑93 (MQ=255)

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

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

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

                                                                                      tGTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAcagcag       >  1:7059976/1‑93 (MQ=255)

                                                                                       gTTGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGa      >  1:3202993/1‑93 (MQ=255)

                                                                                         tGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGaaa    <  1:589564/93‑1 (MQ=255)

                                                                                         tGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGaaa    <  1:6053016/93‑1 (MQ=255)

                                                                                         tGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGaaa    <  1:3156324/93‑1 (MQ=255)

                                                                                         tGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGaaa    <  1:6161744/93‑1 (MQ=255)

                                                                                         tGGTCGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGaaa    <  1:6689320/93‑1 (MQ=255)

                                                                                         tGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGaaa    <  2:5212994/93‑1 (MQ=255)

                                                                                           gACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAACt  <  1:5448439/93‑1 (MQ=255)

CCAACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAACT  >  NC_000913/980974‑981157

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