BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,823,951	T→C	V372A (GTG→GCG)	recG →	DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,823,951	0	T	C	99.7%	995.1 / NA	286	V372A (GTG→GCG)	recG	DNA helicase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (124/161); minor base G (1/0); total (125/161)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 299 total aligned reads displayed.

CACCGCTCGGTATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTT  >  NC_000913/3823861‑3824040

 aCCGCTCGGTATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCa                                                                                        >  2:145142/1‑93 (MQ=255)

 aCCGCTCGGTATCGAAGTGGGCCGGCTCGCCGATAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCa                                                                                        <  2:4202751/93‑1 (MQ=255)

  ccGCTCGGTATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAg                                                                                       >  2:2794693/1‑93 (MQ=255)

   cGCTCGGTATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAga                                                                                      <  2:1252074/93‑1 (MQ=255)

    gCTCGGTATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAgat                                                                                     >  1:1099188/1‑93 (MQ=255)

    gCTCGGTATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAgat                                                                                     <  2:2187420/93‑1 (MQ=255)

    gCTCGGTATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAgat                                                                                     <  1:5433175/93‑1 (MQ=255)

    gCTCGGTATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAgat                                                                                     <  1:3611292/93‑1 (MQ=255)

    gCTCGGTATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAgat                                                                                     >  2:6401335/1‑93 (MQ=255)

     cTCGGTATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAgctg                                                                                    >  1:1867479/1‑93 (MQ=255)

     cTCGGTATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAgatg                                                                                    <  1:4066297/93‑1 (MQ=255)

      tCGGTATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAgatga                                                                                   <  1:2861457/93‑1 (MQ=255)

       cGGTATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAgatgat                                                                                  >  1:6316481/1‑93 (MQ=255)

        ggTATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGAtt                                                                                 <  2:3787240/93‑1 (MQ=255)

        ggTATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGAtt                                                                                 <  1:441849/93‑1 (MQ=255)

        ggTATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGAtt                                                                                 <  1:4209985/93‑1 (MQ=255)

         gTATCGAAGTGGGCTGGCGCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTg                                                                                <  1:1599001/93‑1 (MQ=255)

          tATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGt                                                                               <  1:3237013/93‑1 (MQ=255)

          tATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGt                                                                               >  1:1329030/1‑93 (MQ=255)

           aTCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGGGGGCAGGCGCAGATGATTGtg                                                                              >  2:5356647/1‑92 (MQ=255)

           aTCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTc                                                                              <  2:555978/93‑1 (MQ=255)

           aTCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTc                                                                              <  2:5063578/93‑1 (MQ=255)

           aTCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTCGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTc                                                                              <  2:4173591/93‑1 (MQ=255)

              gAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGt                                                                           <  2:4326901/93‑1 (MQ=255)

              gAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGt                                                                           >  2:5391524/1‑93 (MQ=255)

              gAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGt                                                                           >  1:2360296/1‑93 (MQ=255)

               aaGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTa                                                                          <  2:200651/93‑1 (MQ=255)

               aaGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTa                                                                          >  2:1006415/1‑93 (MQ=255)

                aGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTac                                                                         <  1:6377931/93‑1 (MQ=255)

                 gTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCCCGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTAGGTacc                                                                        >  2:4254498/1‑92 (MQ=255)

                 gTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTaca                                                                        >  2:685383/1‑93 (MQ=255)

                  tGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTacac                                                                       <  1:6185672/93‑1 (MQ=255)

                  ggggCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTacac                                                                       <  2:775040/92‑1 (MQ=255)

                    ggctggctCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTacacac                                                                     <  2:4029417/93‑1 (MQ=255)

                    ggctggctCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTacacac                                                                     <  2:3330161/93‑1 (MQ=255)

                     gctggctCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACg                                                                    >  1:5872115/1‑93 (MQ=255)

                      ctggctCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGc                                                                   <  2:350581/93‑1 (MQ=255)

                      ctggctCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGc                                                                   <  1:116011/93‑1 (MQ=255)

                      ctggctCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGc                                                                   <  2:2704069/93‑1 (MQ=255)

                        ggctCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCa                                                                 <  2:2069183/93‑1 (MQ=255)

                         gctCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCAt                                                                >  1:4899106/1‑93 (MQ=255)

                         gctCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCAt                                                                <  1:5622712/93‑1 (MQ=255)

                         gctCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCAt                                                                <  1:6115515/93‑1 (MQ=255)

                          ctCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATc                                                               <  2:6183339/93‑1 (MQ=255)

                            cGCCGGTAATCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCtt                                                             <  1:849062/93‑1 (MQ=255)

                             gCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTc                                                            >  1:3858701/1‑93 (MQ=255)

                               cGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCCTCCa                                                          >  1:1006416/1‑93 (MQ=255)

                                ggTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAg                                                         <  1:4238866/93‑1 (MQ=255)

                                ggTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAg                                                         <  2:4822045/93‑1 (MQ=255)

                                  tAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGa                                                       >  2:3653818/1‑93 (MQ=255)

                                  tAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGa                                                       <  2:188103/93‑1 (MQ=255)

                                    aGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAAc                                                     >  2:857018/1‑93 (MQ=255)

                                    aGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAAc                                                     <  1:4695047/93‑1 (MQ=255)

                                      cAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAg                                                   <  1:780822/93‑1 (MQ=255)

                                       aGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAgg                                                  >  2:4596216/1‑93 (MQ=255)

                                       aGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAgg                                                  >  1:298489/1‑93 (MQ=255)

                                        gAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAggg                                                 >  2:3758609/1‑92 (MQ=255)

                                        gAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGt                                                 >  1:3415308/1‑93 (MQ=255)

                                          aaGGTAAAGCACGGCTGACACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGc                                               <  1:5545454/93‑1 (MQ=255)

                                           aGGTAAAGCACGGCTGGCACAGCAGGAATCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCa                                              <  2:2588806/93‑1 (MQ=255)

                                            ggTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGCGCAg                                             >  1:2196606/1‑93 (MQ=255)

                                              tAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGtt                                           <  1:2250953/93‑1 (MQ=255)

                                                  gCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAAc                                       <  2:6257540/93‑1 (MQ=255)

                                                  gCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAAc                                       <  2:398166/93‑1 (MQ=255)

                                                      ggCTGGCACAGCAGGAAGCCATCGCCAGGGGTCAGGGGCAGATGATTGTCGGTACACACGCCACCTTCCAGGAACAGGTGCAGTTTAACGGcc                                   >  1:1692304/1‑93 (MQ=255)

                                                      ggCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGcc                                   <  2:6085159/93‑1 (MQ=255)

                                                      ggCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGcc                                   <  2:3861158/93‑1 (MQ=255)

                                                       gCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCt                                  >  2:5182910/1‑93 (MQ=255)

                                                       gCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGAGCAGGTTAACGGCCt                                  >  1:2427001/1‑93 (MQ=255)

                                                            cacaGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGgcgc                             <  2:467853/93‑1 (MQ=255)

                                                             acaGCAGGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCt                            >  1:3557235/1‑93 (MQ=255)

                                                              cagcagGAAGGCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTg                           <  1:5735145/93‑1 (MQ=255)

                                                                gcagTAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGt                         <  1:4763813/93‑1 (MQ=255)

                                                                gcagGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGt                         <  1:2141041/93‑1 (MQ=255)

                                                                  agGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGa                       >  2:4537524/1‑93 (MQ=255)

                                                                  agGAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGa                       <  1:2701782/93‑1 (MQ=255)

                                                                   ggAAGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGAt                      <  2:2626338/93‑1 (MQ=255)

                                                                      aGCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTAt                   >  2:4952789/1‑93 (MQ=255)

                                                                      aGCAATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTAt                   <  2:602797/93‑1 (MQ=255)

                                                                       gCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATc                  <  1:3664539/93‑1 (MQ=255)

                                                                       gCCATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATc                  >  2:4358577/1‑93 (MQ=255)

                                                                        ccATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATcg                 <  2:3054387/93‑1 (MQ=255)

                                                                         cATCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATcga                <  2:22104/93‑1 (MQ=255)

                                                                          aTCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATcgac               <  2:1504739/93‑1 (MQ=255)

                                                                          aTCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATcgac               <  2:3708888/93‑1 (MQ=255)

                                                                           tCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATcgacg              >  1:2511740/1‑93 (MQ=255)

                                                                           tCGCCAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATcgacg              >  2:3596020/1‑93 (MQ=255)

                                                                              ccAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAAc           <  1:463373/93‑1 (MQ=255)

                                                                              ccAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAAc           >  1:923821/1‑93 (MQ=255)

                                                                              ccAGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAAc           >  2:1117425/1‑93 (MQ=255)

                                                                                aGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAg         >  1:528841/1‑93 (MQ=255)

                                                                                aGCGGTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCCACGCACAg         >  2:2936408/1‑93 (MQ=255)

                                                                                  cgatCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACGGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCa       <  2:6321790/90‑1 (MQ=255)

                                                                                   ggTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCAt      <  2:3181309/93‑1 (MQ=255)

                                                                                    gTCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATc     >  1:4972146/1‑93 (MQ=255)

                                                                                     tCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCg    >  1:1437917/1‑93 (MQ=255)

                                                                                     tCAGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCg    <  2:5807260/93‑1 (MQ=255)

                                                                                       aGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGtt  <  1:5800388/93‑1 (MQ=255)

                                                                                       aGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGtt  >  1:6020955/1‑93 (MQ=255)

                                                                                       aGGCGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGtt  <  2:1411094/93‑1 (MQ=255)

CACCGCTCGGTATCGAAGTGGGCTGGCTCGCCGGTAAGCAGAAAGGTAAAGCACGGCTGGCACAGCAGGAAGCCATCGCCAGCGGTCAGGTGCAGATGATTGTCGGTACACACGCCATCTTCCAGGAACAGGTGCAGTTTAACGGCCTGGCGCTGGTGATTATCGACGAACAGCATCGTT  >  NC_000913/3823861‑3824040

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