BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,240,123	T→C	M278V (ATG→GTG)	uxaA ←	altronate hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,240,123	0	T	C	83.6%	593.0 / 113.9	269	M278V (ATG→GTG)	uxaA	altronate hydrolase
Reads supporting (aligned to +/- strand): ref base T (16/28); new base C (117/108); total (133/136)
Fisher's exact test for biased strand distribution p-value = 6.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

Only 100 of 270 total aligned reads displayed.

GCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCATCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCTGGCTCGCGTTTATCGT  >  NC_000913/3240033‑3240215

gCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCg                                                                                            >  2:24104/1‑93 (MQ=255)

gCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCg                                                                                            >  1:1592226/1‑93 (MQ=255)

  aaCTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCATCggg                                                                                          >  1:6690158/1‑93 (MQ=255)

  aaCTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCggg                                                                                          <  2:4668311/93‑1 (MQ=255)

   aCTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGt                                                                                         >  1:1257139/1‑93 (MQ=255)

    cTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGtt                                                                                        >  2:1576062/1‑93 (MQ=255)

     tGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTa                                                                                       <  1:5846358/93‑1 (MQ=255)

      gCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAg                                                                                      >  1:4925900/1‑93 (MQ=255)

        tctGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCa                                                                                    <  1:1301065/93‑1 (MQ=255)

        tctGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCa                                                                                    >  2:2610449/1‑93 (MQ=255)

          tGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCATCGGGTTAGCAGt                                                                                  >  1:7066846/1‑93 (MQ=255)

            cgcCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTaa                                                                                <  2:1953143/93‑1 (MQ=255)

             gcCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAAt                                                                               <  1:2255042/93‑1 (MQ=255)

               cAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATAc                                                                             >  1:4306031/1‑93 (MQ=255)

               cAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATAc                                                                             >  2:2772295/1‑93 (MQ=255)

               cAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATAc                                                                             <  2:2018293/93‑1 (MQ=255)

                aaaCATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATAcc                                                                            >  1:5799852/1‑93 (MQ=255)

                 aaCATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCa                                                                           >  2:1440918/1‑93 (MQ=255)

                  aCATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAg                                                                          <  2:4855908/93‑1 (MQ=255)

                   cATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGa                                                                         >  1:1791343/1‑93 (MQ=255)

                    aTCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCATCGGGTTAGCAGTAATACCAGaa                                                                        <  2:645635/93‑1 (MQ=255)

                     tctcCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGaaa                                                                       >  1:593686/1‑93 (MQ=255)

                       tcCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGa                                                                     >  1:6290763/1‑93 (MQ=255)

                       tcCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGa                                                                     >  1:4648131/1‑93 (MQ=255)

                         cGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCATCGGGTTAGCAGTAATACCAGAAAGAcc                                                                   <  2:7200173/93‑1 (MQ=255)

                         cGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGAcc                                                                   >  2:6234801/1‑93 (MQ=255)

                           gCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGt                                                                 <  2:6763204/93‑1 (MQ=255)

                           gCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGCTAGCAGTAATACCAGAAAGACCGt                                                                 >  2:3303973/1‑93 (MQ=255)

                            cACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTc                                                                <  2:1346448/93‑1 (MQ=255)

                            cACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTc                                                                <  2:941388/93‑1 (MQ=255)

                             aCTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCATCGGGTTAGCAGTAATACCAGAAAGACCGTCa                                                               <  2:3917257/93‑1 (MQ=255)

                             aCTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCa                                                               >  2:5955933/1‑93 (MQ=255)

                             aCTTCGATCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCATCGGGTTAGCAGTAATACCAGAAAGACCGTCa                                                               <  2:4242836/93‑1 (MQ=255)

                               ttCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGa                                                             >  1:4508360/1‑93 (MQ=255)

                               ttCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGa                                                             >  2:4386997/1‑93 (MQ=255)

                                 cGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCATCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAac                                                           <  2:3629446/93‑1 (MQ=255)

                                  ggTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAacc                                                          >  1:1963206/1‑93 (MQ=255)

                                   gTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCATCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAacca                                                         <  2:269191/93‑1 (MQ=255)

                                   gTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAacca                                                         >  2:1209316/1‑93 (MQ=255)

                                   gTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAacca                                                         <  2:1146566/93‑1 (MQ=255)

                                   gTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAacca                                                         <  1:6089149/93‑1 (MQ=255)

                                   gTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAacca                                                         <  1:5491649/93‑1 (MQ=255)

                                    tCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAaccac                                                        >  2:3811179/1‑93 (MQ=255)

                                    tCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAaccac                                                        <  1:3628769/93‑1 (MQ=255)

                                    tCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAaccac                                                        <  2:2153938/93‑1 (MQ=255)

                                    tCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAaccac                                                        <  1:7284751/93‑1 (MQ=255)

                                     cAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAaccacc                                                       <  1:380815/93‑1 (MQ=255)

                                     cAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAaccacc                                                       <  1:4146181/93‑1 (MQ=255)

                                      aGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCg                                                      >  1:5602770/1‑93 (MQ=255)

                                       gTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCATCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGc                                                     >  1:308671/1‑93 (MQ=255)

                                       gTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGc                                                     >  2:6643242/1‑93 (MQ=255)

                                        tACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCa                                                    <  2:5356661/93‑1 (MQ=255)

                                            gtagtaCCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTcc                                                <  2:3094458/93‑1 (MQ=255)

                                               gtaCCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGa                                             >  1:2749537/1‑93 (MQ=255)

                                                taCCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGAc                                            >  2:696605/1‑93 (MQ=255)

                                                     ccgTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAAc                                       >  2:1749255/1‑93 (MQ=255)

                                                      cgTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACt                                      >  2:729620/1‑93 (MQ=255)

                                                      cgTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACt                                      >  1:2966143/1‑93 (MQ=255)

                                                      cgTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACt                                      <  1:3314790/93‑1 (MQ=255)

                                                        ttAGCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTc                                    >  1:3243034/1‑93 (MQ=255)

                                                           gCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGt                                 >  1:6895781/1‑93 (MQ=255)

                                                           gCAATCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGt                                 >  1:2463211/1‑93 (MQ=255)

                                                            cAATCACGTAGTCAGAGAAACGCCCCAGCATCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGtt                                <  1:5513875/93‑1 (MQ=255)

                                                              aTCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCg                              >  1:7427636/1‑93 (MQ=255)

                                                               tCACGTAGTCAGAGAAACGCCCCAGCATCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGc                             <  2:5074242/93‑1 (MQ=255)

                                                               tCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGc                             <  1:1401360/93‑1 (MQ=255)

                                                               tCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGc                             >  2:4611146/1‑93 (MQ=255)

                                                               tCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGc                             <  2:344716/93‑1 (MQ=255)

                                                               tCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGc                             >  2:575225/1‑93 (MQ=255)

                                                               tCACGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTACAGACCAAACTTCAGTTCGc                             <  2:6934246/93‑1 (MQ=255)

                                                                 aCGTAGTCAGAGAAACGCCCCAGCATCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTg                           <  1:1038104/93‑1 (MQ=255)

                                                                 aCGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTg                           <  1:3394825/93‑1 (MQ=255)

                                                                 aCGTAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTg                           <  1:424999/93‑1 (MQ=255)

                                                                    tAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGt                        >  1:2610744/1‑93 (MQ=255)

                                                                    tAGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGt                        <  2:2502200/93‑1 (MQ=255)

                                                                     aGTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGtt                       <  2:5005458/93‑1 (MQ=255)

                                                                      gTCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGttt                      >  1:6116657/1‑93 (MQ=255)

                                                                       tCAGAGAAACGCCCCAGCATCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTg                     >  2:4045142/1‑93 (MQ=255)

                                                                       tCAGAGAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTg                     >  2:529959/1‑93 (MQ=255)

                                                                         agagAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGcc                   <  1:154936/93‑1 (MQ=255)

                                                                         agagAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGcc                   <  1:1110680/93‑1 (MQ=255)

                                                                         agagAAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGcc                   >  2:4095949/1‑93 (MQ=255)

                                                                          gagaAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCt                  <  2:7491663/93‑1 (MQ=255)

                                                                          gagaAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCt                  <  2:379046/93‑1 (MQ=255)

                                                                          gagaAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCt                  >  1:789884/1‑93 (MQ=255)

                                                                          gagaAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCt                  >  1:5233241/1‑93 (MQ=255)

                                                                           agaAACGCCCCAGCATCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCTg                 <  1:1181667/93‑1 (MQ=255)

                                                                           agaAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCTg                 <  2:7398014/93‑1 (MQ=255)

                                                                           agaAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCTg                 <  1:6776387/93‑1 (MQ=255)

                                                                            gaAACGCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCTgg                <  1:318776/93‑1 (MQ=255)

                                                                             aaaCGCCCCAGCACCGGGTTAGCAGTAATTCCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCTGGc               <  2:4686081/93‑1 (MQ=255)

                                                                                 gCCCCAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCTGGCTcgc           <  1:3557642/93‑1 (MQ=255)

                                                                                    ccAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCTGGCTCGCGtt        >  2:2842137/1‑93 (MQ=255)

                                                                                    ccAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCTGGCTCGCGtt        >  1:6474175/1‑93 (MQ=255)

                                                                                     cAGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCTGGCTCGCGttt       >  1:2099797/1‑93 (MQ=255)

                                                                                      aGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCTGGCTCGCGTTTa      <  2:6838102/93‑1 (MQ=255)

                                                                                      aGCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCTGGCTCGCGTTTa      <  2:3339252/93‑1 (MQ=255)

                                                                                       gCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCTGGCTCGCGTTTAt     <  2:5696945/93‑1 (MQ=255)

                                                                                       gCACCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCTGGCTCGCGTTTAt     <  1:3163013/93‑1 (MQ=255)

                                                                                         aTCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCTGGCTCGCGTTTATCg   <  1:4463162/93‑1 (MQ=255)

GCAACTGCTCTGCGCCAAACATCTCCGGCACTTCGGTCAGTACGGTAGTACCGCCGTTAGCAATCACGTAGTCAGAGAAACGCCCCAGCATCGGGTTAGCAGTAATACCAGAAAGACCGTCAGAACCACCGCACTCCAGACCAAACTTCAGTTCGCTGAGTTTGCCTGGCTCGCGTTTATCGT  >  NC_000913/3240033‑3240215

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: `‑`