BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,225,050	A→G	E1215G (GAG→GGG)	metH →	B12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,225,050	0	A	G	100.0%	542.4 / NA	170	E1215G (GAG→GGG)	metH	B12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (87/83); total (87/83)

Only 100 of 174 total aligned reads displayed.

GACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGT  >  NC_000913/4224962‑4225139

 aCAGAAAGTACTAGGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgct                                                                                      <  1:3378674/93‑1 (MQ=255)

  cAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctg                                                                                     <  1:3753214/93‑1 (MQ=255)

   aGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGTgctgg                                                                                    >  1:4684757/1‑93 (MQ=255)

    gaaaGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggc                                                                                   <  1:3602411/91‑1 (MQ=255)

    gCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggc                                                                                   <  2:560150/93‑1 (MQ=255)

      aaGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggctg                                                                                 >  2:2548337/1‑93 (MQ=255)

      aaGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggctg                                                                                 >  1:4636492/1‑93 (MQ=255)

        gTACTCCGCTGTAGCACAAATCCAGCGCGATCAGGTGGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCtggctggc                                                                               <  2:2112534/93‑1 (MQ=255)

         tactacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCa                                                                              >  1:4416039/1‑93 (MQ=255)

           ctacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTGTGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCAcc                                                                            <  1:2520454/93‑1 (MQ=255)

            tacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCg                                                                           >  1:830593/1‑93 (MQ=255)

             acGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGa                                                                          >  2:4491467/1‑93 (MQ=255)

              cGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGaa                                                                         <  1:1931846/93‑1 (MQ=255)

               gCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAAt                                                                        >  1:2624538/1‑93 (MQ=255)

                cTGTACCACAAATTCAGCGCGATCAGGTTAAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATc                                                                       <  2:1029559/93‑1 (MQ=255)

                 tGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCt                                                                      >  1:4529066/1‑93 (MQ=255)

                  taaGCAAAAATTGCGCGCGAACAGGTTGAACAGTAGGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTg                                                                     <  1:4224388/91‑1 (MQ=255)

                    tgCACAAGTCCAGCGCGATCAGGTTGAAGTTTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTggg                                                                   <  1:2829509/92‑1 (MQ=255)

                    aGCACAAATTCAGCGCGATCAGGTTAAAGATTATCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTggg                                                                   <  2:1205849/93‑1 (MQ=255)

                     gCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTgggg                                                                  <  1:2863789/93‑1 (MQ=255)

                     gCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTgggg                                                                  >  2:2939885/1‑93 (MQ=255)

                        caAATTCAGCGCGATCAGGATTAAGCTTATGCCCGGCGTACAGGTACGGGCGTAACCGACGCTGCGTGCTGGCTAGCACCGAATCAGGGGTAt                                                               >  2:1368104/1‑93 (MQ=255)

                        aaaaTTCGGCGGAATCGGTTTAAAGATGATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTAt                                                               <  1:3492139/92‑1 (MQ=255)

                         aaaTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGACGTTGGGCGCTGGCCGGCACCGAATCTGGGGGATg                                                              >  1:3232342/1‑93 (MQ=255)

                          atttCAGCGCGATCAGGTTAAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGa                                                             <  1:4142865/91‑1 (MQ=255)

                           aTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAc                                                            >  1:326433/1‑93 (MQ=255)

                            ttCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAcg                                                           >  1:4601015/1‑93 (MQ=255)

                             tCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAcgc                                                          >  1:1024642/1‑93 (MQ=255)

                               aGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCgg                                                        >  2:4097662/1‑93 (MQ=255)

                                  gcgATCGCGTTAAAGATTGTGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACt                                                     <  2:2847581/93‑1 (MQ=255)

                                     aTCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAt                                                  >  2:3286114/1‑93 (MQ=255)

                                     aTCAGGTGGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAt                                                  <  1:649790/93‑1 (MQ=255)

                                      tCATGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAtt                                                 <  2:329772/93‑1 (MQ=255)

                                      tCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAtt                                                 >  1:4648839/1‑93 (MQ=255)

                                      tCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGGATGACGCGGACTGAtt                                                 >  1:2140381/1‑93 (MQ=255)

                                        aGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTca                                               >  1:2376311/1‑93 (MQ=255)

                                         ggTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTGGGGCGCTGGCTGGCACCCAATCTGGGGTATGACGCGGACTGATTcac                                              >  1:3040448/1‑93 (MQ=255)

                                          gTTGAGGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTcaca                                             <  1:147195/93‑1 (MQ=255)

                                            tGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACaaa                                           >  1:3187655/1‑93 (MQ=255)

                                            ggAAGATTAGCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACaaa                                           <  2:862240/92‑1 (MQ=255)

                                              aaGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGCATGACGCGGACTGATTCACAAATc                                         <  2:4519573/93‑1 (MQ=255)

                                                 aTTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGt                                      >  1:1093729/1‑93 (MQ=255)

                                                  ttATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTc                                     >  1:3787309/1‑93 (MQ=255)

                                                  ttATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTc                                     <  2:1375560/93‑1 (MQ=255)

                                                  ttATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTc                                     >  2:4607195/1‑93 (MQ=255)

                                                  ttAGGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTc                                     <  1:1362872/93‑1 (MQ=255)

                                                   tATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCa                                    <  1:692010/93‑1 (MQ=255)

                                                   tATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGAGGACTGATTCACAAATTTGTCa                                    >  2:536926/1‑93 (MQ=255)

                                                    ctccccGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCAc                                   <  1:1006149/90‑1 (MQ=255)

                                                    aggCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCAc                                   <  1:1861525/91‑1 (MQ=255)

                                                     tGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACt                                  <  1:3008586/93‑1 (MQ=255)

                                                      gCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACtt                                 <  1:1624614/93‑1 (MQ=255)

                                                       cccGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACttt                                >  2:4419208/1‑93 (MQ=255)

                                                       cccGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACttt                                <  2:1919293/93‑1 (MQ=255)

                                                        ccgccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACtttt                               >  2:720691/1‑93 (MQ=255)

                                                        ccgccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACtttt                               <  1:142801/93‑1 (MQ=255)

                                                         cgccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTc                              <  2:2674476/93‑1 (MQ=255)

                                                         cgccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTc                              <  1:3728536/93‑1 (MQ=255)

                                                          gccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTcc                             >  2:1818379/1‑93 (MQ=255)

                                                          gccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTcc                             <  1:2460842/93‑1 (MQ=255)

                                                          gccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTcc                             <  2:2290970/93‑1 (MQ=255)

                                                          gccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTcc                             >  2:1509792/1‑93 (MQ=255)

                                                          gccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTcc                             <  2:3873256/93‑1 (MQ=255)

                                                           ccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCt                            >  2:1048709/1‑93 (MQ=255)

                                                           ccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCt                            <  2:338369/93‑1 (MQ=255)

                                                            cgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCtt                           >  2:679928/1‑93 (MQ=255)

                                                             gTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTa                          >  1:430539/1‑93 (MQ=255)

                                                              tAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTac                         >  2:169910/1‑93 (MQ=255)

                                                                 aGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTacaac                      >  2:4306669/1‑93 (MQ=255)

                                                                 aGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTacaac                      <  2:3027063/93‑1 (MQ=255)

                                                                   gTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTAcaacaa                    <  2:2025565/93‑1 (MQ=255)

                                                                    tATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaa                   >  2:1695910/1‑93 (MQ=255)

                                                                    tATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaa                   >  1:977517/1‑93 (MQ=255)

                                                                    tATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaa                   >  1:2029934/1‑93 (MQ=255)

                                                                     aTGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaac                  <  2:2507548/93‑1 (MQ=255)

                                                                     aTGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaac                  <  2:3630351/93‑1 (MQ=255)

                                                                     aTGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaac                  >  1:1187108/1‑93 (MQ=255)

                                                                      tGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaaca                 >  2:2950850/1‑93 (MQ=255)

                                                                      tGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaaca                 >  2:3580986/1‑93 (MQ=255)

                                                                      tGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaaca                 <  1:563728/93‑1 (MQ=255)

                                                                       gAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAg                >  2:3416188/1‑93 (MQ=255)

                                                                       gAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAg                <  1:1246871/93‑1 (MQ=255)

                                                                        aGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAgg               <  2:3578347/93‑1 (MQ=255)

                                                                        aGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAgg               >  1:1208670/1‑93 (MQ=255)

                                                                        aGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAgg               <  2:4164962/93‑1 (MQ=255)

                                                                          cGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgc             <  2:47010/93‑1 (MQ=255)

                                                                            ttACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcgc           >  2:301509/1‑93 (MQ=255)

                                                                            ttACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcgc           >  2:2844847/1‑93 (MQ=255)

                                                                            ttACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcgc           >  2:1403479/1‑93 (MQ=255)

                                                                              aCCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTc         <  1:4468964/93‑1 (MQ=255)

                                                                               ccGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCa        >  2:2353650/1‑93 (MQ=255)

                                                                               ccGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCa        >  2:3706473/1‑93 (MQ=255)

                                                                                cGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCaa       <  2:898088/93‑1 (MQ=255)

                                                                                   aGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa    <  1:1181755/93‑1 (MQ=255)

                                                                                   aGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa    <  1:2956756/93‑1 (MQ=255)

                                                                                   aGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa    <  1:3604297/93‑1 (MQ=255)

                                                                                    gTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGGGCTCAATGAg   >  2:2477195/1‑93 (MQ=255)

                                                                                    gTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg   <  2:3848834/93‑1 (MQ=255)

                                                                                    gTTGGGCGCTGGCGGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg   <  2:3200644/93‑1 (MQ=255)

                                                                                     ttGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAgc  >  1:2993249/1‑92 (MQ=255)

GACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGT  >  NC_000913/4224962‑4225139

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