BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,147,941	A→G	V356A (GTC→GCC)	b2070 ←	putative chaperonin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,147,941	0	A	G	99.4%	1067.6 / ‑3.0	316	V356A (GTC→GCC)	b2070	putative chaperonin
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (138/176); total (139/177)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-01

Only 100 of 338 total aligned reads displayed.

AAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACAC  >  NC_000913/2147851‑2148030

aaaGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCt                                                                                         >  2:7025259/1‑93 (MQ=255)

 aaGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCtt                                                                                        <  2:6174113/93‑1 (MQ=255)

 aaGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCtt                                                                                        <  1:5334704/93‑1 (MQ=255)

  aGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCttt                                                                                       >  1:4113071/1‑93 (MQ=255)

  aGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCttt                                                                                       <  1:1355002/93‑1 (MQ=255)

   ggAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTc                                                                                      <  1:8677972/93‑1 (MQ=255)

     aaGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGc                                                                                    >  2:6367902/1‑93 (MQ=255)

     aaGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGc                                                                                    >  1:77625/1‑93 (MQ=255)

      aGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGcc                                                                                   >  1:7344159/1‑93 (MQ=255)

          gCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGca                                                                               <  1:2899731/93‑1 (MQ=255)

             tGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACaa                                                                            >  1:6337054/1‑93 (MQ=255)

              gCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACaaa                                                                           >  2:5783717/1‑93 (MQ=255)

               ccAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACaaaa                                                                          <  2:7763666/93‑1 (MQ=255)

               ccAGTCGGTTGGAAATGATCGGGTAAGGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCAGCGCCCAGGCTTTCGCCAGCACaaaa                                                                          <  1:3902131/93‑1 (MQ=255)

                cAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAAt                                                                         <  2:2201185/93‑1 (MQ=255)

                cAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAAt                                                                         >  2:7094551/1‑93 (MQ=255)

                  gTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTc                                                                       >  1:7902426/1‑93 (MQ=255)

                   tCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTcc                                                                      <  2:173739/93‑1 (MQ=255)

                   gctgTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATGcc                                                                      <  2:1017931/90‑1 (MQ=255)

                      gTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCccgc                                                                   <  1:4279984/93‑1 (MQ=255)

                      gTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCccgc                                                                   >  2:8166332/1‑93 (MQ=255)

                       ttGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCccgcc                                                                  >  2:343936/1‑93 (MQ=255)

                         tgAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGTTTCTACTGACGCCATCCCCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCgg                                                                <  2:6489227/92‑1 (MQ=255)

                         ggAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCgg                                                                <  1:4969310/93‑1 (MQ=255)

                         ggAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCgg                                                                <  1:831791/93‑1 (MQ=255)

                         ggAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCgg                                                                >  2:5558225/1‑93 (MQ=255)

                          gAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCggt                                                               <  2:5159281/93‑1 (MQ=255)

                          gAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCggg                                                               >  2:1559876/1‑92 (MQ=255)

                                aTCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAAt                                                         <  2:6632704/93‑1 (MQ=255)

                                aTCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAAt                                                         >  2:4529106/1‑93 (MQ=255)

                                 tCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATg                                                        <  1:1056538/93‑1 (MQ=255)

                                  cGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGa                                                       <  1:3748640/93‑1 (MQ=255)

                                   gggTAATGGCATTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGAc                                                      <  1:6234729/93‑1 (MQ=255)

                                   gggTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGAc                                                      <  2:5628103/93‑1 (MQ=255)

                                   gggTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGAc                                                      >  1:2816131/1‑93 (MQ=255)

                                    ggTAATGGGAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGAcc                                                     <  1:1495619/93‑1 (MQ=255)

                                    ggTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGAcc                                                     >  2:1202943/1‑93 (MQ=255)

                                    ggTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGAcc                                                     <  1:9057/93‑1 (MQ=255)

                                    ggTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGAcc                                                     >  2:5928350/1‑93 (MQ=255)

                                     gTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACct                                                    <  1:5842947/93‑1 (MQ=255)

                                      tAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACctc                                                   <  1:7614870/93‑1 (MQ=255)

                                      tAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACctc                                                   >  2:3306973/1‑93 (MQ=255)

                                       aaTGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACctct                                                  >  2:8448160/1‑93 (MQ=255)

                                       aaTGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACctct                                                  <  1:3270888/93‑1 (MQ=255)

                                         tGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTc                                                <  2:2082799/93‑1 (MQ=255)

                                              gTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCga                                           >  2:2502361/1‑93 (MQ=255)

                                              gTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCga                                           >  1:3556965/1‑93 (MQ=255)

                                              gTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCga                                           <  1:1893881/93‑1 (MQ=255)

                                              gTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCga                                           >  2:3993519/1‑93 (MQ=255)

                                                tGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCgaga                                         <  1:302948/93‑1 (MQ=255)

                                                tGCAATGCTGGCGGGGCAACTGACCCCATCCCCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCgaga                                         <  2:6718284/93‑1 (MQ=255)

                                                 gCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTGCCGCCGGTGGTAATGACCTCTTCGCCgagac                                        >  1:5546111/1‑92 (MQ=255)

                                                 gCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGAt                                        <  2:3820544/93‑1 (MQ=255)

                                                 gCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGAt                                        <  1:5251840/93‑1 (MQ=255)

                                                  cAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATc                                       <  2:876746/93‑1 (MQ=255)

                                                   aaTCCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCa                                      <  1:2581114/93‑1 (MQ=255)

                                                     tGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAAt                                    >  2:4952635/1‑93 (MQ=255)

                                                      gCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATa                                   >  2:7467582/1‑93 (MQ=255)

                                                      gCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATa                                   <  2:7328116/93‑1 (MQ=255)

                                                        tGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATaaa                                 >  1:4055677/1‑93 (MQ=255)

                                                         ggCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGTCCTCTTCGCCGAGATCAATaaaa                                >  2:4751621/1‑93 (MQ=255)

                                                         ggCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATaaaa                                <  1:5958238/93‑1 (MQ=255)

                                                            ggTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAg                             >  1:8236745/1‑93 (MQ=255)

                                                            ggTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAg                             <  1:3402564/93‑1 (MQ=255)

                                                             gTGCTACTGACGCCATCCCCTTCGCCAGGGCTTGCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGa                            <  1:6179386/93‑1 (MQ=255)

                                                             gTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGa                            <  2:8296101/93‑1 (MQ=255)

                                                             gTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGa                            >  1:1963260/1‑93 (MQ=255)

                                                             gTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGa                            >  1:3966495/1‑93 (MQ=255)

                                                               gCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATa                          >  1:2250751/1‑93 (MQ=255)

                                                               gCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATa                          >  2:2726057/1‑93 (MQ=255)

                                                                cTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATaa                         <  2:4312875/93‑1 (MQ=255)

                                                                 tACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAac                        >  1:2146164/1‑93 (MQ=255)

                                                                 tACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAac                        <  2:268983/93‑1 (MQ=255)

                                                                  aCTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAaca                       <  1:1514283/93‑1 (MQ=255)

                                                                   cTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAacac                      <  2:3704855/93‑1 (MQ=255)

                                                                   cTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGCGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAacac                      >  2:3890249/1‑93 (MQ=255)

                                                                     gACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCt                    >  2:3052253/1‑93 (MQ=255)

                                                                       cGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTg                  <  2:163903/93‑1 (MQ=255)

                                                                         ccATCCCCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGAt                <  2:5129681/93‑1 (MQ=255)

                                                                          catcgccatcgcCCAGGCTTCCGCAAGCATAAAATTCCCGCCGGTGTTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGAtt               <  1:5159358/93‑1 (MQ=255)

                                                                            tCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTAt             <  2:2825040/93‑1 (MQ=255)

                                                                             ccccATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATc            <  1:3891635/91‑1 (MQ=255)

                                                                             cACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGCTTATc            >  1:7291535/1‑93 (MQ=255)

                                                                             cACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATc            <  1:1689428/93‑1 (MQ=255)

                                                                              aCCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATcc           >  1:1108898/1‑93 (MQ=255)

                                                                               ccATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCa          <  2:7940658/93‑1 (MQ=255)

                                                                                cATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAg         <  1:68325/93‑1 (MQ=255)

                                                                                 aTCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAgc        <  1:4321120/93‑1 (MQ=255)

                                                                                 aTCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAgc        <  1:8068515/93‑1 (MQ=255)

                                                                                   cGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAgcgc      <  2:2590944/93‑1 (MQ=255)

                                                                                    gCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGca     <  1:5729226/93‑1 (MQ=255)

                                                                                    gCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGca     >  2:3571947/1‑93 (MQ=255)

                                                                                    gCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGca     <  2:5342281/93‑1 (MQ=255)

                                                                                     cccAGGCTTTCGCCAGCACAAAATTCCCGCCTGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGcac    >  2:2389034/1‑93 (MQ=255)

                                                                                      ccAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGcaca   <  2:2760619/93‑1 (MQ=255)

                                                                                      ccAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGcaca   >  2:1906373/1‑93 (MQ=255)

                                                                                       cAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCacac  >  1:23746/1‑93 (MQ=255)

                                                                                       cAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCacac  >  1:695495/1‑93 (MQ=255)

                                                                                       cAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCacac  >  2:4353586/1‑93 (MQ=255)

                                                                                       cAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCacac  >  2:5948541/1‑93 (MQ=255)

AAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACAC  >  NC_000913/2147851‑2148030

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