BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,015,060	G→A	D670N (GAC→AAC)	metE →	cobalamin‑independent homocysteine transmethylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,015,060	0	G	A	100.0%	1610.8 / NA	370	D670N (GAC→AAC)	metE	cobalamin‑independent homocysteine transmethylase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (177/193); total (177/193)

Only 100 of 370 total aligned reads displayed.

AAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGA  >  NC_000913/4014962‑4015157

aaaTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAAc                                                                                                 >  1:7467732/1‑101 (MQ=255)

aaaTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAAc                                                                                                 <  2:7467009/101‑1 (MQ=255)

 aaTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACa                                                                                                <  2:7467267/101‑1 (MQ=255)

 aaTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACa                                                                                                <  2:7466955/101‑1 (MQ=255)

 aaTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACa                                                                                                <  1:7467190/101‑1 (MQ=255)

  aTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACAt                                                                                               >  1:7467736/1‑101 (MQ=255)

  aTCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACAt                                                                                               <  1:7467304/101‑1 (MQ=255)

    ccACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATgg                                                                                             <  1:7467148/101‑1 (MQ=255)

    ccACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATgg                                                                                             <  2:7467260/101‑1 (MQ=255)

     cacaCTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGa                                                                                            >  1:7467741/1‑101 (MQ=255)

     cacaCTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGa                                                                                            >  2:7467739/1‑101 (MQ=255)

       cacTCACGTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGt                                                                                          >  1:7467745/1‑101 (MQ=255)

        acTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGtt                                                                                         <  2:7467121/101‑1 (MQ=255)

        acTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGtt                                                                                         <  2:7466732/101‑1 (MQ=255)

         cTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTg                                                                                        <  2:7467090/101‑1 (MQ=255)

          tCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGc                                                                                       >  2:7467749/1‑101 (MQ=255)

             caTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTgg                                                                                    <  1:7467326/101‑1 (MQ=255)

              aTGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGa                                                                                   >  1:7467753/1‑101 (MQ=255)

                 tgtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTc                                                                                <  1:7467309/101‑1 (MQ=255)

                 tgtTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTc                                                                                >  1:7467757/1‑101 (MQ=255)

                   ttATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGt                                                                              >  2:7467761/1‑101 (MQ=255)

                   ttATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGt                                                                              <  2:7467192/101‑1 (MQ=255)

                    tATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGtt                                                                             <  1:7467113/101‑1 (MQ=255)

                      ttGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAg                                                                                  >  2:7467766/1‑94 (MQ=255)

                      ttGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGt                                                                                 >  1:7467765/1‑95 (MQ=255)

                         cGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAag                                                                        >  2:7467773/1‑101 (MQ=255)

                         cGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAag                                                                        >  1:7467775/1‑101 (MQ=255)

                         cGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAag                                                                        <  1:7467372/101‑1 (MQ=255)

                         cGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAag                                                                        <  2:7467336/101‑1 (MQ=255)

                          gAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAaga                                                                       >  1:7467778/1‑101 (MQ=255)

                          gAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAaga                                                                       <  1:7467252/101‑1 (MQ=255)

                           aGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAagag                                                                      <  1:7467562/101‑1 (MQ=255)

                            gTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGt                                                                     >  1:7467783/1‑101 (MQ=255)

                             ttCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGtt                                                                    <  1:7467400/101‑1 (MQ=255)

                              tCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGttt                                                                   >  2:7467785/1‑101 (MQ=255)

                              tCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGttt                                                                   <  1:7467224/101‑1 (MQ=255)

                                aaCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGa                                                                 >  1:7467789/1‑101 (MQ=255)

                                aaCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGa                                                                 <  2:7467229/101‑1 (MQ=255)

                                 aCGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGAt                                                                <  2:7467543/101‑1 (MQ=255)

                                  cGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGAtt                                                               >  1:7467794/1‑101 (MQ=255)

                                  cGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGAtt                                                               <  2:7467471/101‑1 (MQ=255)

                                   gACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTa                                                              <  1:7467126/101‑1 (MQ=255)

                                   gACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTa                                                              <  1:7467075/101‑1 (MQ=255)

                                     catcatGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATc                                                            >  1:7467798/1‑101 (MQ=255)

                                     catcatGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATc                                                            <  2:7467509/101‑1 (MQ=255)

                                      atcatGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATcc                                                           >  2:7467802/1‑101 (MQ=255)

                                       tcatGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCa                                                          >  1:7467805/1‑101 (MQ=255)

                                       tcatGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCa                                                          <  1:7467157/101‑1 (MQ=255)

                                         atGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCaaa                                                        >  2:7467809/1‑101 (MQ=255)

                                         atGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCaaa                                                        <  1:7466545/101‑1 (MQ=255)

                                          tGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAAt                                                       <  2:7467249/101‑1 (MQ=255)

                                            gATTCGATTGCGGCGCTTGACGCAGACGTCATCTCCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTg                                                                  <  1:7467007/88‑1 (MQ=255)

                                            gATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGa                                                     >  1:7467813/1‑101 (MQ=255)

                                             aTTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGaa                                                    <  2:7467144/101‑1 (MQ=255)

                                                cGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATc                                                 >  2:7467817/1‑101 (MQ=255)

                                                cGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATc                                                 <  1:7467271/101‑1 (MQ=255)

                                                  aTTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCgg                                               >  1:7467820/1‑101 (MQ=255)

                                                  aTTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCgg                                               <  2:7467169/101‑1 (MQ=255)

                                                   ttGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGt                                              >  2:7467824/1‑101 (MQ=255)

                                                    tGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTc                                             <  1:7467346/101‑1 (MQ=255)

                                                    tGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTc                                             >  1:7467826/1‑101 (MQ=255)

                                                      cggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCt                                           >  2:7467832/1‑101 (MQ=255)

                                                       ggcgCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTg                                          >  1:7467835/1‑101 (MQ=255)

                                                         cgcTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGc                                        <  2:7467390/101‑1 (MQ=255)

                                                          gcTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCg                                       >  1:7467838/1‑101 (MQ=255)

                                                             ggACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCt                                    <  1:7467383/101‑1 (MQ=255)

                                                              gACGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTa                                   >  2:7467841/1‑101 (MQ=255)

                                                               aCGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTAt                                  <  1:7467281/101‑1 (MQ=255)

                                                               aCGCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTAt                                  <  1:7467442/101‑1 (MQ=255)

                                                                 gCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGa                                >  2:7467846/1‑101 (MQ=255)

                                                                 gCAGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGa                                >  1:7467848/1‑101 (MQ=255)

                                                                   aGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACa                              >  1:7467852/1‑101 (MQ=255)

                                                                   aGACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACa                              <  2:7467424/101‑1 (MQ=255)

                                                                    gACGTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACAt                             <  1:7467416/101‑1 (MQ=255)

                                                                       gTCATCACCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCa                          >  1:7467859/1‑101 (MQ=255)

                                                                        tcatcaCCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCAc                         >  2:7467861/1‑101 (MQ=255)

                                                                        tcatcaCCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCAc                         <  1:7467176/101‑1 (MQ=255)

                                                                          atcaCCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTc                       >  1:7467864/1‑101 (MQ=255)

                                                                           tcaCCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCg                      <  2:7467398/101‑1 (MQ=255)

                                                                            caCCATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGc                     >  1:7467869/1‑101 (MQ=255)

                                                                              ccATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCa                   <  1:7467200/101‑1 (MQ=255)

                                                                               cATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCaa                  <  1:7467460/101‑1 (MQ=255)

                                                                               cATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCaa                  >  2:7467875/1‑101 (MQ=255)

                                                                               cATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCaa                  <  1:7467634/101‑1 (MQ=255)

                                                                               cATCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCaa                  >  1:7467874/1‑101 (MQ=255)

                                                                                aTCGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCaaa                 <  1:7467501/101‑1 (MQ=255)

                                                                                  cGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACg               <  2:7467362/101‑1 (MQ=255)

                                                                                  cGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGGTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACg               >  1:7467881/1‑101 (MQ=255)

                                                                                  cGAAACCTCGCGTTCCAACATGGAGTTGCTGGAGCCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACg               >  2:7467880/1‑101 (MQ=255)

                                                                                    aaaCCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTa             >  1:7467888/1‑101 (MQ=255)

                                                                                      aCCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTAcc           >  2:7467889/1‑101 (MQ=255)

                                                                                      aCCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTAcc           <  2:7467316/101‑1 (MQ=255)

                                                                                      aCCTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTAcc           <  1:7467215/101‑1 (MQ=255)

                                                                                        cTCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATg                                                      <  2:7467036/56‑1 (MQ=255)

                                                                                         tCGCGTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAg        <  2:7467297/101‑1 (MQ=255)

                                                                                           gcgTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCg      >  1:7467895/1‑101 (MQ=255)

                                                                                             gTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTg    <  1:7467298/101‑1 (MQ=255)

                                                                                             gTTCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTg    >  1:7467898/1‑101 (MQ=255)

                                                                                              ttCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTgg   >  1:7467902/1‑101 (MQ=255)

                                                                                               tCCAACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGa  >  2:7467905/1‑101 (MQ=255)

AAATCCACACTCACATGTGTTATTGCGAGTTCAACGACATCATGGATTCGATTGCGGCGCTGGACGCAGACGTCATCACCATCGAAACCTCGCGTTCCGACATGGAGTTGCTGGAGTCGTTTGAAGAGTTTGATTATCCAAATGAAATCGGTCCTGGCGTCTATGACATTCACTCGCCAAACGTACCGAGCGTGGA  >  NC_000913/4014962‑4015157

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