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	99.5%	745.5 / NA	214	G51G (GGT→GGC)	sbm	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (104/108); total (105/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

TGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTT  >  NC_000913/3058932‑3059111

  tCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCTATCTCGATAATCTGGAGGTGACAGGCAccc                                                                                       <  1:1881877/93‑1 (MQ=255)

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

  gcGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCAccc                                                                                       <  2:662180/92‑1 (MQ=255)

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

     aCTCGCTGGTTCATCAAACCGCGGAAGGTATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTc                                                                                    >  2:4235117/1‑93 (MQ=255)

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

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

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

       tCGCTGGTTCACCAAACCGCTGAAGGGATCCCCATCAAGCCGCTTTATACCGAAGCCCATCTCGATAATCTGGAGGTGACAGGCACCCTTCCt                                                                                  <  1:1735925/93‑1 (MQ=255)

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

           tGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGtt                                                                              >  1:651778/1‑93 (MQ=255)

            ggTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGttt                                                                             >  2:636317/1‑93 (MQ=255)

             gTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTg                                                                            >  1:3184399/1‑93 (MQ=255)

              ttCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTgc                                                                           <  2:1235591/93‑1 (MQ=255)

              taaatcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTgc                                                                           <  1:5817117/90‑1 (MQ=255)

               tcatcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTgcc                                                                          <  2:3844367/93‑1 (MQ=255)

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

                 atcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTgccgc                                                                        <  1:2563201/93‑1 (MQ=255)

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

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

                     aaCCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTa                                                                    <  2:5800883/93‑1 (MQ=255)

                      aCCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTAc                                                                   <  2:5162822/93‑1 (MQ=255)

                      aCCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTAc                                                                   <  1:4332239/93‑1 (MQ=255)

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

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

                        cgcgGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGt                                                                 <  2:6126794/93‑1 (MQ=255)

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

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

                           ggAGGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCg                                                              <  1:5398757/93‑1 (MQ=255)

                           ggAAGGGATCGCCATCAAGCCGCTGTATGCCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCg                                                              <  2:4121555/93‑1 (MQ=255)

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

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

                            gAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGt                                                             >  2:2284246/1‑93 (MQ=255)

                            gAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGt                                                             <  1:1311204/93‑1 (MQ=255)

                            gAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGt                                                             >  2:3016560/1‑93 (MQ=255)

                            gAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGt                                                             >  1:4166052/1‑93 (MQ=255)

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

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

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

                                ggATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTggtc                                                         >  1:228021/1‑91 (MQ=255)

                                ggATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGcc                                                         >  2:1100635/1‑93 (MQ=255)

                                 gATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGccc                                                        >  1:761513/1‑93 (MQ=255)

                                 gATCGCCATCAAGCCGCTGTATACCAAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGccc                                                        <  1:4826117/93‑1 (MQ=255)

                                   tCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCcgc                                                      <  1:3195207/93‑1 (MQ=255)

                                    cGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCcgcg                                                     >  1:3749872/1‑93 (MQ=255)

                                     gCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCcgcgg                                                    >  2:1677176/1‑92 (MQ=255)

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

                                       cATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGc                                                  <  2:3625042/93‑1 (MQ=255)

                                        aTCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGcc                                                 >  2:3040391/1‑93 (MQ=255)

                                        aTCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGcc                                                 <  2:4550379/93‑1 (MQ=255)

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

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

                                           aaGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGAACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCCTGCCACt                                              >  1:2805039/1‑93 (MQ=255)

                                            aGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTa                                             <  1:5876174/93‑1 (MQ=255)

                                              ccGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATg                                           <  1:4703486/93‑1 (MQ=255)

                                               cGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGt                                          <  2:3951987/93‑1 (MQ=255)

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

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

                                                gCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGta                                         >  2:1546828/1‑93 (MQ=255)

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

                                                  tGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACtatatata                                       >  1:5535195/1‑93 (MQ=255)

                                                  tGTATACCGAAGCCGATCTAGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGtata                                       >  1:5448517/1‑93 (MQ=255)

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

                                                     ataCCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTCCCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCg                                    <  1:2191320/93‑1 (MQ=255)

                                                     ataCCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCg                                    >  1:3839491/1‑93 (MQ=255)

                                                       aCCGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGcc                                  >  2:5617267/1‑93 (MQ=255)

                                                        ccGAAGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGccc                                 >  2:254459/1‑93 (MQ=255)

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

                                                            aGCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAAcc                             <  2:3988824/93‑1 (MQ=255)

                                                             gCCGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCg                            >  2:3483489/1‑93 (MQ=255)

                                                             accGATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCg                            <  1:5179086/92‑1 (MQ=255)

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

                                                                gATCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTgg                         >  1:5804201/1‑93 (MQ=255)

                                                                 aTCTCGATAATCTGGAGGTGACAGGCACCCTTCCTGGGTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGa                        >  1:716947/1‑93 (MQ=255)

                                                                    tcGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCa                     <  1:4030022/93‑1 (MQ=255)

                                                                     cGATAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACGGCGGACCAt                    >  2:680477/1‑93 (MQ=255)

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

                                                                      cagcaTCTGGAGGAGACGGGCACCCTTCCTGGTTTTCCGCCCTGCGTACGTGGCCCGCGGGCCACTATGTATACCGCCCAACCGGGGACCATc                   <  1:5086363/89‑1 (MQ=255)

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

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

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

                                                                        tAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCg                 >  2:2381080/1‑93 (MQ=255)

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

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

                                                                        tAATCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCg                 >  2:3187809/1‑93 (MQ=255)

                                                                        tAATCGGGAGGTGACAGGCCCCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCg                 <  1:1270686/93‑1 (MQ=255)

                                                                         aaTCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCgc                >  2:4303350/1‑92 (MQ=255)

                                                                         aaTCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGt                >  1:2692521/1‑93 (MQ=255)

                                                                          aTCTGGAGGTGACAGGCACCCTTCCTGGTTTGACGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATATCGTCCAACCGCGGAACAACCGTc               >  1:3890707/1‑93 (MQ=255)

                                                                           tCTGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCa              <  1:2087549/93‑1 (MQ=255)

                                                                           ccTGGAGGGGACAGGCACCCTTCCTGACTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCa              <  1:1372957/92‑1 (MQ=255)

                                                                             tGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGt            >  2:5024163/1‑93 (MQ=255)

                                                                             tGGAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGt            <  2:118284/93‑1 (MQ=255)

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

                                                                               gAGGTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTAt          >  1:4217970/1‑93 (MQ=255)

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

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

                                                                                  gTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATgat       >  1:2162848/1‑91 (MQ=255)

                                                                                  gTGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCt       >  2:1605234/1‑93 (MQ=255)

                                                                                   tGACAGGCACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCACGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTg      <  1:4492287/93‑1 (MQ=255)

TGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTT  >  NC_000913/3058932‑3059111

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