BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,059,022	T→C	G51G (GGT→GGC)	sbm →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,059,022	0	T	C	100.0%	786.5 / NA	219	G51G (GGT→GGC)	sbm	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (115/104); total (115/104)

Only 100 of 221 total aligned reads displayed.

TGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTT  >  NC_000913/3058932‑3059111

tGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGTAGGTGACAGGCAc                                                                                         >  2:4523197/1‑93 (MQ=255)

tGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCAc                                                                                         <  2:4878735/93‑1 (MQ=255)

  tCGCCTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCAccc                                                                                       <  1:4594079/93‑1 (MQ=255)

  tCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCAccc                                                                                       <  1:3697367/93‑1 (MQ=255)

  tCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCAccc                                                                                       <  2:15731/93‑1 (MQ=255)

   cGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCt                                                                                      >  1:899113/1‑93 (MQ=255)

    gACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCtt                                                                                     >  2:2584909/1‑93 (MQ=255)

     aCTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTc                                                                                    >  2:1299408/1‑93 (MQ=255)

     aCTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTc                                                                                    <  2:2847181/93‑1 (MQ=255)

     aCTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTc                                                                                    >  1:2599748/1‑93 (MQ=255)

      cTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTcc                                                                                   <  2:5717443/93‑1 (MQ=255)

       tCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCt                                                                                  <  2:4006829/93‑1 (MQ=255)

       tCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCt                                                                                  >  1:1731529/1‑93 (MQ=255)

       tCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCt                                                                                  <  1:1408915/93‑1 (MQ=255)

       tCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCt                                                                                  >  1:1303130/1‑93 (MQ=255)

       tCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCt                                                                                  >  2:6154565/1‑93 (MQ=255)

        cGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTg                                                                                 <  2:5518929/93‑1 (MQ=255)

         gCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTgg                                                                                >  1:1153284/1‑93 (MQ=255)

          cTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGt                                                                               >  2:2298471/1‑93 (MQ=255)

            ggTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTACCTGGttt                                                                             <  1:6502682/93‑1 (MQ=255)

                catcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTgccg                                                                         >  1:3329475/1‑93 (MQ=255)

                 atcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTgccgc                                                                        >  2:1671883/1‑93 (MQ=255)

                  tcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTgccgcc                                                                       >  2:5838746/1‑93 (MQ=255)

                  tcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTgccgcc                                                                       >  2:5243110/1‑93 (MQ=255)

                  tcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTgccgcc                                                                       <  2:2709643/93‑1 (MQ=255)

                   caAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGccc                                                                      <  1:2286111/93‑1 (MQ=255)

                    aaaCCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCt                                                                     <  1:5292319/93‑1 (MQ=255)

                    aaaCCGCGGAAGGGATCGCCAACAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCt                                                                     <  1:5943805/93‑1 (MQ=255)

                       ccGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACg                                                                  <  1:4381098/93‑1 (MQ=255)

                         gcgGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGtt                                                                <  2:3798729/93‑1 (MQ=255)

                          cgGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTc                                                               <  2:631273/93‑1 (MQ=255)

                           ggAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCg                                                              <  1:5513172/93‑1 (MQ=255)

                           ggAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCg                                                              >  1:4409103/1‑93 (MQ=255)

                             aaGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTg                                                            <  2:8543/93‑1 (MQ=255)

                              aGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTgg                                                           >  1:2984210/1‑93 (MQ=255)

                              aGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTgg                                                           <  2:2231420/93‑1 (MQ=255)

                               gggATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGc                                                          >  1:6131314/1‑93 (MQ=255)

                                ggATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGcc                                                         <  1:6079956/93‑1 (MQ=255)

                                  aTCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCcg                                                       <  2:4231249/93‑1 (MQ=255)

                                  aTCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCcg                                                       <  1:1854540/93‑1 (MQ=255)

                                   tCGCCGTCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCcgc                                                      >  2:2429145/1‑93 (MQ=255)

                                   tCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACTTTCGTGGCCcgc                                                      <  2:3585918/93‑1 (MQ=255)

                                     gCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGt                                                    >  1:2866608/1‑93 (MQ=255)

                                     gCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGt                                                    >  2:1544947/1‑93 (MQ=255)

                                      ccATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTg                                                   >  2:3110360/1‑93 (MQ=255)

                                       cATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGc                                                  >  2:1213192/1‑93 (MQ=255)

                                         tCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCa                                                <  2:2060731/93‑1 (MQ=255)

                                         tCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCa                                                >  2:3935183/1‑93 (MQ=255)

                                          cAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCAc                                               <  1:1061697/93‑1 (MQ=255)

                                          cAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCAc                                               <  1:5582215/93‑1 (MQ=255)

                                          cAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCAc                                               >  1:6382494/1‑93 (MQ=255)

                                             gCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTAt                                            >  1:5789083/1‑93 (MQ=255)

                                               cGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGt                                          >  1:4785188/1‑93 (MQ=255)

                                               cGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGt                                          >  2:3476578/1‑93 (MQ=255)

                                                gCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGta                                         <  2:5661583/93‑1 (MQ=255)

                                                 cTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGtat                                        <  2:1263756/93‑1 (MQ=255)

                                                 cTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGtat                                        >  2:2360746/1‑93 (MQ=255)

                                                  tGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGtata                                       <  1:5137434/93‑1 (MQ=255)

                                                   gTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATAc                                      <  1:3903260/93‑1 (MQ=255)

                                                   gTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATAc                                      >  1:5453088/1‑93 (MQ=255)

                                                    tataCCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATAcc                                     >  2:3240551/1‑93 (MQ=255)

                                                      taCCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGc                                   <  1:4142181/93‑1 (MQ=255)

                                                       aCCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGcc                                  <  2:623900/93‑1 (MQ=255)

                                                        ccGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGccc                                 <  1:1791739/93‑1 (MQ=255)

                                                          gAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCaa                               <  1:2611527/93‑1 (MQ=255)

                                                           aaGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAAc                              >  2:5800545/1‑93 (MQ=255)

                                                           aaGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAAc                              <  1:3777988/93‑1 (MQ=255)

                                                            aGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAAcc                             >  1:2496541/1‑93 (MQ=255)

                                                              ccGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGt                           <  1:6250872/93‑1 (MQ=255)

                                                               cGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTg                          <  2:436263/93‑1 (MQ=255)

                                                                gATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTgg                         <  2:4603023/93‑1 (MQ=255)

                                                                 aTCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGa                        >  1:168914/1‑93 (MQ=255)

                                                                 aTCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGa                        <  2:5360155/93‑1 (MQ=255)

                                                                 aTCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGa                        <  2:1779203/93‑1 (MQ=255)

                                                                 aTCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGa                        >  2:4730453/1‑93 (MQ=255)

                                                                 aTCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGa                        >  2:427833/1‑93 (MQ=255)

                                                                  tctcGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGAc                       >  1:360540/1‑93 (MQ=255)

                                                                    tcGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGTCCCCGTGCCACTATGTATACCGCCCAACCGTGGAcct                     >  1:592862/1‑92 (MQ=255)

                                                                    tcGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCa                     >  1:1194369/1‑93 (MQ=255)

                                                                    tcGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCa                     >  1:2673245/1‑93 (MQ=255)

                                                                    tcGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATAACGCCCAACCGTGGACCa                     >  2:5964380/1‑93 (MQ=255)

                                                                      gATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATc                   >  1:430998/1‑93 (MQ=255)

                                                                      gATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATc                   <  1:5716127/93‑1 (MQ=255)

                                                                      gATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATc                   <  2:6051463/93‑1 (MQ=255)

                                                                      gATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATc                   >  2:5342029/1‑93 (MQ=255)

                                                                       aTAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATcc                  >  2:3331671/1‑93 (MQ=255)

                                                                       aTAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATcc                  <  1:4546727/93‑1 (MQ=255)

                                                                        tAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCg                 <  2:4452145/93‑1 (MQ=255)

                                                                          aTCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTc               <  2:209023/93‑1 (MQ=255)

                                                                           tCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCa              >  2:746175/1‑93 (MQ=255)

                                                                           tCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCa              >  1:701744/1‑93 (MQ=255)

                                                                              ggAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTa           >  2:4833359/1‑93 (MQ=255)

                                                                               gAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTAt          <  1:659309/93‑1 (MQ=255)

                                                                                 ggTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGc        >  2:1378389/1‑93 (MQ=255)

                                                                                 ggTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGc        <  2:1492246/93‑1 (MQ=255)

                                                                                   tGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTg      <  2:3063258/93‑1 (MQ=255)

                                                                                    gACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTgg     <  1:3089028/93‑1 (MQ=255)

                                                                                      cAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGtt   >  2:6601342/1‑93 (MQ=255)

                                                                                      cAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGtt   >  2:101955/1‑93 (MQ=255)

                                                                                      cAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGtt   <  2:925239/93‑1 (MQ=255)

TGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTT  >  NC_000913/3058932‑3059111

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