BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,648,042	T→C	L59P (CTG→CCG)	arsC →	arsenate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,648,042	0	T	C	96.6%	610.1 / 6.6	205	L59P (CTG→CCG)	arsC	arsenate reductase
Reads supporting (aligned to +/- strand): ref base T (2/5); new base C (76/122); total (78/127)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 228 total aligned reads displayed.

CGACTATTATCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCTGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTCAGCAC  >  NC_000913/3647952‑3648130

cGACTATTATCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCtgccgc                                                                                        <  2:659006/93‑1 (MQ=255)

 gACTATTATCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAGACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCcgcg                                                                                       >  2:3096905/1‑93 (MQ=255)

  aCTATTATCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGt                                                                                      <  2:574962/93‑1 (MQ=255)

   cTATTATCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTa                                                                                     >  2:447389/1‑93 (MQ=255)

    tattatCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTaa                                                                                    >  2:3707711/1‑93 (MQ=255)

    tattatCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTaa                                                                                    >  1:4653285/1‑93 (MQ=255)

    tattatCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTaa                                                                                    >  2:2764895/1‑93 (MQ=255)

    tattatCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTaa                                                                                    >  1:476140/1‑93 (MQ=255)

       tatCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTaaaaa                                                                                 <  2:1621591/93‑1 (MQ=255)

       tatCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTaaaaa                                                                                 <  2:542970/93‑1 (MQ=255)

        atCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGTTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAAc                                                                                <  2:812219/93‑1 (MQ=255)

        atCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAAc                                                                                >  2:3976670/1‑93 (MQ=255)

        atCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAAc                                                                                <  2:3541151/93‑1 (MQ=255)

        atCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAAc                                                                                <  1:1751531/93‑1 (MQ=255)

        atCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAAc                                                                                >  2:3326818/1‑93 (MQ=255)

         tCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACg                                                                               >  1:3966597/1‑93 (MQ=255)

           cATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTc                                                                             >  1:2747003/1‑93 (MQ=255)

            aTTATCTGGAAACTCCGCCAACGCGCGATGAACTGTTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCg                                                                            >  2:508852/1‑93 (MQ=255)

            aTTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCg                                                                            <  2:102599/93‑1 (MQ=255)

             ttATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGa                                                                           >  2:4313716/1‑93 (MQ=255)

              tATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGaa                                                                          <  1:3505976/93‑1 (MQ=255)

               aTCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAAc                                                                         >  2:460659/1‑93 (MQ=255)

               aTCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAAc                                                                         <  1:13858/93‑1 (MQ=255)

                tCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAAcc                                                                        <  1:3148490/93‑1 (MQ=255)

                   ggAAACTCCGCTAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTa                                                                     <  2:1238312/93‑1 (MQ=255)

                   ggAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTa                                                                     <  1:2057915/93‑1 (MQ=255)

                   ggAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTa                                                                     <  1:4586177/93‑1 (MQ=255)

                    gAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTAt                                                                    >  2:4498948/1‑93 (MQ=255)

                     aaaCTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATg                                                                   <  1:4713958/93‑1 (MQ=255)

                       aCTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATgag                                                                 <  1:390984/93‑1 (MQ=255)

                         tCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATgagga                                                               >  2:250037/1‑93 (MQ=255)

                         tCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATgagga                                                               >  2:4218266/1‑93 (MQ=255)

                         tCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATgagga                                                               <  1:153062/93‑1 (MQ=255)

                          ccGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATgaggag                                                              <  2:3673036/93‑1 (MQ=255)

                           cGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGc                                                             <  2:1534678/93‑1 (MQ=255)

                            gCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCTGCGTAAAAACGTCGAACCGTATGAGGAGCt                                                            <  1:4279582/93‑1 (MQ=255)

                            gCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCt                                                            <  1:2305338/93‑1 (MQ=255)

                            cccAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCt                                                            <  2:1478472/92‑1 (MQ=255)

                              cAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTgg                                                          >  2:2613270/1‑93 (MQ=255)

                               aaCGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTggg                                                         <  2:68953/93‑1 (MQ=255)

                               aaCGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTggg                                                         >  1:3254147/1‑93 (MQ=255)

                               aaCGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTggc                                                         >  1:1501423/1‑92 (MQ=255)

                                aCGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTgggg                                                        >  2:384744/1‑92 (MQ=255)

                                 cgcgcgTTGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTAAGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGcc                                                       >  2:2650136/1‑93 (MQ=255)

                                 cgcgcgATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGcc                                                       >  1:438896/1‑93 (MQ=255)

                                  gcgcgATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCt                                                      <  2:1429129/93‑1 (MQ=255)

                                   cgcgATGAACTGGTAAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCtt                                                     <  2:3608694/93‑1 (MQ=255)

                                    gcgATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTg                                                    >  1:363712/1‑93 (MQ=255)

                                     cgATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGc                                                   <  1:2175826/93‑1 (MQ=255)

                                     cgATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGc                                                   <  1:1820060/93‑1 (MQ=255)

                                         gAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGaa                                               <  2:289258/93‑1 (MQ=255)

                                            cTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGAt                                            >  1:2358223/1‑93 (MQ=255)

                                            cTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGAt                                            >  2:4636335/1‑93 (MQ=255)

                                             tGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGGATGAGGAGCTGGGCCTTGCGGACGATa                                           >  1:3565194/1‑93 (MQ=255)

                                              ggTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCTCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATaa                                          <  1:4241289/93‑1 (MQ=255)

                                              ggTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATaa                                          <  1:929384/93‑1 (MQ=255)

                                              ggTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATaa                                          <  1:3271190/93‑1 (MQ=255)

                                              ggTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATaa                                          <  2:2471154/93‑1 (MQ=255)

                                                tCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAAt                                        <  2:733812/93‑1 (MQ=255)

                                                 cAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAAtt                                       <  2:892222/93‑1 (MQ=255)

                                                 cAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAAtt                                       <  2:2088237/93‑1 (MQ=255)

                                                   aaCTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTa                                     >  1:27017/1‑93 (MQ=255)

                                                   aaCTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTa                                     <  1:3767533/93‑1 (MQ=255)

                                                    aCTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTAc                                    >  2:338443/1‑93 (MQ=255)

                                                     cTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACt                                   <  2:718338/93‑1 (MQ=255)

                                                     cTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACt                                   <  2:3792892/93‑1 (MQ=255)

                                                     cTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACt                                   >  2:2916109/1‑93 (MQ=255)

                                                     cTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACt                                   >  1:4513935/1‑93 (MQ=255)

                                                      ccATTGCCGATATGGGGA‑TTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGa                                 <  1:736696/92‑1 (MQ=255)

                                                       cATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGa                                 <  1:291973/93‑1 (MQ=255)

                                                       cATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGa                                 <  1:2979290/93‑1 (MQ=255)

                                                       cATTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGa                                 <  1:4136524/93‑1 (MQ=255)

                                                        aTTGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGAc                                <  2:2241370/93‑1 (MQ=255)

                                                         ttGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACg                               <  1:1110189/93‑1 (MQ=255)

                                                          tGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGa                              <  2:1975313/93‑1 (MQ=255)

                                                          tGCCGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACATCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGa                              <  2:2999558/93‑1 (MQ=255)

                                                            ccGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATc                            <  2:1788790/93‑1 (MQ=255)

                                                            ccGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATc                            <  2:3889118/93‑1 (MQ=255)

                                                            ccGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATc                            >  1:3942674/1‑93 (MQ=255)

                                                             cGATATGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCg                           <  1:523701/93‑1 (MQ=255)

                                                               atatGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGt                         <  1:4436652/93‑1 (MQ=255)

                                                                tatGGGGATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGtt                        >  2:32290/1‑93 (MQ=255)

                                                                 atGGGGATTTCCGTACGCGCGCTGCTGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTa                       >  2:694168/1‑93 (MQ=255)

                                                                    gggATTTCCGTACGCGCGCTGCTGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATc                    <  2:4367141/93‑1 (MQ=255)

                                                                    gggATTTCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATc                    <  2:1046814/93‑1 (MQ=255)

                                                                         ttCCGTACGCGCGCTGCCGCGTAAAAGCGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACtt               <  2:233069/93‑1 (MQ=255)

                                                                         ttCCGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACtt               >  1:817820/1‑93 (MQ=255)

                                                                            cGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTAt            <  1:263932/93‑1 (MQ=255)

                                                                            cGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTAt            >  1:2585453/1‑93 (MQ=255)

                                                                            cGTACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTAt            <  1:4675113/93‑1 (MQ=255)

                                                                              tACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGc          >  2:1303105/1‑93 (MQ=255)

                                                                              tACGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGc          <  1:1236392/93‑1 (MQ=255)

                                                                               aCGCGCGCTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCt         >  1:4543432/1‑93 (MQ=255)

                                                                                 gcgcgcTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTc       <  1:2990475/93‑1 (MQ=255)

                                                                                 gcgcgcTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTc       <  1:86096/93‑1 (MQ=255)

                                                                                 gcgcgcTGCCGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTc       <  1:1205724/93‑1 (MQ=255)

                                                                                     gctgccGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTCAGCa   <  1:104471/93‑1 (MQ=255)

                                                                                     gctgccGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTCAGCa   <  1:1572235/93‑1 (MQ=255)

                                                                                      ctgccgCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTCAGCAc  >  1:2284210/1‑93 (MQ=255)

                                                                                      ctgccgCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTCAGCAc  <  1:43351/93‑1 (MQ=255)

CGACTATTATCCATTATCTGGAAACTCCGCCAACGCGCGATGAACTGGTCAAACTCATTGCCGATATGGGGATTTCCGTACGCGCGCTGCTGCGTAAAAACGTCGAACCGTATGAGGAGCTGGGCCTTGCGGAAGATAAATTTACTGACGATCGGTTAATCGACTTTATGCTTCAGCAC  >  NC_000913/3647952‑3648130

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