BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,221,526	A→G	E40E (GAA→GAG)	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,221,526	0	A	G	100.0%	672.8 / NA	214	E40E (GAA→GAG)	metH	B12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (97/117); total (97/117)

Only 100 of 228 total aligned reads displayed.

GCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCCACAACGC  >  NC_000913/4221436‑4221615

 cAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAgcgc                                                                                        <  2:1314822/93‑1 (MQ=255)

 cAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAgcgc                                                                                        <  1:372040/93‑1 (MQ=255)

 cAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAgcgc                                                                                        >  1:3607414/1‑93 (MQ=255)

  aGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCt                                                                                       >  1:4191642/1‑93 (MQ=255)

   gTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCtt                                                                                      <  2:1538437/93‑1 (MQ=255)

   gTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCtt                                                                                      >  2:2822734/1‑93 (MQ=255)

   gTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCtt                                                                                      >  2:4138746/1‑93 (MQ=255)

    ttAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCttt                                                                                     <  1:1406882/93‑1 (MQ=255)

    ttAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCttt                                                                                     >  1:967936/1‑93 (MQ=255)

    ttAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGGTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCttt                                                                                     >  2:1684910/1‑93 (MQ=255)

       aaTGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGcc                                                                                  <  2:3413824/93‑1 (MQ=255)

       aaTGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGcc                                                                                  >  1:1544362/1‑93 (MQ=255)

       aaTGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGcc                                                                                  <  2:758647/93‑1 (MQ=255)

       aaTGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCAGGAGTTATCGACTGACCGAAGCCGATTTTCGTGGTGAGCGCTTTGcc                                                                                  <  2:4281154/93‑1 (MQ=255)

         tGAACGTATTCTGGTGCTGGACGGCGTTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGa                                                                                <  1:95698/93‑1 (MQ=255)

         tGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGa                                                                                <  1:1893297/93‑1 (MQ=255)

          gAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGAc                                                                               <  2:1296208/93‑1 (MQ=255)

           aaCGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACt                                                                              >  1:3797431/1‑93 (MQ=255)

              gTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGc                                                                           <  1:1570971/93‑1 (MQ=255)

                aTTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCa                                                                         <  2:1028175/93‑1 (MQ=255)

                aTTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCGTTACCGACTGGCCa                                                                         >  2:2727925/1‑93 (MQ=255)

                 ttCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCAt                                                                        <  2:2125317/93‑1 (MQ=255)

                 ttCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCAt                                                                        >  1:3453946/1‑93 (MQ=255)

                  tCTGGTGCTGGACGGCGTTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATg                                                                       <  2:1362204/93‑1 (MQ=255)

                  tCTGGTGCTGGACGGCGGTATGGGCACCAGGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATg                                                                       <  2:2179847/93‑1 (MQ=255)

                   cTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGc                                                                      >  1:1290569/1‑93 (MQ=255)

                   cTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGc                                                                      <  1:1931509/93‑1 (MQ=255)

                     ggTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGa                                                                    <  2:3043635/93‑1 (MQ=255)

                       tGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGAcc                                                                  >  2:1395828/1‑93 (MQ=255)

                       tGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGAcc                                                                  >  1:2460822/1‑93 (MQ=255)

                        gCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCt                                                                 <  1:4462460/93‑1 (MQ=255)

                        gCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCt                                                                 <  2:4214134/93‑1 (MQ=255)

                            gACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCaaa                                                             <  1:3821375/93‑1 (MQ=255)

                            gACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCaaa                                                             <  1:722392/93‑1 (MQ=255)

                            gACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCaaa                                                             <  2:2606170/93‑1 (MQ=255)

                             aCGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAg                                                            >  1:3932548/1‑93 (MQ=255)

                             aCGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAg                                                            <  2:2038676/93‑1 (MQ=255)

                              cggcggTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGGCTGGCCATGCGACCTCAAAgg                                                           >  2:2270622/1‑93 (MQ=255)

                              cggcggTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAgg                                                           <  1:2791990/93‑1 (MQ=255)

                              cggcggTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAgg                                                           >  1:1659366/1‑93 (MQ=255)

                                gcggTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGca                                                         <  2:3975237/93‑1 (MQ=255)

                                 cggTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGcaa                                                        <  1:2049157/93‑1 (MQ=255)

                                  ggTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGcaac                                                       <  1:2680692/93‑1 (MQ=255)

                                   gTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGcaaca                                                      >  1:471645/1‑93 (MQ=255)

                                    tATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGcaacaa                                                     >  1:3368035/1‑93 (MQ=255)

                                    tATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGcaacaa                                                     <  1:1770410/93‑1 (MQ=255)

                                    tATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGcaacaa                                                     >  1:2807653/1‑93 (MQ=255)

                                    tATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGcaacaa                                                     <  2:426952/93‑1 (MQ=255)

                                    tATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGcaacaa                                                     <  2:35230/93‑1 (MQ=255)

                                      tGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGTCCTCAAAGGCAACAACg                                                   >  2:2323807/1‑93 (MQ=255)

                                        ggCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGAc                                                 <  1:4158567/93‑1 (MQ=255)

                                        ggCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGAc                                                 >  1:3061562/1‑93 (MQ=255)

                                        ggCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGAc                                                 <  2:4155751/93‑1 (MQ=255)

                                        ggCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGAc                                                 <  1:2755537/93‑1 (MQ=255)

                                          cACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACct                                               <  1:2395672/93‑1 (MQ=255)

                                           aCCATGATCAAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGGAAGGGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACctg                                              >  2:4557675/1‑93 (MQ=255)

                                            ccATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACctgc                                             >  2:646829/1‑93 (MQ=255)

                                             cATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACctgct                                            >  2:3847979/1‑93 (MQ=255)

                                              aTGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACctgctg                                           >  2:3400838/1‑93 (MQ=255)

                                              aTGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGGGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACctgctg                                           <  1:644279/93‑1 (MQ=255)

                                                  tCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTAc                                       >  2:335208/1‑93 (MQ=255)

                                                   ccAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACt                                      >  1:562828/1‑93 (MQ=255)

                                                      gagTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAg                                   >  1:708970/1‑93 (MQ=255)

                                                        gTTATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTa                                 >  2:1241632/1‑93 (MQ=255)

                                                         ttATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTaa                                <  1:3215570/93‑1 (MQ=255)

                                                         ttATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTaa                                <  1:4386130/93‑1 (MQ=255)

                                                         ttATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTaa                                <  2:2575574/93‑1 (MQ=255)

                                                          tATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTaaa                               >  2:4487466/1‑93 (MQ=255)

                                                          tATCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTaaa                               <  2:2486552/93‑1 (MQ=255)

                                                            tCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAAcc                             <  1:3259127/93‑1 (MQ=255)

                                                            tCGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACATGCTGGTACTCAGTAAAcc                             >  1:899078/1‑93 (MQ=255)

                                                             cGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCg                            <  2:1932878/93‑1 (MQ=255)

                                                             cGACTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAAACAACAACGACCTGCTGGTACTCAGTAAACCg                            >  2:1194244/1‑93 (MQ=255)

                                                                cTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGaa                         >  1:1045039/1‑93 (MQ=255)

                                                                cTGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGaa                         <  1:2937853/93‑1 (MQ=255)

                                                                 tGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAg                        >  2:540422/1‑93 (MQ=255)

                                                                 tGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAg                        <  1:2260250/93‑1 (MQ=255)

                                                                 tGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAg                        <  2:3141572/93‑1 (MQ=255)

                                                                 tGAACGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAg                        <  1:2350501/93‑1 (MQ=255)

                                                                    aCGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGa                     <  1:4257705/93‑1 (MQ=255)

                                                                    aCGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGa                     >  2:3612208/1‑93 (MQ=255)

                                                                     cGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGAt                    <  1:3493688/93‑1 (MQ=255)

                                                                     cGAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAAAAACGACCTGCTGGTACTCAGTAAACCGGAAGTGag                    >  1:593958/1‑92 (MQ=255)

                                                                      gAAGCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAACGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATc                   <  1:2506223/93‑1 (MQ=255)

                                                                         gCCGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATcgcc                >  1:2120837/1‑93 (MQ=255)

                                                                          ccGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATcgccg               >  2:2215841/1‑93 (MQ=255)

                                                                           cGATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATcgccgc              >  1:4333377/1‑93 (MQ=255)

                                                                           cGATTTTCGTGGTGAGCCCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATcgccgc              <  2:3688300/93‑1 (MQ=255)

                                                                            gATTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCt             >  2:2953281/1‑93 (MQ=255)

                                                                             aTTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTa            >  1:3431789/1‑93 (MQ=255)

                                                                             aTTTTCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTa            <  2:4446702/93‑1 (MQ=255)

                                                                              ttttCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTac           >  1:1007765/1‑92 (MQ=255)

                                                                              ttttCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTAt           >  1:2969509/1‑93 (MQ=255)

                                                                              ttttCGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTAt           >  2:1742403/1‑93 (MQ=255)

                                                                                  cGTGGTGAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCcac       <  1:3293151/93‑1 (MQ=255)

                                                                                   gtggtgAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCcaca      >  1:2617927/1‑93 (MQ=255)

                                                                                    tggtgAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTATACCGGAAGTGATCGCCGCTATCcacac     >  2:2219744/1‑92 (MQ=255)

                                                                                    tggtgAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCCACaa     <  2:835718/93‑1 (MQ=255)

                                                                                    tggtgAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCCACaa     <  2:4090799/93‑1 (MQ=255)

                                                                                       tgAGCGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCCACAACGc  <  1:4064086/93‑1 (MQ=255)

GCAGTTAAATGAACGTATTCTGGTGCTGGACGGCGGTATGGGCACCATGATCCAGAGTTATCGACTGAACGAAGCCGATTTTCGTGGTGAACGCTTTGCCGACTGGCCATGCGACCTCAAAGGCAACAACGACCTGCTGGTACTCAGTAAACCGGAAGTGATCGCCGCTATCCACAACGC  >  NC_000913/4221436‑4221615

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