BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,059,787	A→G	A306A (GCA→GCG)	sbm →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,059,787	0	A	G	100.0%	711.2 / NA	219	A306A (GCA→GCG)	sbm	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (111/108); total (111/108)

Only 100 of 225 total aligned reads displayed.

TTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCG  >  NC_000913/3059698‑3059876

ttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAg                                                                                        >  1:4642276/1‑93 (MQ=255)

  cggcatcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGa                                                                                      <  1:2390352/93‑1 (MQ=255)

  cggcatcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGAGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGa                                                                                      <  1:905847/93‑1 (MQ=255)

   ggcatcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGAc                                                                                     <  2:4934692/93‑1 (MQ=255)

    gcatcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGAcc                                                                                    <  2:5408289/93‑1 (MQ=255)

     catcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGAccc                                                                                   <  1:4792066/93‑1 (MQ=255)

     catcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGAccc                                                                                   >  1:4657884/1‑93 (MQ=255)

      atcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCg                                                                                  >  1:5592204/1‑93 (MQ=255)

      atcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCg                                                                                  >  1:1148286/1‑93 (MQ=255)

      atcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCg                                                                                  >  1:242985/1‑93 (MQ=255)

       tcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGa                                                                                 >  1:5400037/1‑93 (MQ=255)

       tcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGa                                                                                 >  2:2795906/1‑93 (MQ=255)

         ggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGaaa                                                                               >  1:3675833/1‑93 (MQ=255)

          gcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAAt                                                                              >  2:894716/1‑93 (MQ=255)

          gcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAAt                                                                              >  2:2281954/1‑93 (MQ=255)

            atGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCa                                                                            >  1:3296775/1‑93 (MQ=255)

            atGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCa                                                                            >  2:5529826/1‑93 (MQ=255)

             tGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCAc                                                                           >  1:4462714/1‑93 (MQ=255)

              ggATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACt                                                                          <  2:4313744/93‑1 (MQ=255)

              ggATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACt                                                                          <  2:565397/93‑1 (MQ=255)

              ggATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACt                                                                          <  2:3632025/93‑1 (MQ=255)

               gATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTg                                                                         <  1:3099753/93‑1 (MQ=255)

               gATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTg                                                                         >  1:2005832/1‑93 (MQ=255)

                aTCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTgg                                                                        >  2:5432756/1‑93 (MQ=255)

                aTCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTgg                                                                        >  2:61054/1‑93 (MQ=255)

                 tCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGgc                                                                       >  2:4288898/1‑93 (MQ=255)

                 tCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGgc                                                                       >  2:6235709/1‑93 (MQ=255)

                  cTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGgcg                                                                      >  2:5358213/1‑93 (MQ=255)

                   tGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGgcgc                                                                     <  2:2806415/93‑1 (MQ=255)

                     tttATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTg                                                                   >  2:5768219/1‑93 (MQ=255)

                       tATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCg                                                                 <  1:3501451/93‑1 (MQ=255)

                         tGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTa                                                               >  2:3627233/1‑93 (MQ=255)

                          gAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTAc                                                              >  2:5606837/1‑93 (MQ=255)

                          gAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTAc                                                              <  1:1801297/93‑1 (MQ=255)

                           aaCGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTAcc                                                             <  2:4229861/93‑1 (MQ=255)

                            aCGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTAccc                                                            >  1:3792165/1‑93 (MQ=255)

                             cGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCa                                                           <  2:3121433/93‑1 (MQ=255)

                             cGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCa                                                           <  1:4989622/93‑1 (MQ=255)

                             cGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCa                                                           >  1:6116160/1‑93 (MQ=255)

                              gTCGCCATGTTGCGTGCGGCGCGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCAc                                                          <  2:1609734/93‑1 (MQ=255)

                                  ccATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGcc                                                      >  1:1092704/1‑93 (MQ=255)

                                   cATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCa                                                     <  2:5108144/93‑1 (MQ=255)

                                   cATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAATGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCa                                                     <  1:6258125/93‑1 (MQ=255)

                                      gTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGAc                                                  >  2:4692807/1‑93 (MQ=255)

                                       ttGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGAcc                                                 >  1:2054028/1‑93 (MQ=255)

                                       ttGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGAcc                                                 <  2:163677/93‑1 (MQ=255)

                                        tgcgtgcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCt                                                >  1:4162842/1‑93 (MQ=255)

                                         gcgtgcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTc                                               <  1:2160431/93‑1 (MQ=255)

                                          cgtgcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCa                                              <  1:4665151/93‑1 (MQ=255)

                                          cgtgcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCa                                              <  2:1725936/93‑1 (MQ=255)

                                             gcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGc                                           <  2:6250608/93‑1 (MQ=255)

                                             gcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGc                                           <  2:1193429/93‑1 (MQ=255)

                                             gcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGc                                           <  2:2954419/93‑1 (MQ=255)

                                              cgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGTCGCTGCGTACCCACTGCCAGACCTCAGGCt                                          <  2:5254796/93‑1 (MQ=255)

                                              cgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGTCCTCAGGCt                                          <  1:5976655/93‑1 (MQ=255)

                                               ggCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTg                                         >  1:5655473/1‑93 (MQ=255)

                                               ggCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTg                                         >  2:3180740/1‑93 (MQ=255)

                                                gCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTgg                                        >  2:106298/1‑93 (MQ=255)

                                                 cACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGa                                       >  2:2137454/1‑93 (MQ=255)

                                                   cGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGc                                     >  1:741455/1‑93 (MQ=255)

                                                   cGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGc                                     >  1:455885/1‑93 (MQ=255)

                                                      tatttatGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctg                                  >  1:3202992/1‑93 (MQ=255)

                                                       atttatGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctga                                 <  2:1959682/93‑1 (MQ=255)

                                                       atttatGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctga                                 <  2:3107082/93‑1 (MQ=255)

                                                        tttatGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctgac                                >  1:1232175/1‑93 (MQ=255)

                                                         ttatGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctgact                               <  1:1040931/93‑1 (MQ=255)

                                                         ttatGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctgact                               >  1:2206610/1‑93 (MQ=255)

                                                          tatGGAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctgactg                              >  1:3423124/1‑93 (MQ=255)

                                                             ggAGCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAAc                           <  2:4869932/93‑1 (MQ=255)

                                                                gCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAgg                        >  2:2528408/1‑93 (MQ=255)

                                                                gCGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAgg                        <  1:1499993/93‑1 (MQ=255)

                                                                 cGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGa                       <  2:3279894/93‑1 (MQ=255)

                                                                 cGAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGa                       <  1:1832898/93‑1 (MQ=255)

                                                                  gAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGAt                      <  2:2433756/93‑1 (MQ=255)

                                                                  gAAGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGAt                      >  2:1837499/1‑93 (MQ=255)

                                                                   aaGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATc                     >  2:5870249/1‑93 (MQ=255)

                                                                   aaGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATc                     <  1:3950055/93‑1 (MQ=255)

                                                                   aaGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATc                     <  2:3567526/93‑1 (MQ=255)

                                                                   aaGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATc                     <  1:5520314/93‑1 (MQ=255)

                                                                   aaGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATc                     >  2:4456039/1‑93 (MQ=255)

                                                                    aGCGGTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATcc                    >  1:1975517/1‑93 (MQ=255)

                                                                       ggTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGta                 >  1:1568498/1‑93 (MQ=255)

                                                                       ggTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGta                 >  1:4679857/1‑93 (MQ=255)

                                                                        gTCAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGtat                >  1:6295930/1‑93 (MQ=255)

                                                                         gcAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGtata               <  2:1291727/92‑1 (MQ=255)

                                                                          cAGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaa              <  1:519573/93‑1 (MQ=255)

                                                                           aGTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaac             >  2:1252211/1‑93 (MQ=255)

                                                                            gTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaaca            >  2:4539928/1‑93 (MQ=255)

                                                                            gTGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaaca            >  2:4110682/1‑93 (MQ=255)

                                                                             tGGATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaacaa           <  1:5234349/93‑1 (MQ=255)

                                                                             tGGATTGGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTATCCACTGCCAGACCTCAGTCTGGAGCCTGACTGAACAGGATCCGTATaacaa           <  2:2594155/93‑1 (MQ=255)

                                                                               gATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACg         >  1:5775849/1‑93 (MQ=255)

                                                                               gATTTGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACg         >  2:2199577/1‑93 (MQ=255)

                                                                                  ttGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTa      <  2:3918072/93‑1 (MQ=255)

                                                                                  ttGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTa      <  1:388791/93‑1 (MQ=255)

                                                                                   tGGCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTAt     >  1:2781396/1‑93 (MQ=255)

                                                                                    ggCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATc    <  1:2862129/93‑1 (MQ=255)

                                                                                    ggCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATc    <  2:447269/93‑1 (MQ=255)

                                                                                    ggCGCGCAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATc    <  1:3576447/93‑1 (MQ=255)

                                                                                     gcgcgcAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATcc   <  1:4377637/93‑1 (MQ=255)

TTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCG  >  NC_000913/3059698‑3059876

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