BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,650,575	G→A	pseudogene (498/921 nt)	ydfE →	Qin prophage; protein YdfE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,650,575	0	G	A	100.0%	1038.8 / NA	250	pseudogene (498/921 nt)	ydfE	Qin prophage; protein YdfE
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (124/126); total (124/126)

Only 100 of 250 total aligned reads displayed.

AGCGCACAGGAATTGTTGCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAGCTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATGCGCATCTGGAAAAAACGTC  >  NC_000913/1650477‑1650672

aGCGCACAGGAATTGTTGCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACt                                                                                                 <  2:2157423/101‑1 (MQ=255)

 gcgcACAGGAATTGTTGCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACtt                                                                                                <  2:2157472/101‑1 (MQ=255)

  cgcACAGGAATTGTTGCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTa                                                                                               <  1:2157278/101‑1 (MQ=255)

  cgcACAGGAATTGTTGCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTa                                                                                               <  2:2157271/101‑1 (MQ=255)

   gcACAGGAATTGTTGCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTaa                                                                                              <  1:2157383/101‑1 (MQ=255)

     acaGGAATTGTTGCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGt                                                                                            >  2:2157747/1‑101 (MQ=255)

      caGGAATTGTTGCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGtt                                                                                           >  2:2157749/1‑101 (MQ=255)

        ggAATTGTTGCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGc                                                                                         >  2:2157752/1‑101 (MQ=255)

        ggAATTGTTGCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGc                                                                                         <  2:2157294/101‑1 (MQ=255)

           aTTGTTGCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAg                                                                                             <  1:2157440/94‑1 (MQ=255)

           aTTGTTGCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACg                                                                                      >  1:2157757/1‑101 (MQ=255)

             tgttgCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGcc                                                                                    >  2:2157758/1‑101 (MQ=255)

              gttgCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCa                                                                                   <  1:2157546/101‑1 (MQ=255)

                tgCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCAtt                                                                                 >  2:2157761/1‑101 (MQ=255)

                tgCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCAtt                                                                                 <  2:2157483/101‑1 (MQ=255)

                 gCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATtg                                                                                <  1:2157444/101‑1 (MQ=255)

                 gCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATtg                                                                                <  2:2157223/101‑1 (MQ=255)

                   aGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATtggg                                                                              <  2:2157432/101‑1 (MQ=255)

                   aGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATtggg                                                                              >  1:2157763/1‑101 (MQ=255)

                     gAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATtgggtg                                                                            >  2:2157766/1‑101 (MQ=255)

                     gAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATtgggtg                                                                            >  1:2157767/1‑101 (MQ=255)

                     gAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATtgggtg                                                                            <  1:2157108/101‑1 (MQ=255)

                      aaGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATtgggtgg                                                                           <  1:2157584/101‑1 (MQ=255)

                      aaGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATtgggtgg                                                                           >  2:2157770/1‑101 (MQ=255)

                      aaGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATtgggtgg                                                                           <  1:2157429/101‑1 (MQ=255)

                      aaGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATtgggtgg                                                                           <  1:2157398/101‑1 (MQ=255)

                        gACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTgggtgggt                                                                         <  1:2157334/101‑1 (MQ=255)

                        gACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTgggtgggt                                                                         >  1:2157772/1‑101 (MQ=255)

                        gACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTgggtgggt                                                                         <  1:2157470/101‑1 (MQ=255)

                         aCATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGtt                                                                        <  1:2157535/101‑1 (MQ=255)

                         aCATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGtt                                                                        <  2:2157321/101‑1 (MQ=255)

                          cATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTa                                                                       <  1:2157511/101‑1 (MQ=255)

                            tGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAg                                                                     >  2:2157775/1‑101 (MQ=255)

                             ggCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGc                                                                    <  1:2157415/101‑1 (MQ=255)

                                aaCATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATa                                                                 >  2:2157780/1‑101 (MQ=255)

                                aaCATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATa                                                                 >  1:2157779/1‑101 (MQ=255)

                                 aCATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATaa                                                                >  1:2157781/1‑101 (MQ=255)

                                 aCATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATaa                                                                <  2:2157379/101‑1 (MQ=255)

                                  cATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATtgg                                                                               >  2:2157783/1‑85 (MQ=255)

                                   aTTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAAt                                                              <  1:2157424/101‑1 (MQ=255)

                                   aTTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAAt                                                              <  1:2157617/101‑1 (MQ=255)

                                     tCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGt                                                            >  2:2157786/1‑101 (MQ=255)

                                      cGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTa                                                           <  2:2157454/101‑1 (MQ=255)

                                        aCGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTAtt                                                         >  2:2157790/1‑101 (MQ=255)

                                        aCGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTAtt                                                         >  1:2157791/1‑101 (MQ=255)

                                            aCATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTg                                                     <  1:2157352/101‑1 (MQ=255)

                                            aCATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTg                                                     <  1:2157433/101‑1 (MQ=255)

                                             cATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTgg                                                    <  2:2157394/101‑1 (MQ=255)

                                              atatCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGt                                                   >  1:2157794/1‑101 (MQ=255)

                                              atatCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGt                                                   <  2:2157371/101‑1 (MQ=255)

                                               tatCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTg                                                  >  2:2157797/1‑101 (MQ=255)

                                               tatCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTg                                                  <  2:2157510/101‑1 (MQ=255)

                                               tatCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTg                                                  <  2:2157627/101‑1 (MQ=255)

                                                  cGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTa                                               >  2:2157799/1‑101 (MQ=255)

                                                   gTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTaa                                              <  2:2157632/101‑1 (MQ=255)

                                                    ttGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTaaa                                             >  1:2157801/1‑101 (MQ=255)

                                                    ttGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCAAAAATGTATTCCTGGTGCTaaa                                             <  1:2157449/101‑1 (MQ=255)

                                                     tGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAAt                                            >  2:2157805/1‑101 (MQ=255)

                                                      gAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATg                                           <  1:2157475/101‑1 (MQ=255)

                                                      gAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATg                                           <  1:2157452/101‑1 (MQ=255)

                                                      gAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATg                                           >  1:2157806/1‑101 (MQ=255)

                                                       aaGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATgg                                          <  2:2157499/101‑1 (MQ=255)

                                                           gTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCg                                      <  2:2157189/101‑1 (MQ=255)

                                                           gTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCg                                      <  2:2157531/101‑1 (MQ=255)

                                                           gTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCg                                      >  2:2157810/1‑101 (MQ=255)

                                                             tgatgaAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGaa                                    <  1:2157689/101‑1 (MQ=255)

                                                              gatgaAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGaaa                                   >  1:2157814/1‑101 (MQ=255)

                                                              gatgaAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGaaa                                   >  2:2157813/1‑101 (MQ=255)

                                                               atgaAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAAt                                  >  2:2157816/1‑101 (MQ=255)

                                                                tgaAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAAtt                                 <  1:2157250/101‑1 (MQ=255)

                                                                tgaAAATGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAAtt                                 <  1:2157495/101‑1 (MQ=255)

                                                                  aaaaTGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGt                                                                                            >  2:2157820/1‑40 (MQ=255)

                                                                   aaaTGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAg                              <  2:2157670/101‑1 (MQ=255)

                                                                   aaaTGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAg                              >  1:2157821/1‑101 (MQ=255)

                                                                    aaTGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAgg                             >  2:2157824/1‑101 (MQ=255)

                                                                     aTGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGc                            >  2:2157829/1‑101 (MQ=255)

                                                                     aTGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGc                            >  1:2157828/1‑101 (MQ=255)

                                                                     aTGCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGc                            <  1:2157309/101‑1 (MQ=255)

                                                                       gCCAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCag                          >  2:2157831/1‑101 (MQ=255)

                                                                        ccAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCaga                         <  2:2157462/101‑1 (MQ=255)

                                                                         cAGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCagag                        >  2:2157834/1‑101 (MQ=255)

                                                                          aGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTTCTAAATGGCCCGAAATTCAGGCAgaga                       <  2:2157552/101‑1 (MQ=255)

                                                                          aGAAGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAgaga                       >  1:2157836/1‑101 (MQ=255)

                                                                             aGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATgc                    <  1:2157502/101‑1 (MQ=255)

                                                                             aGTTAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATgc                    >  2:2157839/1‑101 (MQ=255)

                                                                              ggtAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATgcg                   >  1:2157842/3‑101 (MQ=255)

                                                                               ttAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATgcgc                  <  1:2157522/101‑1 (MQ=255)

                                                                                tAACGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCag                          >  2:2157844/1‑92 (MQ=255)

                                                                                 aaCGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGAGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATGCGCAt                >  2:2157847/1‑101 (MQ=255)

                                                                                  aCGCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATGCGCATc               >  1:2157848/1‑101 (MQ=255)

                                                                                    gCCATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATGCGCATCTg             >  2:2157853/1‑101 (MQ=255)

                                                                                     ccATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATGCGCATCTgg            <  2:2157407/101‑1 (MQ=255)

                                                                                      cATGTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATGCGCATCTGGa           >  2:2157855/1‑101 (MQ=255)

                                                                                         gTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATGCGCATCTGGaaaa        <  1:2157629/101‑1 (MQ=255)

                                                                                         gTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATGCGCATCTGGaaaa        >  1:2157859/1‑101 (MQ=255)

                                                                                         gTATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATGCGCATCTGGaaaa        >  2:2157858/1‑101 (MQ=255)

                                                                                          tATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATGCGCATCTGGaaaaa       >  1:2157861/1‑101 (MQ=255)

                                                                                          tATCCGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATGCGCATCTGGaaaaa       >  2:2157863/1‑101 (MQ=255)

                                                                                              cGGAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATGCGCATCTGGAAAAAACGt   >  1:2157867/1‑101 (MQ=255)

                                                                                               ggAACTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATGCGCATCTGGAAAAAACGTc  <  2:2157692/101‑1 (MQ=255)

AGCGCACAGGAATTGTTGCAGGAAGACATGGCAACATTCGACGGACATATCGTTGAAGCGTTGATGAAAATGCCAGAAGTTAACGCCATGTATCCGGAGCTTAAGTTGCACGCCATTGGGTGGGTTAAGCATAAATGTATTCCTGGTGCTAAATGGCCCGAAATTCAGGCAGAGATGCGCATCTGGAAAAAACGTC  >  NC_000913/1650477‑1650672

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