BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,982,242	T→C	N460S (AAC→AGC)	aslA ←	arylsulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,982,242	0	T	C	87.3%	810.3 / 90.9	292	N460S (AAC→AGC)	aslA	arylsulfatase
Reads supporting (aligned to +/- strand): ref base T (19/18); new base C (119/136); total (138/154)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 298 total aligned reads displayed.

ACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTT  >  NC_000913/3982151‑3982331

aCCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGtt                                                                                          >  1:5270218/1‑93 (MQ=255)

 cccTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTg                                                                                         >  2:889223/1‑93 (MQ=255)

  ccTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGa                                                                                        <  2:2025965/93‑1 (MQ=255)

  ccTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGa                                                                                        <  1:5387442/93‑1 (MQ=255)

  catGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGACTCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGa                                                                                        <  2:5941914/91‑1 (MQ=255)

   cTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAg                                                                                       >  1:1430414/1‑93 (MQ=255)

     gATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGa                                                                                     >  2:2968607/1‑93 (MQ=255)

      atatCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGaa                                                                                    >  1:1625632/1‑93 (MQ=255)

      atatCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGaa                                                                                    <  2:2455589/93‑1 (MQ=255)

      atatCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGaa                                                                                    <  2:2732409/93‑1 (MQ=255)

       tatCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAag                                                                                   >  2:5500325/1‑93 (MQ=255)

       tatCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAag                                                                                   <  2:2106230/93‑1 (MQ=255)

       tatCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAag                                                                                   >  1:232879/1‑93 (MQ=255)

        atCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAagt                                                                                  >  1:2496465/1‑93 (MQ=255)

        atCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAagt                                                                                  <  2:6554864/93‑1 (MQ=255)

        atCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAagt                                                                                  <  2:260701/93‑1 (MQ=255)

         tCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAagta                                                                                 >  2:4392979/1‑93 (MQ=255)

          ccGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAagtag                                                                                <  2:808192/93‑1 (MQ=255)

           cGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAagtagt                                                                               >  2:5263102/1‑93 (MQ=255)

           cGCTCTGGGTATAAGCGTAAGGTTACTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAagtagt                                                                               <  1:4206825/93‑1 (MQ=255)

             ctctGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGc                                                                             >  2:3347464/1‑93 (MQ=255)

              tctGGGTATAAGCGTAATGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCt                                                                            >  1:3892522/1‑93 (MQ=255)

              tctGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCt                                                                            <  2:4353562/93‑1 (MQ=255)

               ctGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTc                                                                           <  1:3107411/93‑1 (MQ=255)

               ctGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTc                                                                           <  1:5754979/93‑1 (MQ=255)

               ctGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGAGCTc                                                                           >  1:4132034/1‑93 (MQ=255)

                tGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCg                                                                          >  1:3607420/1‑93 (MQ=255)

                tGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCg                                                                          >  1:4609067/1‑93 (MQ=255)

                tGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCg                                                                          >  1:4698671/1‑93 (MQ=255)

                tGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCg                                                                          <  1:1507997/93‑1 (MQ=255)

                tGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCg                                                                          >  1:4546812/1‑93 (MQ=255)

                 gggTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCgg                                                                         <  1:3010500/93‑1 (MQ=255)

                   gTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGcc                                                                       <  2:6226557/93‑1 (MQ=255)

                    tataAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCt                                                                      <  2:1923063/93‑1 (MQ=255)

                    tataAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCt                                                                      >  1:3915116/1‑93 (MQ=255)

                    tataAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCt                                                                      >  2:4169088/1‑93 (MQ=255)

                      taAGCGTAAGTTTGCTAAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTa                                                                    <  1:4792745/93‑1 (MQ=255)

                      taAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTa                                                                    >  2:1201157/1‑93 (MQ=255)

                       aaGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTAc                                                                   >  1:6258918/1‑93 (MQ=255)

                         gCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACgg                                                                 >  1:1703056/1‑93 (MQ=255)

                           gTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGtt                                                               <  1:5096983/93‑1 (MQ=255)

                            tAAGGTTGCTGAATCAGGACGTTATACGTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTa                                                              <  2:6290409/93‑1 (MQ=255)

                            tAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTa                                                              <  1:4343027/93‑1 (MQ=255)

                              aGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGa                                                            >  1:1797449/1‑93 (MQ=255)

                               ggTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGAc                                                           >  2:5167844/1‑93 (MQ=255)

                               ggTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGAc                                                           <  1:2714628/93‑1 (MQ=255)

                                gTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACt                                                          >  2:5445167/1‑93 (MQ=255)

                                 ttGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTg                                                         >  1:3317048/1‑93 (MQ=255)

                                  tGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGa                                                        <  2:69625/93‑1 (MQ=255)

                                   gCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGAc                                                       <  2:4805968/93‑1 (MQ=255)

                                    cTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGAcc                                                      >  2:5226643/1‑93 (MQ=255)

                                    cTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGAcc                                                      <  1:6483303/93‑1 (MQ=255)

                                    cTGAATCAGGACGTGATACTGGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGAcc                                                      <  1:1157601/93‑1 (MQ=255)

                                     tGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCa                                                     <  2:3457935/93‑1 (MQ=255)

                                      gAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCAt                                                    <  2:1402552/93‑1 (MQ=255)

                                      gAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCAt                                                    >  2:1107351/1‑93 (MQ=255)

                                       aaTCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCAtt                                                   <  1:1889434/93‑1 (MQ=255)

                                       aaTCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCAtt                                                   >  1:5925392/1‑93 (MQ=255)

                                          cAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGt                                                <  2:2782154/93‑1 (MQ=255)

                                           aGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGtt                                               >  2:39139/1‑93 (MQ=255)

                                               cGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCa                                           <  1:71904/93‑1 (MQ=255)

                                                 tGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAgg                                         <  1:4833240/93‑1 (MQ=255)

                                                 tGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAgg                                         >  1:2086379/1‑93 (MQ=255)

                                                  gATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGga                                        >  2:5638540/1‑93 (MQ=255)

                                                   aTACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGgaa                                       <  2:3622937/93‑1 (MQ=255)

                                                    tACTTGACCTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGgaag                                      >  2:1885462/1‑93 (MQ=255)

                                                    tACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGgaag                                      >  2:1648497/1‑93 (MQ=255)

                                                    tACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGgaag                                      <  1:68672/93‑1 (MQ=255)

                                                      cTTGACCTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGgaagaa                                    <  1:805171/93‑1 (MQ=255)

                                                      cTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGgaagaa                                    >  2:5820919/1‑93 (MQ=255)

                                                      cTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGgaagaa                                    <  1:3822671/93‑1 (MQ=255)

                                                       ttGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGaagaag                                   >  1:2743699/1‑93 (MQ=255)

                                                       ttGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGaagaag                                   >  1:1017015/1‑93 (MQ=255)

                                                        tGAACTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAgg                                  <  2:2671239/93‑1 (MQ=255)

                                                         gAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGa                                 <  1:4045407/93‑1 (MQ=255)

                                                          aaCTCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGAt                                <  1:1564321/93‑1 (MQ=255)

                                                             tCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTc                             <  2:1819644/93‑1 (MQ=255)

                                                             tCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTc                             >  1:1077766/1‑93 (MQ=255)

                                                             tCATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTc                             >  1:5044761/1‑93 (MQ=255)

                                                              cATCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCt                            <  2:2183371/93‑1 (MQ=255)

                                                                tCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTgg                          >  2:4559615/1‑93 (MQ=255)

                                                                tCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTgg                          <  1:3695714/93‑1 (MQ=255)

                                                                tCCATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTgg                          >  2:6388173/1‑93 (MQ=255)

                                                                  cATACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTc                        >  1:2437852/1‑93 (MQ=255)

                                                                   aTACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTcc                       >  2:3226897/1‑93 (MQ=255)

                                                                    tACGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCca                      >  2:6005591/1‑93 (MQ=255)

                                                                      cGCACAGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCcaca                    >  2:1847632/1‑93 (MQ=255)

                                                                         acaGCAGCGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCa                 >  1:550845/1‑93 (MQ=255)

                                                                             cagcGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGa             <  2:4007386/93‑1 (MQ=255)

                                                                              agcGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGAt            <  1:2952507/93‑1 (MQ=255)

                                                                               gcGAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATa           >  1:6171949/1‑93 (MQ=255)

                                                                                 gAGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATaaa         >  1:3456384/1‑93 (MQ=255)

                                                                                  aGTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAg        <  1:5657274/93‑1 (MQ=255)

                                                                                   gTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAgg       <  2:4609614/93‑1 (MQ=255)

                                                                                   gTTTACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAgg       >  2:3137583/1‑93 (MQ=255)

                                                                                    tttACCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAggt      <  2:3804300/93‑1 (MQ=255)

                                                                                      tACCGCTGAGGCAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTGTGGTCCACACCATCGATATAggtgg    >  1:3629046/1‑93 (MQ=255)

                                                                                       aCCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAggtggt   <  1:602832/93‑1 (MQ=255)

                                                                                       aCCGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAggtggt   <  1:4985763/93‑1 (MQ=255)

                                                                                        ccGCTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGtt  <  2:5679943/93‑1 (MQ=255)

ACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTT  >  NC_000913/3982151‑3982331

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