BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,075,081	G→A	P39S (CCG→TCG)	fumE ←	fumarase E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,075,081	0	G	A	100.0%	1623.6 / NA	378	P39S (CCG→TCG)	fumE	fumarase E
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (168/210); total (168/210)

Only 100 of 383 total aligned reads displayed.

CAGTGCAACATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGTTGCTCAAGCACATCATCT  >  NC_000913/3074983‑3075178

  gTGCAACATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGGCGACAGa                                                                                               <  2:5655812/101‑1 (MQ=255)

    gCAACATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGaaa                                                                                             >  2:5656334/1‑101 (MQ=255)

    gCAACATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGaa                                                                                              <  2:5655868/100‑1 (MQ=255)

     cAACATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAAc                                                                                            >  1:5656336/1‑101 (MQ=255)

      aaCATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACt                                                                                           <  2:5655655/101‑1 (MQ=255)

       aCATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTg                                                                                          <  1:5655412/101‑1 (MQ=255)

        cATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTgg                                                                                         <  2:5655885/101‑1 (MQ=255)

         aTCGATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGc                                                                                        >  1:5656344/1‑101 (MQ=255)

         aTCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGc                                                                                        <  1:5655663/101‑1 (MQ=255)

           cAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACCGAAACTGGCGt                                                                                      >  1:5656350/1‑101 (MQ=255)

           cAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGt                                                                                      >  2:5656351/1‑101 (MQ=255)

            aaTATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTa                                                                                     <  1:5655959/101‑1 (MQ=255)

              tatCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAAc                                                                                            >  2:5656359/1‑92 (MQ=255)

               atCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATcc                                                                                  >  1:5656362/1‑101 (MQ=255)

                tCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGTCAAACAGCGACGACAGAAACTGGCGTATccc                                                                                 <  1:5655930/101‑1 (MQ=255)

                  gTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACACCGACGACAGaa                                                                                              <  2:5655791/86‑1 (MQ=255)

                   tCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTg                                                                              <  1:5655754/101‑1 (MQ=255)

                    cAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTgg                                                                             >  2:5656370/1‑101 (MQ=255)

                     aaGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGa                                                                            <  2:5655607/101‑1 (MQ=255)

                     aaGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGa                                                                            >  1:5656374/1‑101 (MQ=255)

                      aGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAg                                                                           <  1:5655920/101‑1 (MQ=255)

                      aGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAg                                                                           >  2:5656376/1‑101 (MQ=255)

                       gCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGc                                                                          <  1:5655817/101‑1 (MQ=255)

                        cGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGAAAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCa                                                                         >  2:5656379/1‑101 (MQ=255)

                         gggCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGa                                                                            >  1:5656383/1‑97 (MQ=255)

                         gggCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCaa                                                                        <  1:5655892/101‑1 (MQ=255)

                          ggCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCaaa                                                                       >  2:5656385/1‑101 (MQ=255)

                          ggCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCaaa                                                                       <  2:5655749/101‑1 (MQ=255)

                            ccGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAAt                                                                     <  2:5655760/101‑1 (MQ=255)

                             cGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATg                                                                    <  1:5656023/101‑1 (MQ=255)

                             cGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATg                                                                    >  2:5656391/1‑101 (MQ=255)

                              gCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATgg                                                                   >  1:5656395/1‑101 (MQ=255)

                              gCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATgg                                                                   <  1:5655509/101‑1 (MQ=255)

                               cTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGt                                                                  <  1:5655579/101‑1 (MQ=255)

                                tttGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTa                                                                 >  2:5656399/1‑101 (MQ=255)

                                 ttGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTAt                                                                <  1:5655800/101‑1 (MQ=255)

                                 ttGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTAt                                                                <  2:5655991/101‑1 (MQ=255)

                                   gggCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATgc                                                              <  2:5655860/101‑1 (MQ=255)

                                     gCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATgcgc                                                            <  1:5655740/101‑1 (MQ=255)

                                     gCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATgcgc                                                            <  2:5656128/101‑1 (MQ=255)

                                      cGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCa                                                           <  2:5655497/101‑1 (MQ=255)

                                        aGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGt                                                         >  2:5656406/1‑101 (MQ=255)

                                         gTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGtt                                                        >  1:5656409/1‑101 (MQ=255)

                                            aCGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAAt                                                                     <  1:5655193/85‑1 (MQ=255)

                                            aCGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTc                                                     >  1:5656415/1‑101 (MQ=255)

                                            aCGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTc                                                     >  2:5656413/1‑101 (MQ=255)

                                             cGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCa                                                    <  2:5655777/101‑1 (MQ=255)

                                             cGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCa                                                    >  2:5656418/1‑101 (MQ=255)

                                                cTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATaa                                                 >  2:5656424/1‑101 (MQ=255)

                                                cTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATaa                                                 >  1:5656426/1‑101 (MQ=255)

                                                  tCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAAt                                               <  2:5655975/101‑1 (MQ=255)

                                                   cACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATg                                              <  2:5656045/101‑1 (MQ=255)

                                                   cACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATg                                              <  2:5655927/101‑1 (MQ=255)

                                                   cACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATg                                              >  1:5656433/1‑101 (MQ=255)

                                                   cACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATg                                              <  2:5655824/101‑1 (MQ=255)

                                                   cACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATg                                              <  1:5656074/101‑1 (MQ=255)

                                                    aCTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGa                                                                            <  1:5655977/70‑1 (MQ=255)

                                                     cTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGaa                                            <  1:5655765/101‑1 (MQ=255)

                                                     cTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGaa                                            <  2:5656114/101‑1 (MQ=255)

                                                      tGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGaaa                                           <  2:5655633/101‑1 (MQ=255)

                                                      tGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGaaa                                           >  1:5656437/1‑101 (MQ=255)

                                                        cATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAaaat                                         >  2:5656439/1‑101 (MQ=255)

                                                         atatTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAaaata                                        >  1:5656442/1‑101 (MQ=255)

                                                          tatTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAaaataa                                       <  2:5655539/101‑1 (MQ=255)

                                                           atTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAaaataaa                                      <  2:5655730/101‑1 (MQ=255)

                                                           atTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAaaataaa                                      <  1:5655703/101‑1 (MQ=255)

                                                            ttCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAaaataaat                                     <  2:5655280/101‑1 (MQ=255)

                                                             tCCACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAAtt                                    <  2:5656025/101‑1 (MQ=255)

                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTg                                   <  2:5655783/101‑1 (MQ=255)

                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTg                                   <  2:5655909/101‑1 (MQ=255)

                                                              ccACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTg                                   >  1:5656449/1‑101 (MQ=255)

                                                               cACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGt                                  <  2:5655954/101‑1 (MQ=255)

                                                               cACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGt                                  <  2:5656195/101‑1 (MQ=255)

                                                               cACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGt                                  <  1:5656203/101‑1 (MQ=255)

                                                               cACGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGt                                  >  2:5656450/1‑101 (MQ=255)

                                                                aCGATCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTc                                 >  1:5656455/1‑101 (MQ=255)

                                                                   aTCTCTTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCtg                              >  2:5656458/1‑101 (MQ=255)

                                                                      tctTCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGc                           >  2:5656464/1‑101 (MQ=255)

                                                                        ttCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGt                                                                                      <  1:5655680/40‑1 (MQ=255)

                                                                        ttCATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCAt                         <  1:5656121/101‑1 (MQ=255)

                                                                          cATCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATcc                       >  2:5656468/1‑101 (MQ=255)

                                                                           aTCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCa                      >  1:5656469/1‑101 (MQ=255)

                                                                           aTCGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCa                      <  1:5655863/101‑1 (MQ=255)

                                                                             cGTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGt                    >  1:5656476/1‑101 (MQ=255)

                                                                              gTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGtt                   >  1:5656480/1‑101 (MQ=255)

                                                                              gTTATCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGtt                   >  2:5656479/1‑101 (MQ=255)

                                                                                 aTCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGTTGCt                >  1:5656488/1‑101 (MQ=255)

                                                                                 aTCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGTTGCt                >  2:5656487/1‑101 (MQ=255)

                                                                                  tCGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGTTGCTc               <  2:5655691/101‑1 (MQ=255)

                                                                                   cGACAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGTTGCTCa              >  2:5656492/1‑101 (MQ=255)

                                                                                     aCAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGTTGCTCAAg            >  2:5656497/1‑101 (MQ=255)

                                                                                     aCAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGTTGCTCAAg            <  1:5655839/101‑1 (MQ=255)

                                                                                      cAAACAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGTTGCTCAAGc           <  1:5656008/101‑1 (MQ=255)

                                                                                       aaaCAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGTTGCTCAAGca          >  1:5656500/1‑101 (MQ=255)

                                                                                          cAGCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGTTGCTCAAGCAcat       <  2:5655701/101‑1 (MQ=255)

                                                                                            gCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGTTGCTCAAGCAcatca     <  2:5656319/101‑1 (MQ=255)

                                                                                            gCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGTTGCTCAAGCAcatca     >  1:5656506/1‑101 (MQ=255)

                                                                                            gCGACGACAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGTTGCTCAAGCAcatca     <  2:5655802/101‑1 (MQ=255)

                                                                                              gacgacAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGTTGCTCAAGCAcatcat    <  2:5655715/100‑1 (MQ=255)

                                                                                              gacgacAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGTTGCTCAAGCACatcatc   <  1:5655875/101‑1 (MQ=255)

CAGTGCAACATCAATATCGTCAAGCGGGCCGCTTTGGGCGAGTAACGGCTTCACTGCATATTCCACGATCTCTTCATCGTTATCGACAAACAGCGACGGCAGAAACTGGCGTATCCCCTGGAGCAAAATGGTATGCGCAGTTTTCATAAATGAAAATAAATTGTCCTGTGCATCCAGTTGCTCAAGCACATCATCT  >  NC_000913/3074983‑3075178

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